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Research  on  convective  heat  transfer  on  internal  microtube  and  microchannel  has  been  extensively 
conducted  in  the  past  decade.  This  review  summarizes  numerous  researches  on  two  topics;  the  first 
section  focuses  on  studying  the  fluid  flow  and  heat  transfer  behavior  of  different  types  of  microtubes 
(MT)  and  microchannel  (MC)  at  different  orientations.  The  second  section  concentrates  on  nanofluids;  its 
preparation,  properties,  behavior,  and  many  others.  The  purpose  of  this  article  is  to  get  a  clear  view  and 
detailed  summary  of  the  influence  of  several  parameters  such  as  the  geometrical  specifications, 
boundary  conditions,  and  type  of  fluids.  The  maximum  Nusselt  number  is  the  main  target  of  such 
research  where  correlation  equations  were  developed  in  experimental  and  numerical  studies  are 
reported.  The  heat  transfer  enhancement  of  nanofluids  along  with  the  nanofluids  preparation  technique, 
types  and  shapes  of  nanoparticles,  base  fluids  and  additives,  transport  mechanisms,  and  stability  of  the 
suspension  are  also  discussed. 
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1.  Introduction 

The  characteristics  of  flow  and  heat  transfer  in  microchannels 
and  microtubes  have  attracted  much  attention  of  researchers 
because  of  the  rapid  developments  of  micro-electromechanical 
systems  (MEMS)  and  micrototal  analysis  system.  These  develop¬ 
ments  have  great  impacts  on  the  microelectronic  cooling  techni¬ 
ques,  the  microheat  exchanger,  bioengineering,  human  genome 
project,  medicine  engineering,  etc.  It  is  evident  that  the  under¬ 
standing  of  the  microscale  phenomena  is  very  important  for 
designing  efficient  micro-devices.  To  explore  the  fundamental 
physical  mechanisms  of  fluid  flow  and  heat  transfer  in  microchan¬ 
nels,  many  effects,  including  the  size  effect,  rarefaction  effect, 
surface  roughness,  viscous  effect,  electrostatic  force  effect,  axial 
heat  conduction  in  the  channel  wall,  surface  geometry,  the 
measurement  errors,  etc.  should  be  taken  into  account  [1].  A  large 
number  of  experimental  and  numerical  studies  focus  on  flow  and 
heat  transfer  behavior  in  microtube  and  microchannel  have  been 
reported.  These  vast  results  and  information  have  dealt  with 
different  conditions,  parameters,  geometry  dimensions,  and  instru¬ 
mentation,  which  indeed  used  for  comparison  purposes  to  indicate 
more  accurate  methodology  for  solving  the  case  studies. 

The  purpose  of  this  paper  is  threefold.  The  primary  purpose 
is  to  review  the  available  results  of  the  flow  and  heat  transfer 
presented  in  previously  published  data  for  microtubes  and  micro¬ 
channels.  The  results  of  massive  number  of  different  parameters 
affect  the  fluid  flow  and  heat  transfer  characteristics  are  summar¬ 
ized  and  presented  to  predict  the  peaks  of  the  flow  and  heat 
transfer.  The  second  purpose  is  to  understand  the  characteristics 
and  functions  of  nanofluids,  to  expect  their  effects  and  heat 
transfer  enhancement  in  such  geometries.  The  third  purpose  is 
to  open  a  gate  for  new  researches  and  propose  suggestions 
that  could  lead  to  improve  the  ability  to  enhance  the  heat  transfer 
in  MT  and  MC  using  nanofluids.  The  review  starts  with  an 
extensive  review  on  the  flow  and  heat  transfer  in  microtube  and 
microchannel.  Then  finish  with  a  comprehensive  review  for 
nanofluids  and  its  characteristics. 


2.  Microtube  (MT) 

2.1.  Experimental  studies  on  MT 

In  last  few  years,  large  number  of  researchers  intend  to  study  the 
heat  transfer  and  fluid  flow  in  MT.  Zhou  et  al.  [1]  investigated 
experimentally  and  numerically  the  flow  and  heat  transfer  char¬ 
acteristics  of  liquid  laminar  flow  in  MT.  They  used  the  smooth  fused 
silica  and  rough  stainless  steel  microtubes  with  the  hydraulic 
diameters  of  50-100  pm  and  373-1570  pm,  respectively.  Deionized 
water  was  used  as  the  working  fluid.  The  dimensions  and  surface 
roughness  for  the  fused  silica  tubes  and  stainless  steel  are  ranged 
50  pm,  75  pm  and  100  pm.  The  Reynolds  numbers  was  ranged  from 
20  to  2400.  The  friction  factors  for  smooth  fused  silica  tubes  are  in 
good  agreement  with  the  prediction  of  the  conventional  theory.  For 


the  rough  stainless  steel  tubes,  the  friction  factor  was  higher  than 
the  prediction  of  the  conventional  theory.  The  results  indicated  that 
the  conventional  friction  prediction  is  valid  for  water  flow  through 
MT  with  a  relative  surface  roughness  less  than  1.5%.  The  experi¬ 
mental  results  show  that  the  Nusselt  number  along  the  axial 
direction  do  not  accord  with  the  conventional  results  especially 
when  the  Reynolds  number  is  low  and  the  relative  tube  wall 
thickness  is  high. 

Lelea  et  al.  [2]  investigated  experimentally  and  numerically  the 
heat  transfer  and  fluid  flow  characteristics  in  a  microchannel.  The 
diameters  of  the  microchannel  used  were  0.1  mm,  0.3  mm  and 
0.5  mm,  the  flow  regime  was  laminar  and  the  Reynolds  number 
range  was  up  to  800.  The  tube  material  was  stainless  steel  SUS 
304.  The  distilled  water  was  used  as  a  working  fluid.  They  found 
that  the  results  confirm  that  the  conventional  or  classical  theories 
are  applicable  for  water  flow  through  microchannel  of  above  sizes. 

Ding  et  al.  [3]  investigated  experimentally  the  heat  transfer 
behavior  of  aqueous  suspensions  of  multi-walled  carbon  nano¬ 
tubes  (CNT  nanofluids)  flowing  through  a  horizontal  tube.  A 
straight  copper  tube  with  970  mm  length,  4.5  mm  inner  diameter, 
and  6.35  mm  outer  diameter  was  used  as  the  test  section.  Distilled 
water  and  multi-walled  carbon  nanotubes  were  used  to  produce 
nanofluids.  Fig.  1  shows  a  SEM  image  of  the  dispersed  sample. 
Significant  enhancement  was  observed  of  the  convective  heat 
transfer  in  compared  with  pure  water  as  the  working  fluid.  They 
found  this  enhancement  depends  on  the  flow  condition,  CNT 
concentration  and  the  pH  level.  The  axial  position  of  the  maximum 
enhancement  increases  with  the  CNT  concentration  and  the 
Reynolds  number.  The  results  indicate  that  the  Reynolds  number 
causes  a  significant  increase  in  the  convective  heat  transfer 
coefficient. 

Yang  and  Lin  [4]  investigated  experimentally  the  forced  con¬ 
vective  heat  transfer  performance  of  water  flow  through  six  micro 
stainless  steel  tubes  with  inner  diameters  ranging  from  123  pm  to 
962  pm.  The  results  indicate  that  the  conventional  heat  transfer 
correlations  for  laminar  and  turbulent  flows  can  be  well  applied 
for  predicting  the  fully  developed  heat  transfer  performance  in  the 
MT.  The  transition  occurs  at  Reynolds  number  from  2300  to  3000. 
The  surface  temperature  fluctuation  in  transition  regime  is  much 
larger  than  that  in  laminar  and  turbulent  flow  regimes.  They 
observed  that  the  heat  loss  is  less  than  4%  and  1%  of  the  total  heat 
input  for  123  pm  and  962  pm  tubes,  respectively.  At  the  lowest 
Reynolds  number  in  the  123  pm  tube,  a  4%  heat  loss  cause  25% 
heat  transfer  coefficient  overestimate. 

Celata  et  al.  [5]  investigated  experimentally  the  influence  of 
channel  wall  roughness  and  the  channel  wall  hydrophobicity  on 
adiabatic  flow  in  circular  microchannels.  Different  diameters  of  MT 
ranges  from  70  pm  to  326  pm  were  used  with  Reynolds  number 
less  than  300  for  all  diameters.  The  experiments  were  carried  out 
on  sooth  glass/fused  silica  and  a  Teflon  capillary  tubes  used  with 
water.  The  results  indicate  that,  with  degassed  water,  there  was  no 
effect  of  slip  flow  noted  due  to  hydrophobic  channel  walls  even  at 
70  pm  inner  diameter.  For  roughened  glass  channels,  an  increased 
in  friction  factor  above  64/Re  was  observed  only  at  the  smallest 
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Nomenclature 

U,  V 

Dimensionless  velocities  in  x  and  y  directions 

A1203 

Aluminum  oxide 

U0 

V 

Average  jet  velocity  at  the  entrance  (m/s) 

Axial  velocity  (m/s) 

cP 

Specific  heat  of  the  fluid  (J/kg  I<) 

ZnO 

Zinc  oxide 

CuO 

df 

Copper  oxide 

Diameter  of  base  fluid  molecule 

Creek  symbols 

dp 

g 

Nanoparticle  diameter  (m) 

Gravitational  acceleration  (m/s2) 

Meff 

Effective  viscosity 

k 

Thermal  conductivity  (W/m  1<) 

a 

Thermal  diffusivity  (m2/s) 

^eff 

Effective  thermal  conductivity 

p 

Dynamic  viscosity  of  fluid  (kg/ms) 

k  sf 

wall-to-fluid  thermal  conductivity  ratio 

V 

Kinematic  viscosity  of  fluid  (m2/s) 

L 

Length  of  the  tube  (m) 

0 

Inclination  of  tilted  wall 

M 

Molecular  weight  of  base  fluid 

p 

Fluid  density  (kg/m3) 

N 

Avogadro  number  (N= 6.022  x  1023  mob1) 

pi 

Nanofluid  density  (kg/m3) 

Nu 

Nusselt  number  ( Nu  =  hD/K ) 

4> 

Volume  fraction  of  nanoparticles 

P 

Pressure  of  fluid  (Pa) 

8 

Non-dimensional  wall  thickness 

Pr 

Prandtl  number  (Pr=Cp^0/k) 

£ 

Relative  roughness  height 

<7x 

Re 

Fleat  flux  (W/m2) 

Reynolds  number  ( Re=p0V0Dlfi0 ) 

Subscript 

Rnp 

Si02 

Nanoparticle  radius  (m) 

Silicon  oxide 

bf 

nf 

np 

eff 

base-fluid 

T 

Temperature  of  fluid  (K) 

nanofluids 

T 

To 

Bulk  temperature  (K) 

Reference  temperature 

nanoparticle 

effective 

U',  V 

Velocities  in  x'  and  y'  directions  (m/s) 

diameter  of  126  pm.  The  higher  factor  was  caused  by  actual 
deformation  of  channel  circularity  rather  than  increased  friction 
at  the  rougher  wall,  and  similar  behavior  was  observed  in  a 
Teflon  tube. 

Peng  et  al.  [6]  investigated  water  flow  in  a  fused  glass  microtube 
with  inner  diameter  of  230  pm  and  the  outer  diameter  is  about 
500  gin.  The  Reynolds  number  was  ranged  from  1540  to  2960.  The 
results  indicates  that,  the  flow  transition  from  laminar  to  turbulent 
occurs  at  Reynolds  number  1700  to  1900,  and  the  turbulence 
becomes  fully  developed  at  Re  >  2500.  The  details  of  the  flow 
structure  such  as  stream  wise  streaks  and  transverse  vortices  were 
observed  in  the  transitional  flow  fields  for  1800  <Re<  2100. 
The  instantaneous  velocity  which  observed  by  micro-PIV  reveal  large 
scale  coherent  structures  in  the  turbulent  flow  in  MT. 


Fig.  1.  SEM  image  of  a  dispersed  sample  of  carbon  nanotubes  [3], 


Celata  et  al.  [7]  investigated  experimentally  the  characteristics  of 
the  adiabatic  behavior  of  single-phase  laminar  flow  in  circular  micro¬ 
tube.  The  diameters  were  used  from  528  down  to  120  gm.  The 
experiments  were  carried  out  using  water,  which  was  degassed  by 
passing  through  a  very  small  and  continuous  of  helium.  The  result 
indicates  that,  the  Nusselt  number  decreased  with  decreasing  the 
diameter.  The  smaller  MT  shows  progressively  less  evidence  of  flow  in 
the  course  of  thermal  development.  The  decreased  in  heat  transfer 
performance  with  respect  to  the  conventional  value  of  Nusselt  number 
(Nu= 4.36)  was  more  marked  for  decreasing  the  Reynolds  number. 

Sara  et  al.  [8]  investigated  experimentally  the  laminar  con¬ 
vective  mass  transfer  and  friction  factor  in  MT  with  diameter  of 
0.20  mm  and  L/D  in  the  range  of  100  to  500.  The  Reynolds  number 
was  ranged  from  40  to  1400.  The  distilled  water  was  used  as  a 
working  fluid.  The  results  reveal  that  the  friction  factor  was  in  a 
very  good  agreement  with  the  Hagen-Poiseuille  theory  in  the 
laminar  regime  and  the  conventional  correlations  can  be  used  to 
determine  the  friction  factor  for  MT  with  ID  >  0.20  mm.  For  mass 
transfer,  the  Sherwood  number  for  MT  was  smaller  than  those 
obtained  from  the  correlations  developed  for  macrotubes. 

Zhang  and  Fu  [9]  investigated  experimentally  the  characteris¬ 
tics  of  two-phase  flow  of  liquid  nitrogen  in  vertically  upward 
micro-tubes  with  inner  diameters  0.531  and  1.042  mm.  The  results 
indicate  that,  the  main  two-phase  flow  patterns  of  liquid  nitrogen 
in  MT  were  bubbly  flow,  slug  flow,  churn  flow  and  annular  flow. 
The  confined  bubble  flow  was  observed  in  0.531  mm  and 
1.042  mm  MT.  The  dry-out  easily  happened  in  0.531  mm  diameter 
tube  and  led  to  high  vapor  quality  and  mist  flow  appeared.  The 
flow  regime  transition  lines  moved  to  lower  superficial  vapor 
velocity  as  the  diameter  of  MT  decreased.  In  addition,  there  was 
quite  large  divergence  between  the  experimental  results  and  the 
existed  theoretical  models,  and  it  seemed  that  there  was  not  yet  a 
theoretical  model  which  can  predict  the  flow  regime  transition 
well  for  the  two-phase  flow.  The  slip  ratio  between  two  phases 
became  larger  as  the  inner  diameter  of  the  MT  decreased. 

Liu  et  al.  [10]  investigated  experimentally  the  forced  convective 
heat  transfer  characteristics  in  quartz  microtubes  with  inner 
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diameters  of  242  pm,  315  |rm  and  520  pm.  The  de-ionized  water 
was  used  as  the  working  fluid.  The  Reynolds  number  was  ranged 
from  100  to  7000.  The  results  indicate  that,  the  experimental 
Nusselt  number  is  a  good  agreement  with  that  of  the  laminar 
correlations  when  the  flow  state  was  laminar.  The  Reynolds 
number  of  the  laminar  to  turbulent  transition  region  may  be 
about  1500-4000,  1600-4500  and  1900-5500  for  MT  with  inner 
diameters  of  242  pm,  315  pm  and  520  pm,  respectively. 

Celata  et  al.  [11]  investigated  experimentally  the  effects  of 
compressibility  for  helium  flow  through  fused  silica  micro-tubes. 
The  MT  diameter  was  ranged  from  30  pm  to  254  pm.  Reynolds 
number  was  ranged  from  0.8  to  500.  Two  hydrodynamic  models 
for  compressible  flow  (isothermal  and  adiabatic)  were  confronted 
at  normal  conditions.  Local  friction  factor  was  devised  to  compare 
with  the  global  friction  factor  to  quantify  the  effect  of  compres¬ 
sibility.  The  results  reveal  that,  the  friction  factor  has  remarkably 
accurate  prediction  by  the  reference  curve  of  Hagen-Poiseuille  and 
a  very  close  agreement  between  the  incompressible  and  quasi- 
compressible  correlations  was  obtained.  The  effect  of  density- 
change-induced  acceleration  was  extremely  limited  on  the  particle 
scale  and  that  an  incompressible  characterization  of  flow. 

Qi  et  al.  [12,13]  investigated  experimentally  the  flow  boiling  of 
liquid  nitrogen  in  the  MT  with  the  diameter  range  from  0.531  mm 
to  1.931  mm.  They  investigated  the  heat  transfer  characteristics 
and  the  critical  heat  flux  (CHF).  The  results  reveal  that,  the  wall 
increased  in  the  single-phase  flow  region,  drops  abruptly  at  ONB 
and  then  decreased  monotonically  in  the  two-phase  flow  region 
and  this  decreased  was  a  result  of  large  two-phase  flow  pressure 
drop.  For  low  mass  quality  region,  heat  transfer  was  governed  by 
nucleate  boiling  and  heat  transfer  coefficient  was  dependent  on 
heat  flux  and  pressure.  In  addition,  the  heat  transfer  coefficient 
increases  with  mass  quality  at  low  and  medium  heat  fluxes.  CHF 
increased  with  the  mass  flux  whereas  the  critical  mass  quality 
decreased  with  mass  flux.  The  CHF  was  governed  by  the  liquid  film 
near  the  inner  wall. 

Wen  and  Ding  [14]  investigated  experimentally  the  effects  of 
convective  heat  transfer  of  nanofluids  at  the  entrance  region  of  the 
copper  tube  with  970  mm  length.  The  de-ionized  water  with  A1203 
nanoparticles  were  used  as  working  fluids.  Reynolds  number  was 
ranged  from  500  to  2100.  The  results  indicate  that,  the  use  of  A1203 
nanoparticles  as  the  dispersed  phase  in  water  significantly 
enhanced  the  convective  heat  transfer  in  the  laminar  flow  regime 
and  this  enhancement  increased  with  Reynolds  number  and  the 
particle  concentration.  The  enhancement  was  particularly  signifi¬ 
cant  in  the  entrance  region  and  decreased  with  the  axial  distance. 
The  thermal  developing  length  of  nanofluids  was  greater  than  that 
of  pure  base  liquid  and  increased  with  an  increase  in  particle 
concentration. 

2.2.  Numerical  studies  on  MT 

Xiong  and  Chung  [15]  proposed  a  novel  bottom-up  approach  to 
generate  a  three-dimensional  random  roughness  on  a  MT  surface. 
They  used  Gaussian  random  number  generator  to  generate  the 
corner  points  in  the  space  with  random  heights  and  the  bi-cubic 
Coons  patch  to  create  the  boundary  curves  and  surface  based  on 
the  corner  points.  Many  tiny  generated  surfaces  were  used  to  form 
the  random  rough  of  MT.  Computational  fluid  dynamics  (CFD) 
solver  was  used  to  isolate  the  roughness  effect  and  to  solve  the 
three  dimensional  equations  for  the  water  flow  through  generated 
the  rough  MT  with  50  pm  diameter  and  100  pm  length.  They 
assumed  N  points  in  space,  and  it  is  connected  one  by  one  using 
straight  lines  to  form  the  bones  of  the  rough  surface  in  step  two.  At 
the  third  step,  they  created  the  curve  surface  and  in  the  last  step, 
they  created  two  more  surfaces  for  the  inlet  and  outlet  and  then 
they  used  it  to  form  a  volume  as  shown  in  Fig.  2.  Reynolds 


numbers  ranged  from  40  to  2000  and  the  viscous  dissipation  is 
not  considered.  They  observed  that  the  temperature  has  large 
values  at  the  peaks  and  small  values  at  the  valleys  of  the  MT 
surface.  The  roughness  has  a  negligible  on  the  averaged  Nusselt 
number  because  the  effects  of  peaks  and  valleys  counteract  each 
other  for  the  whole  MT  surface.  The  local  Nusselt  numbers  does 
not  have  a  large  fluctuation  away  from  the  theoretical  value  and  it 
randomly  distributes  along  the  MT.  For  the  hydraulic  diameter 
they  noticed,  the  improper  use  of  the  hydraulic  diameter  can 
result  in  increasing  or  decreasing  the  friction  factor  and  Nu 
number. 

Zhang  et  al.  [16]  investigated  numerically  the  effects  of  wall  axial 
heat  conduction  and  fluid  axial  conduction  for  simultaneously 
developing  the  laminar  flow  and  heat  transfer  in  the  wall  of  MT 
with  constant  outside  wall  temperature.  The  results  indicated  that 
the  heat  transfer  process  is  most  sensitive  to  wall-to-fluid  con¬ 
ductivity  ration  kSf,  and  in  condition  when  30  <RePr<  11,560  and 
ksf>  25;  But  when  the  conditions  were  used  30  <  RePr  <  11,560  and 
ksf<  25,  the  increasing  tube  thickness  and  the  decreasing  ksf  could 
make  the  inner  wall  surface  approaching  the  uniform  heat  flux 
condition.  The  maximum  Nusselt  number  they  obtained  is  about 
4.3.  It  turns  out  that  the  basic  function  of  the  wall  axial  heat 
conduction  is  to  unify  the  inner  wall  surface  heat  flux. 

Giulio  and  Paola  [17]  investigated  numerically  the  roughness 
effects  on  the  heat  transfer  and  pressure  loss  performances  of 
microscale  tubes  and  channels.  The  finite  element  was  used  based 
on  the  fractional  step  or  projection  approach.  The  surface  rough¬ 
ness  was  modeled  through  a  number  of  peaks  and  valleys  along 
the  ideal  smooth  surface.  Different  peak  shapes  and  distribution 
was  considered  to  arrive  to  the  real  roughness.  The  relative 
roughness  height  (e)  ranged  from  0.0%  to  5.3%  was  used  and  the 
geometrical  parameter  which  representative  of  the  tubes  is  the 
diameter  ranged  from  50  pm  to  150  pm.  The  results  indicated  that 
for  both  configurations,  a  significant  increase  in  Poiseuille  number 
was  detected.  While  the  influence  of  roughness  on  the  heat 
transfer  is  smaller  and  more  sensitive  to  the  geometry  of  both 
duct  and  roughness  obstruction  and  it  is  highly  dependent  on  the 
tube  shape. 

Kamali  and  Binesh  [18]  investigated  numerically  the  effects  of 
heat  transfer  for  multi-wall  carbon  nanotube  (MWCNT)  based 
nanofluids  in  a  straight  tube  under  constant  wall  heat  flux  condi¬ 
tion.  The  tube  length  was  914.4  mm  and  the  tube  diameter  was 
1.55  mm.  Both  water  base  fluid  and  aqueous  1%  CNT  nanofluid  were 
used  as  fluids  in  simulations.  The  single  phase  model  was  used  of 
two  dimensional  axisymmetric  steady,  forced  laminar  convection 
flows  of  nanofluids  inside  a  straight  circular  tube.  The  results 
indicate  that  the  heat  transfer  coefficient  is  dominated  by  the  wall 
region  due  to  non-Newtonian  behavior  of  CNT  nanofluid. 

Hong  and  Asako  [19]  investigated  the  heat  transfer  character¬ 
istics  of  gaseous  flows  in  a  MT  with  constant  heat  flux  whose  value 
is  positive  or  negative  on  2D  compressible  laminar  flow  for  no-slip 
regime.  The  tube  diameter  range  was  from  10  pm  to  100  pm  and  the 
length  to  diameter  ratio  was  200.  The  constant  heat  flux  range  was 
from  10~4W/m2  to  104W/m2.  The  numerical  methodology  was 
based  on  the  Arbitrary-Lagranging-Eulerian  (ALE)  method.  They 
selected  stagnation  pressure  in  such  a  way  that  the  Mach  number  at 
exit  ranged  from  0.1  to  0.7  and  the  outlet  pressure  was  fixed  at  the 
atmosphere  pressure.  The  results  indicated  that,  in  the  case  of  slow 
flow,  the  wall  and  bulk  temperatures  normalized  by  heat  fluxes 
increase  for  both  positive  and  negative  heat  flux.  A  correlation  for 
the  prediction  of  the  wall  temperature  of  the  gaseous  flow  in  the 
MT  was  proposed.  Supplementary  runs  with  slip  boundary  condi¬ 
tions  for  the  case  of  D=10pm  were  conducted  and  rarefaction 
effect  was  discussed.  With  increased  Ma  number,  the  compressi¬ 
bility  effects  are  more  dominant  and  the  rarefaction  effect  is 
relatively  insignificant  where  Kn  number  is  less  than  0.0096. 
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Fig.  2.  Four  steps  to  form  the  rough  microtube:  (a)  points  (b)  lines  (c)  surfaces  and  (d)  volume  [15], 


Bianco  et  al.  [20]  investigated  numerically  the  development  of 
laminar  forced  convection  flow  of  a  water-Al203  nanofluid  in  a 
circular  tube  with  a  constant  heat  flux  at  the  wall.  Single  and  two- 
phase  model  were  employed  with  either  constant  or  temperature- 
dependent  properties.  The  Re  number  was  ranged  from  250  to 
1050.  The  heat  flux  was  ranged  from  5000  W/m2  to  10,000  W/m2 
and,  0  (the  nanoparticle  volume  fraction)  was  ranged  from  1%  to 
4%  for  both  constant  and  temperature-dependent  properties.  In  all 
cases  the  size  for  the  spherical  particles  was  considered  to  equal 
100  nm.  The  results  indicate  that  the  inclusion  of  nanoparticles 
produced  a  considerable  increase  of  the  heat  transfer  with  respect 
to  that  of  the  base  liquid.  The  heat  transfer  enhancement  was 
increased  with  the  concentration  of  particle  volume  and  this 
increase  was  accompanied  by  increase  wall  shear  stress  values. 
The  results  of  single  and  two-phase  models  were  quite  similar,  the 
maximum  difference  detected  was  11%  for  the  average  heat 
transfer  coefficient  especially  for  0=4%. 

Xiao  et  al.  [21]  investigated  numerically  the  effects  of  gas  flow 
in  microtube.  The  results  indicate  that,  the  heat  transfer  effects 
associated  with  rarefied  flow  were  reduced  for  the  second-order 
model  and  the  effects  of  the  second  order  terms  was  most 
significant  at  the  upper  of  the  slip  regime.  For  the  airflow  at 
standard  conditions,  the  maximum  second-order  changes  of 
Nusselt  number  was  on  the  order  of  15%.  Nusselt  number  was 
increased  relative  to  the  values  of  no-slip  when  the  temperature 
jump  effect  was  small. 

Koo  and  ICleinstreuer  [22]  investigated  the  effects  of  viscous 
dissipation  on  the  temperature  field  and  the  friction  factor  in 
microtube  (MT)  and  microchannel  (MC).  Three  working  fluids  were 
used,  water,  methanol  and  iso-propanol.  The  results  indicate  that, 
the  viscous  dissipation  effect  on  the  friction  factor  was  increased  as 
the  system  size  decreases.  For  water  flow  in  a  tube  with  D  <  50  pm, 
viscous  dissipation  becomes  significant  and  should  be  considered  for 
all  experimental  and  computational  analyses.  In  addition,  the  aspect 
ratio  of  a  channel  plays  an  important  role  in  viscous  dissipation,  as 
the  aspect  ratio  deviates  from  unity,  the  viscous  dissipation  effects 
increase.  Ignoring  the  viscous  dissipation  effects  could  ultimately 
affect  the  friction  factor  measurements  for  flow  in  MT. 

Lelea  and  Cioabla  [23,24]  investigated  numerically  the  effects  of 
viscous  dissipation  on  heat  transfer  and  fluid  flow  in  micro-tubes. 


Three  different  fluids,  water  and  two  dielectric  fluids,  HFE-7600  and 
FC-70  with  temperature  dependent  fluid  properties  were  used.  The 
diameter  ratio  of  the  MT  was  D,/Do=(0.1/0.33)mm  with  a  tube 
length  L=100mm.  Two  different  heat  transfer  conditions,  cooling 
and  heating  and  three  different  Brinkman  number  0.01,  0.1  and 
0.5  was  used.  Non-dimensional  parameters  were  obtained,  for 
water  (Re  =  258-1684,  Pr  =  6.4-6.9),  HFE-7600  (Re  =  186-1130, 
Pr  =  29.3-32.2)  and  FC-70  (Re  =  7.2-43.2,  Pr  =  189.8-317.4). 
The  results  reveal  that,  the  friction  factor  and  Poiseuille  constant 
PD  are  affected  at  Br  =  0.5  for  water  and  start  with  Br  =  0.1  for  highly 
viscous  fluids  HFE-7600  and  FC-70.  The  Nusselt  number  for  the 
cooling  case  for  Br  =  0.5  was  about  three  times  higher  than  for 
heating  case,  regardless  the  fluid  type. 

Satapathy  [25]  investigated  forced  convection  for  gaseous  slip 
flow  in  an  infinite  microtube  with  mixed  boundary  conditions. 
The  results  indicate  that  for  both  the  hot  and  cold  wall  situation  in 
the  thermal  entrance  region,  the  local  bulk-mean  temperature 
increased  as  Peclet  number  or  Knudsen  number  increases  and 
vice-versa.  The  local  Nusselt  number  increased  with  the  increase 
in  Peclet  number  and  decreased  with  the  increase  in  Knudsen 
number.  The  fully-developed  Nusselt  number  decreased  with  the 
increase  in  Knudsen  number.  For  a  fixed  value  of  Knudsen  number, 
it  reaches  to  constant  value  for  all  Peclet  numbers. 

Wang  and  Wang  [26]  investigated  numerically  the  influences  of 
three-dimensional  wall  roughness  on  the  laminar  flow  in  micro¬ 
tube.  The  relative  roughness  was  ranged  from  0  to  0.05  and  the 
wave  number  from  0  to  30.  Reynolds  number  was  ranged  from 
1  to  500.  The  results  reveal  that  if  the  wave  number  of  rough 
function  or  Reynolds  number  decreased,  the  disturbed  areas  of 
flow  field  becomes  larger.  The  variation  of  disturbed  area  was 
independent  of  the  variation  of  relative  roughness  if  the  relative 
roughness  is  a  small  quantity.  In  addition,  the  dimensionless 
pressure  drop  increased  with  the  increases  of  relative  roughness, 
Reynolds  number  or  wave  number  of  wall  function.  The  influence 
of  wall  roughness  on  the  flow  pattern  was  different  from  those  on 
the  pressure  drop. 

Aziz  and  Niedbalski  [27]  investigated  numerically  the  first  and 
second  order  slip  flow  effects  on  the  thermal  development  of 
dilute  gas  flow  in  a  micro-tube  with  axial  conduction  and  viscous 
dissipation.  The  results  indicate  that  for  the  fully  developed 


Table  1 

Summary  of  experimental  and  numerical  studies  for  MT. 


Researchers  Diameter/flow  nature  Reynolds  Fluid/gas  Findings 

number 


Zhou  et  al.  [1] 

50-1570  pm/laminar 

Lelea  et  al.  [2] 

0.1 -0.5  mm/laminar 

Ding  et  al.  [3] 

6.35  mm/laminar 

Yang  and  Lin 

123-962  pm/laminar  and 

[4] 

turbulent 

Celata  et  al.  [5] 

70-326  pm/laminar 

Peng  et  al.  [6] 

230  pm/laminar  and 
turbulent 

Celata  et  al.  [7] 

50-528  pm/laminar  and 
turbulent 

Sara  et  al.  [8] 

0.20  mm/laminar 

Zhang  and  Fu 

0.531  and  1.042  mm 

[91 

Liu  et  al.  [10] 

242,  315  and  520  pm/ 
laminar  and  turbulent 

Celata  et  al. 

30-254  pm/laminar 

ini 

Qi  et  al.  [12,13] 

0.531-1.931  mm 

Wen  and  Ding 

4.5  mm/laminar 

[14] 

Xiong  and 

50  pm/laminar 

Chung  [15] 

Zhang  et  al. 

120-528  pm/laminar 

[16] 

Giulio  and 

50-150  pm/laminar 

Paola  [17] 

Kamali  and 

1.55  mm/laminar 

Binesh  [18] 

Hong  and 

10-100  pm/laminar 

Asako  [19] 

Bianco  et  al. 

0.01  m/laminar 

[20] 

Xiao  et  al.  [21] 

- 

Koo  and 

50  pm/laminar 

Kleinstreuer 

[22] 

Lelea  and 

0.1  mm/laminar 

Cioabla 

[23,24] 

Satapathy  [25] 

- 

Wang  and 

- 

Wang  [26] 

Aziz  and 

- 

Niedbalski 

[27] 

Salman  et  al. 

50  pm/laminar 

[28] 

20-2400 
Up  to  800 
700-1300 

100- 

10,000 

Less  than 

300 

1540- 

2960 

50-3138 

40-1400 


100-7000 

0.8-500 

500-2100 

40-2000 

10-2280 

100-1700 

600-1200 

14-1829 

250-1050 

20-2000 

7.2-1684 

1-500 


10-1500 


Deionized  water 
Distilled  water 
Distilled  water 
/0.5  wt%  CNT 
Water 


The  Nusselt  number  along  the  axial  direction  did  not  accord  with  the  conventional  results 

Results  confirm  that  the  conventional  or  classical  theories  are  applicable  for  water  flow  through  microchannel  of  above  sizes 
At  Re =80  Reynolds  number  causes  a  significant  increase  350%  in  the  convective  heat  transfer  coefficient 

Heat  loss  of  the  total  heat  input  depends  on  the  inner  diameter  and  the  Reynolds  number  values 


Demineralized  water 
Deionized  water 


No  effect  of  slip  flow  noted  due  to  hydrophobic  channel  walls  even  at  70  pm  inner  diameter.  For  roughened  glass  channels,  an  increased  in 
friction  factor  above  64/Re  was  observed  only  at  the  smallest  diameter  of  126  pm 
The  flow  transition  from  laminar  to  turbulent  occurs  at  Reynolds  number  1700-1900 


Demineralized  water  The  Nusselt  number  decreased  with  decreasing  the  diameter 


Distilled  water 
Nitrogen  liquid 


Friction  factor  was  in  a  very  good  agreement  with  the  Hagen-Poiseuille  theory.  The  Sherwood  number  for  MT  was  smaller  than  those 
obtained  from  the  correlations  developed  for  macrotubes 

There  was  not  yet  a  theoretical  model  which  can  predict  the  flow  regime  transition  well  for  the  two-phase  flow 


Deionized  water  Nusselt  number  is  rough  agreement  with  that  of  the  laminar  correlations 


Helium 

Hydrogen  liquid 
Al203/water 

Water 


Friction  factor  has  remarkably  accurate  prediction  by  Hagen-Poiseuille  and  effect  of  density-change-induced  acceleration  was  extremely 
limited  on  the  particle  scale  and  that  an  incompressible  characterization  of  flow 
Heat  transfer  coefficient  increased  with  mass  quality  at  low  and  medium  heat  fluxes 

Significant  enhancement  the  convective  heat  transfer  in  the  laminar  flow  regime  and  this  enhancement  increased  with  Reynolds  number  and 
the  particle  concentration 

The  hydraulic  diameter  can  result  in  increasing  or  decreasing  the  friction  factor  and  Nu  number 


Water  Heat  transfer  process  is  most  sensitive  to  wall-to-fluid  conductivity  ration  ksf. 

The  influence  of  roughness  on  the  heat  transfer  is  smaller  and  more  sensitive  to  the  geometry  of  both  duct  and  roughness  obstruction  and  it  is 
highly  dependent  on  the  tube  shape 

Water/aqueous  1%  Heat  transfer  coefficient  was  dominated  by  the  wall  region  due  to  non-Newtonian  behavior  of  CNT  nanofluid 
CNT 

Gas  A  correlation  for  the  prediction  of  the  wall  temperature  of  the  gaseous  flow  in  the  MT  was  proposed 


Al203/water  Maximum  difference  detected  was  11%  for  the  average  heat  transfer  coefficient  especially  for  0=4% 


Gas  The  maximum  second-order  changes  of  Nusselt  number  was  on  the  order  of  15% 

Water,  methanol  and  Viscous  dissipation  effects  could  ultimately  affect  the  friction  factor  measurements  for  flow  in  MT 
iso-propanol 

Water,  HFE-7600  and  Friction  factor  and  Poiseuille  constant  P0  affected  on  Nu  at  different  values  of  Br 
FC-70 


Gas  The  dimensionless  pressure  drop  increased  with  the  increases  of  relative  roughness,  Reynolds  number  or  wave  number  of  wall  function 

The  influence  of  wall  roughness  on  the  flow  pattern  was  different  from  those  on  the  pressure  drop 

Gas  As  the  viscous  dissipation  gathers  strength,  the  Nusselt  number  reaches  a  minimum  and  then  increases  again  as  the  effect  of  viscous 

dissipation  continues  to  prevail  over  the  cooling  effect 


A1203/EG 

CuO/EG 

Si02/EG 

ZnO/EG 


Si02-EG  nanofluid  has  the  highest  Nusselt  number,  followed  by  ZnO-EG,  CuO-EG,  A1203-EG,  and  lastly  pure  EG.  The  Nusselt  number  for  all 
cases  increases  with  the  volume  fraction  but  it  decreases  with  the  rise  in  the  diameter  of  nanoparticles.  In  all  configurations,  the  Nusselt 
number  increases  with  Reynolds  number 


B.H.  Salman  et  al.  /  Renewable  and  Sustainable  Energy  Reviews  28  (2013)  848-880  853 


854 


B.H.  Salman  et  al.  /  Renewable  and  Sustainable  Energy  Reviews  28  (2013)  848-880 


hydrodynamic  condition,  the  analytical  solution  for  the  velocity  with 
the  second-order  slip  model  predicts  higher  velocities  in  the  central 
region  of  the  flow  but  the  trend  was  reversed  near  the  wall  region.  In 
addition,  the  local  Nusselt  number  decreased  as  the  fluid  flows 
through  the  cooled  section  of  the  tube.  However,  as  the  viscous 
dissipation  gathers  strength,  the  Nusselt  number  reaches  a  minimum 
and  then  increases  again  as  the  effect  of  viscous  dissipation  continues 
to  prevail  over  the  cooling  effect.  Salman  et  al.  [28]  investigated 
numerically  laminar  convective  heat  transfer  in  a  microtube  (MT)  with 
50  pm  diameter  and  250  pm  length  with  constant  heat  flux  using 
finite  volume  method  (FVM).  Different  types  of  nanofluids  A1203,  CuO, 
Si02  and  ZnO,  with  different  nanoparticle  sizes  25,  45,  65  and  80nm, 
and  different  volume  fractions  ranged  from  1%  to  4%  using  ethylene 
glycol  as  base  fluids  were  used.  This  investigation  covers  Reynolds 
number  in  the  range  of  10  to  1500.  The  results  have  shown  that  Si02- 
EG  nanofluid  has  the  highest  Nusselt  number,  followed  by  ZnO-EG, 
CuO-EG,  A1203-EG,  and  lastly  pure  EG.  The  Nusselt  number  for  all 
cases  increases  with  the  volume  fraction  but  it  decreases  with  the  rise 
in  the  diameter  of  nanoparticles.  In  all  configurations,  the  Nusselt 
number  increases  with  Reynolds  number.  A  summary  for  the  experi¬ 
mental  and  numerical  studies  in  MT  is  shown  in  Table  1. 


3.  MicroChannel  (MC) 

3.1.  Experimental  studies  on  MC 

Jung  et  al.  [29]  measured  experimentally  the  heat  transfer  coeffi¬ 
cient  and  friction  factor  of  nanofluids  in  rectangular  microchannel 
(MC).  Aluminum  dioxide  (A1203)  with  diameter  of  170  nm  nanofluids 
with  various  particle  volume  fractions  of  0.6%,  1.2%  and  1.8%  were  used 
as  working  fluids.  The  Reynolds  number  was  tasted  in  the  range  of 


Square:  17  channels,  Dh  =  0.762  mm 


Rectangular:  20  channels,  Dh  =  0,424  mm _ 
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W  shape  insert:  19  channels,  Dh  =  0.732  mm _ 
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N  shape:  19  channels,  Dh  =  0.536  mm 

Fig.  4.  Tube  shapes  [31], 


5  and  300.  The  results  indicated  that  the  convective  heat  transfer 
coefficient  of  A1203  nanofluid  in  laminar  flow  regime  was  increased  up 
to  32%  compared  to  the  distilled  water  at  a  volume  fraction  of  1.8% 
without  major  friction  loss.  The  Nusselt  number  increases  with 
increasing  the  Reynolds  number  in  laminar  flow  regime.  The  Nusselt 
number  was  successfully  correlated  with  Reynolds  number  and 
Prandtl  number  based  on  the  thermal  conductivity  of  nanofluids. 

Cheng  et  al.  [30]  reported  a  novel  concept  to  control  the 
interface  location  of  the  pressure  driven  multi-phase  flow  in 
microchannel.  An  H-shape  microfluidic  device  was  used  to 
demonstrate  the  concept  of  interface  control  as  shown  in  Fig.  3. 
Two  fluids,  aqueous  NaCl  solution  and  diluted  glycerol  were 
introduced  through  the  inlets  by  a  syringe  pump.  The  flow  rates 
for  NaCl  solution  and  diluted  glycerol  was  maintained  at  the  same 
value.  For  a  purely  pressure  driven  two-fluid  flow,  the  aqueous 
NaCl  solution  with  a  lower  viscosity  flows  faster  and  occupied  a 
smaller  width  of  the  microchannel.  The  results  show  that  this 
method  is  effective  for  low  flow  rates,  and  small  viscosity  ratio 
close  to  1.  The  experimental  results  demonstrate  a  new  method  to 
solve  the  unmatched  viscosity  problem  of  the  two-fluid  flow  in 
microchannel. 

Agarwal  et  al.  [31]  investigated  experimentally  the  heat  trans¬ 
fer  coefficients  in  six  non-circular  horizontal  MC  with  diameters 
ranged  from  0.424  mm  to  0.939  mm  of  different  shapes  during 
condensation  of  refrigerant  R134a  over  the  mass  flux  in  the  range 
of  150  kg/m2/s  to  750  kg/m2/s.  The  channels  included  barrel- 
shape,  N-shape,  rectangular,  square,  and  triangular  extruded  tubes 
as  shown  in  Fig.  4.  The  results  reveal  that  the  models  developed 
for  larger  tubes  result  in  significant  deviations  from  the  measured 
data,  because  the  models  were  not  intended  for  and  could  not 
account  the  flow  phenomena  and  interfacial  shear  behavior 
specific  to  MC.  For  square,  rectangular,  and  barrel-shaped  channels 
an  annular-film  flow  based  model  was  developed.  For  triangular, 
N-shaped,  and  W-insert  channels,  it  was  assumed  that  the  liquid 
fraction  was  retained  at  the  corners,  with  the  flow  at  the  rest  of  the 
wall  and  in  the  core  being  approximated  as  mist  flow.  The  annular- 
flow-based  model  recommended  to  be  used  for  the  square,  barrel¬ 
shaped  and  rectangular  channels,  while  the  mist-flow  correlation 
should  be  used  for  channels  with  sharp  corners. 

Roland  et  al.  [32]  investigated  experimentally  and  numerically  the 
hydrodynamics  in  two-dimensional  MC  having  700  pm  and  200  pm 
height.  Reynolds  number  used  was  ranged  from  200  to  8000.  The 
tested  fluids  were  demineralized  water  with  an  electrical  conductivity 
of  300  pS/cm.  In  addition,  two-dimensional  numerical  model  corre¬ 
sponding  to  the  symmetry  plane  of  the  test  section  was  developed. 
The  model  accounted  for  free  convection  heat  transfer  with  surround¬ 
ing  air  at  the  outside  walls  of  the  test  section.  The  results  indicate  that 
a  reduction  in  the  Nusselt  number  for  laminar  flow  at  the  microscale  is 
occurred.  In  the  laminar  regime,  the  validity  of  the  conventional 
transport  phenomena  theory  may  be  extended  to  the  small  scales 
(200  pm).  Transition  to  turbulence  was  observed  both  in  the  global 
pressure  losses  measurements  and  in  the  Nu  with  Re  evolutions.  In 
addition,  the  critical  Reynolds  numbers  were  found  to  vary  from  3500 
to  4500  which  are  consistent  with  the  values  observed  at  the  larger 
scales.  The  experimental  corrected  data  matched  both  with  the  Blasius 
hydrodynamic  law  and  the  Dittus  and  Boelter  heat  transfer  correlation. 

Wu  and  Cheng  [33]  investigated  experimentally  laminar  heat 
transfer  and  pressure  drop  of  water  through  13  different  trape¬ 
zoidal  silicon  micro-channels  with  different  geometry,  surface 
roughness  and  surface  hydrophilic  properties.  The  results  indicate 
that  the  bottom  to  top  width  ratio,  the  height  to  top  width  ratio 
and  the  length  to  diameter  ratio  have  a  great  effect  on  the  laminar 
Nusselt  number  and  apparent  friction  constant  of  the  trapezoidal 
MCs.  The  laminar  Nusselt  number  and  apparent  friction  constant 
of  the  trapezoidal  MCs  increased  with  the  increase  of  surface 
roughness  and  this  increase  is  more  obvious  at  large  Reynolds 
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numbers.  The  Nusselt  Number  and  apparent  friction  constant  are 
higher  when  using  strong  hydrophilic  surfaces  than  those  with 
weak  hydrophilic  surfaces  (silicon  surfaces).  In  addition,  the 
laminar  convective  heat  transfer  shows  two  different  character¬ 
istics  at  low  and  high  Reynolds  number  ranges.  For  very  low 
Reynolds  number  flow  (0<Re<100),  the  Nusselt  number 
increased  almost  linearly  with  the  rise  of  Reynolds  number. 
However,  the  increased  of  the  Nusselt  number  after  Re=100  was 
gentle  with  the  rice  of  Reynolds  number. 

Park  and  Punch  [34]  investigated  experimentally  the  friction 
factor  and  heat  transfer  in  multiple  MCs  with  hydraulic  diameters 
from  106  pm  to  307  pm  for  single  phase  liquid  flow.  The  Reynolds 
number  used  was  ranged  from  69  to  800.  Deionized  water  was 
used  as  a  working  fluid.  The  results  indicate  that  the  friction 
factors  obtained  by  experiments  on  the  MC  prove  that  the 
conventional  theory  for  fully-developed  flow  was  applicable 
within  the  range  of  the  values  of  parameters  they  used.  In  order 
to  predict  heat  transfer  rate  accurately,  they  proposed  an  empirical 
correlation  in  terms  of  Nu/(Re0  62Pr033) and  Brinkman  number 
confined  to  the  experimental  range.  The  correlation  was  expected 
to  be  useful  to  design  the  MC  devices  related  to  heat  transfer. 

Quan  et  al.  [35]  investigated  experimentally  in  determination 
annular  condensation  heat  transfer  coefficient  of  steam  in  two  silicon 
microchannels  with  hydraulic  diameters  of  127  pm  and  173  pm  with 
trapezoidal  cross  sections.  The  air  was  used  as  the  coolant  for  the 
steam  flowing  in  the  MC.  The  results  indicate  that  the  annular 
condensation  heat  transfer  coefficient  in  MC  increased  with  mass  flux 
and  quality  and  decreased  with  the  hydraulic  diameter.  The  annular 
condensation  heat  transfer  coefficient  based  on  a  semi-analytical 
approach  of  turbulent  liquid  film  with  correlations  for  two-phase 
pressure  drop  and  void  fraction  are  valid  in  MC. 

Lee  et  al.  [36]  investigated  experimentally  and  numerically  the 
behavior  in  single-phase  flow  through  rectangular  microchannels  with 
width  range  from  194  pm  to  534  pm.  Deionized  water  was  used  as  a 
working  fluid.  Reynolds  number  was  ranged  from  300  to  3500.  The 
result  indicates  that  the  heat  transfer  coefficient  increased  with  the 
decrease  of  the  channel  size  depending  on  the  flow  rate.  The 
experimental  results  were  compared  against  conventional  correlations 
to  evaluate  their  applicability  in  predicting  MC  heat  transfer.  The  wide 
disparities  revealed  that  the  mismatch  in  the  boundaiy  and  inlet 
conditions  between  the  MC  experiments  and  the  conventional  corre¬ 
lations  precluded  their  use  for  predictions.  The  numerical  results 
agreed  with  the  experimental  data,  suggesting  that  such  approaches 
when  coupled  with  carefully  matched  entrance  and  boundaiy  condi¬ 
tions  can  be  employed  with  confidence  for  predicting  heat  transfer 
behavior  in  MC.  The  results  also  confirm  that  the  simplified  thin  wall 
analysis  can  be  used  as  a  computationally  economical  alternative  to  a 
full  3D  conjugate  analysis. 

Harirchian  and  Garimella  [37]  investigated  experimentally  the  local 
flow  boiling  heat  transfer  of  a  dielectric  fluid  (Fluorinert  FC-77)  in 
microchannel  heat  sinks  as  shown  in  Fig.  5.  The  mass  flux  used  was 
ranged  from  250  kg/m2/s  to  1600  kg/m2/s.  Seven  different  test  pieces 
made  from  silicon  and  consisting  of  parallel  MCs  with  nominal  width 
ranged  from  100  pm  to  5850  pm  with  a  nominal  depth  of  400  pm 
were  obtained.  The  results  indicate  that  for  MCs  of  width  400  pm  and 
greater,  the  heat  transfer  coefficients  corresponded  to  a  fixed  wall  heat 
flux  as  well  as  the  boiling  curves  was  independent  of  channel  size.  The 
heat  transfer  coefficients  and  boiling  curves  was  independent  of  mass 
flux  in  the  nucleate  boiling  region  for  a  fixed  channel  size,  but  affected 
by  mass  flux  as  convective  boiling  dominates.  In  addition,  a  strong 
dependence  of  pressure  drop  on  both  channel  size  and  mass  flux  was 
observed. 

Ergu  et  al.  [38]  investigated  experimentally  the  pressure  drop 
and  point  mass  transfer  in  a  rectangular  microchannel  with  a 
width  of  3.70  mm,  height  of  0.107  mm  and  length  of  35  mm. 
Distilled  water  was  used  as  a  working  fluid.  Reynolds  number 


used  was  ranged  from  100  to  845.  The  results  indicate  that  the 
friction  factor  values  were  slightly  higher  than  those  calculated  by 
theoretical  correlation  and  these  values  increased  with  increasing 
Reynolds  number.  The  classical  laminar  flow  equations  can  be  used 
for  the  calculation  of  the  friction  factor  in  MCs. 

Owhaib  et  al.  [39]  investigated  experimentally  saturated  flow 
boiling  in  MCs.  Heat  transfer  coefficients  were  measured  for 
saturated  boiling  of  working  fluid  (R  134a)  in  vertical  circular 
tubes  with  internal  diameters  ranged  from  0.862  mm  to  1.7  mm, 
with  a  uniform  heating  length  of  220  mm.  Heat  transfer  coeffi¬ 
cients  were  obtained  for  a  heat  flux  ranged  from  3  kW/m2  to 
34  kW/m2  and  a  mass  flux  was  ranged  from  50  l<W/m2/s  to 
400  kW/m2/s  and  vapor  qualities  up  to  0.6.  The  results  reveal  that 
the  saturated  flow  boiling  heat  transfer  coefficient  was  more  or 
less  independent  of  vapor  quality  and  mass  flux,  but  it  is  highly 
dependent  on  the  heat  flux  for  vapor  fraction  x  <  0.6.  The  con¬ 
tribution  of  forced  convection  heat  transfer  to  the  overall  transfer 
rate  was  small  and  that  heat  transfer  by  a  mechanism  similar  to 
that  in  nucleate  boiling  was  dominant.  In  addition,  flow  boiling 
heat  transfer  coefficients  was  higher  for  smaller  diameter  tubes. 
Also,  an  increase  in  the  system  pressure  improves  the  heat  transfer 
performance. 

Chiu  et  al.  [40]  investigated  experimentally  and  numerically 
the  heat  transfer  performance  and  characteristics  of  liquid  cooling 
heat  sink  microchnnels.  The  height  of  MCs  was  constant. 
The  aspect  ratio  was  set  from  1.67  to  14.29  and  the  porosity  was 
from  25%  to  85%.  The  imposed  pressure  drop  used  was  ranged 
from  490  Pa  and  2940  Pa.  The  result  indicate  that  the  overall 
averaged  Reynolds  number  in  MCs  decreased  with  both  increasing 
the  aspect  ratio  and  cross-sectional  porosity  under  same  pressure 
drop.  The  aspect  ratio  for  the  lowest  effective  thermal  resistance 
was  dependent  on  the  pressure  drop  between  the  inlet  and  the 
outlet  of  heatsink.  In  addition,  the  heat  transfer  enhancement 
due  to  raising  pressure  drop  was  more  significant  on  MCs  with 
high  aspect  ratio  hence  thinner  channels  are  more  advantageous 
under  high  pressure  drop.  The  local  Nusselt  number  in  a  channel 
decreased  with  increasing  the  aspect  ratio.  The  value  of  Nu 
monotonically  decreased  along  the  channel  for  high  aspect  ratio. 

Shen  et  al.  [41  ]  investigated  experimentally  the  single  phase  heat 
transfer  in  a  compact  heat  sink  consisting  of  2  rectangular  MCs  of 
300  pm  width  and  800  pm  depths.  The  relative  roughness  was  ranged 
from  4%  to  6%.  Deinonized  water  was  used  as  the  working  fluid.  The 
Reynolds  number  used  was  ranged  from  162  to  1257.  The  inlet  liquid 
temperature  used  was  ranged  from  30  C  to  70  °C  and  the  heating 
powers  were  used  from  140  W  to  450  W.  The  results  indicate  that  the 
surface  roughness  has  a  great  effect  on  the  laminar  flow  in  rough  MCs. 
For  developed  flow,  the  Poiseuille  number  in  the  regime  of  high 
Reynolds  numbers  was  higher  than  the  conventional  theory  predic¬ 
tions  and  increase  with  increasing  the  Reynolds  number,  rather  than 
remaining  constant.  In  addition,  the  average  Nusselt  number  increased 
with  the  increase  of  Reynolds  number  and  Prandtl  number.  Both  the 
local  and  the  average  Nusselt  numbers  were  significantly  lower  than 
the  conventional  theory  predictions  for  all  flow  rates  and  this  was 
attributed  to  two  main  factors,  the  cross-sectional  aspect  ratio  and  the 
surface  roughness. 

Barber  et  al.  [42]  investigated  experimentally  the  effects  of 
hydrodynamics  and  heat  transfer  during  flow  boiling  instabilities 
in  a  single  MC  with  a  hydraulic  diameter  of  727  pm.  The  uniform 
heat  flux  used  was  4.26  kW/m2  and  the  inlet  liquid  mass  flow  rate 
held  constant  at  1.13  x  10~5  kg/s.  The  results  reveal  that  the  flow, 
pressure  and  temperature  conditions  leading  to  nucleate  boiling, 
flow  instabilities  and  transition  regimes  during  flow  boiling  in  a 
MC.  A  quantify  and  characterize  the  timescales  of  various  observed 
instabilities  during  flow  boiling  in  a  MC.  High  speed  imaging 
revealed  some  of  the  controlling  physical  mechanisms  responsible 
for  the  observed  instabilities. 
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Fig.  5.  (a)  MicroChannel  test  section  and  (b)  integrated  heaters  and  temperature  sensors  in  the  microchannel  test  piece  [37], 


Lee  and  Garimella  [43]  investigated  experimentally  the  flow 
boiling  in  arrays  of  parallel  MCs  by  using  a  silicon  test  piece  with 
imbedded  discrete  heat  sources  and  integrated  local  temperature 
sensors.  The  MCs  width  used  was  ranged  from  102  pm  to  997  pm 
and  the  width  was  400  pm.  Deionized  water  was  used  as  a 
working  fluid.  The  results  reveal  that  the  pressure  drop  across 
the  MCs  increased  rapidly  with  the  heat  flux  when  the  incipience 
heat  flux  (for  the  onset  of  nucleate  boiling)  was  exceeded.  At  low 
to  medium  heat  fluxes,  the  local  heat  coefficient  increased  almost 
linearly  with  heat  flux.  While,  at  higher  heat  fluxes,  the  saturated 
heat  transfer  coefficient  becomes  largely  insensitive  to  heat  flux 
for  the  range  tested. 

Celata  et  al.  [44]  investigated  experimentally  heat  transfer 
characteristics  of  subcooled  flow  boiling  of  FC-72  in  a  single 
horizontal  circular  cross  section  MC  with  480  pm  inner  diameter, 
800  pm  outer  diameter  and  102  mm  long.  Different  flow  patterns, 
both  in  the  stable  and  unstable  flow  boiling  regimes  were  con¬ 
sidered.  The  results  indicate  that  the  stable  flow  boiling  with 
alternating  bubbly/slug  flow,  slug/annular  flow  and  annular/mist 


flow  was  observed  for  heat  flux  of  150  kW/m2  or  higher  and  mass 
flux  of  1500  kg/m2  s  or  higher.  Back  and  forth  oscillations  with  flow 
instabilities  were  observed  in  cases  of  lower  heat  and  mass  fluxes. 
No  complete  reverse  flow  was  observed  in  upstream  direction.  In 
addition,  the  heat  transfer  coefficient  increased  with  the  increases 
of  vapor  quality  at  higher  vapor  quality,  and  at  lower  vapor  quality, 
the  heat  transfer  coefficient  was  rather  independent  of  heat  flux. 

Megahed  [45]  investigated  experimentally  flow  boiling  charac¬ 
teristics  in  a  cross-linked  MC  heat  sink  at  low  mass  fluxes  and  high 
heat  fluxes.  The  heat  sink  used  consists  of  45  straight  MCs  with  a 
hydraulic  diameter  of  248  pm  and  heated  length  of  16  mm.  The 
dielectric  coolant  FC-72  was  used  as  a  working  fluid.  The  heat  flux 
used  was  ranged  from  7.2  kW/m2  tol04.2  kW/m2,  mass  flux  from 
99  kg/m2/s  to  290  kg/m2/s,  and  exit  quality  from  0.01  to  0.71.  The 
results  indicate  that  the  instability  increased  with  decreasing  the 
mass  flux.  The  inlet  pressure  and  outlet  saturation  temperature 
oscillate  with  higher  amplitudes  as  the  mass  flux  decreases.  The 
two-phase  pressure  drop  through  MC  heat  sinks  with  cross-links 
was  much  higher  compared  to  the  straight  MC  heat  sinks  due  to  the 
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cross-link  effect.  The  two-phase  heat  transfer  coefficient  was 
strongly  increased  with  the  exit  quality  because  of  the  domination 
of  the  nucleation  boiling  mechanism  in  the  cross-link  region. 

Diaz  and  Schmidt  [46]  investigated  experimentally  flow  boiling 
heat  transfer  in  a  single  rectangular  MC  made  of  the  nickel  alloy 
Inconel  600  with  0.3  mm  height,  12.7  mm  width  and  200  mm 
length.  Water  and  ethanol  were  used  as  test  fluids.  The  mass  fluxes 
used  ranged  between  50  kg/m2/s  and  500  kg/m2/s  and  heat  fluxes 
up  to  400  l<W/m2  at  an  outlet  pressure  of  0.1  MPa.  The  results 
reveal  that,  high  amplitude  oscillations  for  water  were  found  in 
the  region  of  subcooled  boiling  at  qualities  slightly  below  zero.  The 
amplitude  decreased  as  quality  increases  and  no  characteristic 
frequency  was  observed.  Flow  boiling  of  ethanol  yielded  low 
amplitude  oscillations  of  the  wall  temperature,  even  within  the 
subcooled  boiling  region.  Increase  in  quality  resulted  in  increasing 
oscillation  amplitude.  In  addition,  the  heat  transfer  coefficient 
decreased  with  increasing  quality  for  water.  Different  observed  for 
ethanol  after  the  low  quality  region.  The  heat  transfer  coefficient 
increased  with  increasing  quality  at  low  heat  fluxes.  With  increas¬ 
ing  heat  flux,  the  heat  transfer  coefficient  decreased  with  quality. 

Chen  et  al.  [47]  investigated  experimentally  the  condensation 
of  steam  in  a  series  of  triangular  silicon  MCs  with  hydraulic 
diameters  ranged  from  100  to  250  pm.  The  number  of  channels 
in  each  silicon  chip  was  10  and  the  length  for  each  channel  was 
56.7  mm.  The  results  indicate  that  the  droplet,  annular,  injection 
and  slug-bubbly  flow  are  the  dominant  flow  patterns  in  triangular 
silicon  MCs.  The  flow  of  condensation  instability  was  higher  with 
increased  inlet  vapor  Reynolds  number,  condensate  Weber 
number  and  the  prolongation  of  the  injection  location,  or  with  a 
decrease  in  the  hydraulic  diameter  of  the  channel.  In  addition,  the 
wall  temperature  of  the  channel  decreased  along  the  condensation 
stream.  The  total  pressure  drop,  the  average  condensation  heat 
transfer  coefficient  and  the  average  Nusselt  number  were 
observed  to  be  larger  with  increased  inlet  vapor  Reynolds  number. 
Moreover,  it  was  found  that  the  condensation  heat  transfer  was 
enhanced  by  a  reduction  in  the  channel  scale. 

Betz  and  Attinger  [48]  investigated  experimentally  heat  transfer 
enhancement  in  MC  heat  sink  with  the  introduction  of  segment 
flow.  The  segmented  flow  pattern  was  created  by  the  periodic 
injection  of  air  bubbles  into  water-filled  channels.  Reynolds  number 
used  was  ranged  from  160  to  1580.  The  results  reveal  that  the 
segmented  flow  would  provide  an  intermediate  step  between 
single-phase  and  boiling  flow  for  the  purpose  of  electronic  cooling. 
The  segmented  flow  increases  the  Nusselt  number  of  laminar  flow 
by  more  than  100%.  The  heat  transfer  enhancement  only  occurs  for 
a  specific  range  of  flow  rates  and  capillary  numbers.  At  lower  or 
higher  capillary  numbers,  there  was  no  significant  heat  transfer 
enhancement  observed  because  segmented  flow  was  replaced  by 
bubbly  or  churn  flow,  respectively. 

Wojtan  et  al.  [49]  investigated  experimentally  the  saturated 
critical  heat  flux  (CHF)  in  0.5  mm  and  0.8  mm  internal  diameter 
MC  tubes  as  a  function  of  refrigerant  mass  velocity,  heated  length, 
saturation  temperature  and  inlet  liquid  subcooling.  The  tested 
refrigerants  were  R-134a  and  R-245fa  and  the  heated  length  of  MC 
was  varied  between  20  mm  and  70  mm.  The  results  reveal  that  a 
strong  dependence  of  CHF  on  mass  velocity,  heated  length  and  MC 
diameter  but  no  influence  of  liquid  subcooling  (2-15  °C)  was 
observed.  A  new  MC  version  of  the  Katto-Ohno  correlation  was 
developed  to  predict  the  CHF  in  circular,  uniformly  heated  MCs. 

Asthana  et  al.  [50]  investigated  experimentally  convective  heat 
transfer  in  serpentine  MCs  cross  section  100  pm  with  a  segmented 
flow  of  two  immiscible  liquids  for  flow  rates  up  to  130  pl/min.  The 
results  reveal  that  for  slug  flow,  the  movement  of  the  oil  slugs 
induces  intensive  fluid  recirculation  in  the  carrier  water  phase 
between  the  slugs,  which  disrupts  the  formation  of  a  boundary 
layer  and  improves  mixing,  more  so  for  increasing  oil  volume 


fractions.  In  addition,  the  local  temperature  of  the  water  region 
separating  the  droplets  was  higher  than  in  the  case  of  pure  water 
flow.  For  segmented  flow,  up  to  a  four-fold  increased  the  Nusselt 
number  compared  to  pure  water  flow  was  observed.  Increased 
heat  removal  increased  the  pressure  drop,  which  was  strongly 
affected  by  the  fluid  temperature. 

Wang  and  Cheng  [51]  investigated  experimentally  stable  and 
unstable  flow  boiling  phenomena  of  deionized  water  in  a  single 
MC  with  hydraulic  diameter  of  155  pm.  The  results  reveal  that 
stable  flow  boiling  occurred  in  a  single  MC  if  exit  vapor  quality 
xe  <  0.013,  where  vapor  bubbles  generated  and  pushed  away  by 
the  incoming  subcooled  liquid.  Single-phase  liquid  flow,  bubbly 
flow  and  elongated  bubbly/slug  flow  were  observed  sequentially 
in  the  stable  flow  boiling  mode.  In  addition,  unstable  flow  boiling 
occurred  in  a  single  MC  if  xe  >  0.013,  where  vapor  plug  expanded 
upstream  resulting  in  oscillations  of  temperature  and  pressure.  In 
this  unstable  flow  boiling  mode,  bubbly  flow,  elongated  bubbly/ 
slug  flow,  semiannular  flow  and  annular/mist  flow  were  observed 
at  different  times  near  the  outlet  depending  on  the  mass  flux. 

Wu  et  al.  [52]  investigated  experimentally  steam  condensation 
flow  in  wide  rectangular  silicon  MCs  with  the  hydraulic  diameter 
of  90.6  pm,  483.4  pm  width  and  50  pm  depth.  The  results  reveal 
that  droplet-annular  compound  flow,  injection  flow,  and  vapor 
slug-bubbly  flow  were  observed  along  the  channel,  which  differ 
from  that  in  other  cross-sectional  shape  MCs.  In  the  droplet- 
annular  compound  flow  region,  the  vertical  walls  (short  side)  of 
the  channel  were  completely  covered  by  the  condensate,  while 
droplet  condensation  still  exists  on  the  horizontal  wall  (long  side) 
of  the  channel.  The  location  of  the  injection  flow  postponed  with 
the  increasing  inlet  vapor  Reynolds  number.  The  injection 
frequency  increased  with  increasing  the  inlet  vapor  Reynolds 
number  and  condensate  Weber  number.  In  addition,  the  frequency 
in  the  wide  rectangular  MCs  was  lower  than  that  in  triangular  MCs 
with  the  same  hydraulic  diameter.  It  is  confirmed  that  the  cross- 
sectional  shape  of  the  MC  plays  a  significant  role  on  the  instability 
of  condensation  flow.  The  correlation  of  Nusselt  number  was  also 
obtained. 

Steinke  and  Kandlikar  [53]  investigated  experimentally  the 
different  components  of  the  total  pressure  drop  measurement 
for  adiabatic  and  single-phase  degassed  water  in  MC  with  a  width 
of  200  pm,  a  depth  of  250  pm  and  a  length  of  10  mm.  An  in-depth 
comparison  of  previous  experimental  data  was  performed  to 
identify  the  discrepancies  in  reported  literature.  The  Reynolds 
number  used  was  ranged  from  14  to  789.  The  results  reveal  that, 
the  procedure  for  correcting  the  measured  pressure  drop,  the  data 
was  corrected  for  inlet  and  exit  losses  and  then  it  was  corrected  for 
developing  flows.  The  corrected  data  shows  good  preliminary 
agreement  in  value  and  trend  with  the  conventional  theory  for 
laminar  fluid  flow.  The  conventional  Stokes  and  Poiseuille  flow 
theories  apply  for  single-phase  liquid  flow  in  MC  flows. 

Wu  et  al.  [54]  investigated  experimentally  static  and  dynamic 
flow  instability  of  a  parallel  MC  heat  sink  consisting  of  26 
rectangular  MCs  with  300  pm  width,  800  pm  depth,  at  high  heat 
fluxes.  Water  and  methanol  were  used  as  a  working  fluid  with 
mass  fluxes  of  20-1200  kg/m2/s.  The  inlet  temperature  used  was 
ranged  from  300  °C  to  70  °C.  The  results  indicate  that  the  onset  of 
flow  instability  occurs  at  the  outlet  fluid  temperature  of  93-96  °C, 
which  were  several  degrees  lower  than  the  saturated  temperature 
corresponding  to  the  exit  pressure  of  the  MCs.  In  addition,  the 
outlet  temperature  and  inlet  pressure  were  always  in  phase,  but 
they  were  out  of  phase  with  the  pressure  drop  across  the  MCs.  The 
wall  temperatures  generally  match  the  oscillation  trends  of  the 
inlet  pressure  and  the  outlet  temperature,  but  a  time  delay  occurs 
due  to  the  thermal  inertia  of  the  copper  block. 

Zhang  et  al.  [55]  investigated  experimentally  the  control  of 
pressure-drop  flow  instabilities  in  boiling  microchannel  systems  with 
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61  urn  width,  272  pm  height  and  15  mm  length.  Deionized  water  was 
used  as  a  working  fluid.  The  results  reveal  that  the  predictions  from 
the  flow  oscillation  model  were  agreed  well  with  the  experimental 
pressure-drop  observations  across  a  flow  meter  and  a  MCHS.  In 
addition,  a  family  of  state  and  dynamic  output-feedback  active  flow 
controllers  was  developed  and  evaluated  for  the  dynamic  pressure- 
drop  instability  suppression.  The  flow  oscillation  amplitude  can  be 
regulated  to  whatever  level  is  desired  under  certain  conditions. 

Kohl  et  al.  [56]  investigated  experimentally  the  flow  with 
internal  pressure  measurements  in  straight  microchannel  test 
section  with  hydraulic  diameters  ranged  from  25  pm  to  100  pm. 
The  Reynolds  number  used  was  ranged  from  4.9  to  18,814.  The 
results  suggest  that  friction  factors  for  MCs  can  be  accurately 
determined  from  data  for  standard  large  channels. 

Rosengarten  et  al.  [57]  investigated  experimentally  the  effect  of 
contract  angle  (surface  wettability)  on  the  convective  heat  transfer 
coefficient  in  microchannels  using  thin  film  coating  technology. 
Distilled  water  was  used  as  a  working  fluid.  The  results  reveal  that, 
the  pressure  drop  decreases  relative  to  that  expected  with  non-slip 
theory  above  a  critical  shear  rate  and  this  decrease  was  dependent 
on  the  contact  angle.  Associated  with  the  decrease  in  pressure 
drop  was  a  concurrent  decrease  in  the  convective  heat  transfer 
coefficient  which  again  was  more  pronounced  for  hydrophobic 
surfaces.  Slip  flow  was  an  extremely  important  for  micro  and  in 
particular  nano  heat  exchanger  designs  as  heat  transfer  rates  may 
decrease  more  than  10%  at  high  shear  rates. 

Tang  et  al.  [58]  investigated  experimentally  compressibility,  rough¬ 
ness  and  rarefaction  influences  for  nitrogen  and  helium  in  stainless 
steel  MTs  with  diameters  ranged  from  119  pm  to  300  pm,  fused  silica 
MTs  with  hydraulic  diameter  ranged  from  52  pm  to  201  pm  and  fused 
silica  square  MCs  with  diameters  ranged  from  52  pm  to  100  pm.  The 
Reynolds  number  used  was  ranged  from  3  to  6300.  The  tested  outlet 
Knudsen  number  used  was  ranged  from  1.4  x  10~4  to  4  x  10~3,  and 
the  tested  outlet  Mach  number  was  ranged  from  0.01  to  0.6.  The 
results  indicate  that  the  gaseous  flow  friction  factor  in  stainless  steel 
tube  (SST)  was  observed  to  be  much  higher  than  the  theoretical 
prediction  for  conventional  size  tubes.  The  friction  factor  in  the  fused 
silica  tubes  and  fused  silica  channels  (FST,  FSC)  agreed  well  with  the 
theoretical  prediction.  In  addition,  the  increase  in  friction  factor  was 
attributed  to  the  dense  roughness  distribution  in  the  inner  surface  of 
the  stainless  steel  tubes.  Surface  roughness  in  MCs  was  found  to  affect 
the  friction  factor.  It  is  suggested  that  for  gaseous  flow  in  MCs  with  a 
relative  surface  roughness  less  than  1%,  the  conventional  laminar 
prediction  should  still  be  applied. 

Huh  et  al.  [59]  investigated  experimentally  the  unique  character¬ 
istics  of  flow  boiling  in  a  polydimethylsiloxane  (PDMS)  rectangular 
single  microchannel  with  a  hydraulic  diameter  of  103.5  pm  and  length 
of  40  mm.  Deionized  water  was  used  as  a  working  fluid.  The  mass 
fluxes  used  were  170  kg/m2  and  360  kg/m2  and  the  heat  fluxes  were 
ranged  from  200  kW/m2  to  530  kW/m2.  The  results  indicate  that  flow 
boiling  fluctuations  with  a  very  long  period  and  large  amplitude  were 
observed.  Both  the  pressure  drop  and  mass  flux  had  these  fluctuations, 
but  a  phase  shift  occurred  between  them  that  controlled  the  flow 
patterns.  The  fluctuation  of  wall  temperature  and  heat  flux  had  the 
same  period  as  that  of  the  mass  flux.  In  addition,  the  variation  in 
temperatures,  pressure,  and  mass  flux  exactly  matched  the  flow 
pattern  transitions  between  bubbly/slug  flow  and  elongated  slug/ 
semiannular  flow.  The  flow  pattern  transition  instabilities  in  the  single 
MC  caused  the  cyclic  behavior  of  the  wall  temperature,  pressure  drop, 
and  mass  flux.  As  the  supplied  heat  flux  and  vapor  quality  increased, 
the  period  of  the  unstable  fluctuations  lengthened  and  the  amplitude 
also  increased. 

Mokrani  et  al.  [60]  investigated  experimentally  the  fluid  flow 
and  convective  heat  transfer  in  flat  microchannels  with  constant 
width,  height  ranged  between  50  pm  and  500  pm  and  the  hydrau¬ 
lic  diameter  was  varying  between  1  pm  and  100  pm.  The  tap  water 


was  used  as  working  fluid  with  flow  rate  ranged  from  0.04  ml/min 
to  240  ml/min.  Reynolds  number  used  was  ranged  from  100  to 
5000.  The  results  reveal  that  the  conventional  laws  and  correla¬ 
tions  describing  the  flow  and  convective  heat  transfer  in  ducts  of 
large  dimension  were  directly  applicable  to  the  MCs  of  heights 
between  500  pm  and  50  pm.  In  addition,  the  transition  from  the 
laminar  to  turbulent  regime  occurred  in  a  range  of  Reynolds 
numbers  similar  to  that  observed  in  ducts  of  large  sizes. 

Bertsch  et  al.  [61]  investigated  experimentally  the  flow  boiling 
heat  transfer  with  the  refrigerants  R-134a  and  R-245fa  in  copper 
microchannel  cold  plate  evaporators.  Array  of  MCs  were  used  with 
hydraulic  diameter  1.09  mm  and  0.54  mm.  The  aspect  ratio  of  the 
rectangular  cross  section  of  the  channels  in  both  test  section  was 
2.5.  The  saturation  temperatures  used  were  ranged  from  8  "C  to 
30  °C,  mass  flux  from  20  kg/m2/s  to  350  kg/m2/s  and  the  heat  flux 
was  ranged  from  0  W/cm2  to  22  W/cm2.  The  results  indicate  that 
the  heat  transfer  coefficient  for  R-245fa  in  comparison  with  R-134a 
in  single-phase  flow  was  higher.  The  pool  boiling  heat  transfer  for 
R-245fa  was  lower  as  a  result  of  its  higher  molecular  mass  and 
surface  tension.  In  addition,  the  flow  boiling  heat  transfer  coeffi¬ 
cient  strongly  increased  with  the  increase  of  the  heat  flux  and  was 
dominated  by  nucleate  boiling.  It  also  shows  a  strong  dependence 
on  thermodynamic  vapor  quality  with  a  rapid  decrease  at  qualities 
above  0.5.  The  maximum  heat  transfer  coefficient  occurs  at  vapor 
qualities  between  0.1  and  0.5  depending  on  fluid,  geometry  and 
flow  conditions.  The  heat  transfer  coefficient  increased  weakly 
with  increases  mass  flux  and  was  almost  constant.  However,  a 
comparison  of  the  measurements  against  predictions  from  several 
correlations  showed  reasonable  agreement  only  with  a  very  few 
correlations  except  of  a  pool  boiling  equation  that  were  developed 
specifically  for  small  channels  did  not  predict  the  heat  transfer 
coefficient  better. 

Ngo  et  al.  [62]  investigated  experimentally  the  thermal-hydraulic 
characteristics  of  microchannel  heat  exchangers  (MCHEs)  with 
S-shaped  and  zigzag  fins  by  varying  Prandtl  number  widely  from 
0.75  to  2.2.  The  results  indicate  that  the  pressure-drop  factor  of  the 
MCHE  with  S-shaped  fins  was  4-5  times  less  than  that  of  MCHE 
with  zigzag  fins,  although  Nu  was  24-34%  less,  depending  on  the 
Reynolds  number  within  its  range.  However,  the  Nu  correlation  of 
the  MCHE  with  S-shaped  fins  reproduces  the  experimental  data  of 
overall  heat  transfer  coefficients  with  a  standard  deviation  of 
+  2.3%,  although  it  is  +  3.0%  for  the  MCHE  with  zigzag  fins.  The 
calculated  pressure  drops  obtained  from  pressure-drop  factor  corre¬ 
lations  agreed  with  the  experimental  data  within  a  standard  devia¬ 
tion  of  + 16.6%  and  + 13.5%  for  the  MCHEs  with  S-shaped  and 
zigzag  fins,  respectively. 

Wei  and  Joshi  [63]  investigated  experimentally  and  numerically 
the  sidewall  profile  effects  on  flow  and  heat  transfer  inside  micro¬ 
channels  with  diameter  ranged  from  53  pm  to  112  pm.  Numerical 
modeling  was  conducted  for  conjugate  heat  transfer  inside  MCs.  To 
simulate  the  practical  application  conditions  such  as  those  in 
electronic  cooling,  uniform  heat  flux  condition  was  applied  at  the 
bottom  region  while  the  interface  of  the  channel  was  assumed 
adiabatic.  Three  different  cases  were  considered  in  which  the 
channel  profiles  were  rectangular,  positively  tapered  and  re¬ 
entrant  respectively.  The  results  reveal  that  the  location  of  the 
maximum  velocity  deviates  from  the  mid-plane  along  the  depth 
direction  due  to  the  wall  taper.  In  addition,  for  constant  heat  flux 
condition,  serious  degradation  in  overall  heat  transfer  occurs  for 
channel  with  5  slope.  For  conjugate  heat  transfer,  the  degradation 
in  heat  transfer  was  much  less  and  re-entrant  channel  has  better 
heat  transfer  performance  than  positively  tapered  channel.  It  was 
concluded  that  for  MCs  with  large  aspect  ratio  the  transport 
characteristics  very  sensitive  to  the  channel  profiles,  therefore, 
careful  design  and  control  of  the  micro-fabrication  process  were 
crucial  to  the  performance  of  the  micro-scale  devices. 
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Fang  et  al.  [64]  investigated  experimentally  the  influences  of 
film  thickness  and  cross-sectional  geometry  on  hydrophilic  micro- 
channel  condensation,  with  aspect  ratios  ranged  from  1  to  5  and 
hydraulic  diameters  ranged  front  100  pm  through  300  pm.  The 
results  reveal  that  the  experimental  measurement  qualitatively 
matches  the  prediction  of  previous  theoretical  models  accounting 
for  the  surface  tension  effect,  which  highlights  the  importance  of 
surface  tension  force  and  channel  geometry  in  the  MC  condensa¬ 
tion.  In  addition,  the  pressure  drop  increased  linearly  with  both 
inlet  vapor  mass  flux  and  inlet  vapor  mass  flowrate.  Under  the 
same  vapor  flow  rate  and  aspect  ratio,  the  channel  with  the 
smaller  hydraulic  diameter  exhibits  the  higher  pressure  drop. 
While  the  pressure  drop  reduced  with  increases  the  hydraulic 
diameters  so  does  the  mean  heat  flux.  Therefore,  the  optimization 
of  hydrophilic  channel  geometry  for  film  condensation  enhance¬ 
ment  lies  in  the  trade-off  between  these  two  factors.  If  the 
pumping  power  of  the  compact  heat  exchanger  is  primary  con¬ 
cern,  a  larger  channel  was  more  desirable  considering  its  small 
pressure  drop.  In  contrast,  if  heat  exchanger  size  or  channel  length 
is  the  primary  limit,  a  smaller  channel  was  more  favorable  due  to 
its  high  mean  heat  flux. 

Bogojevic  et  al.  [65]  investigated  experimentally  the  effects  of  non- 
uniform  heating  on  flow  boiling  instabilities  in  microchannel-based 
heat  sink  using  water  as  the  cooling  liquid.  The  results  indicate  that 
the  onset  of  boiling  caused  metastable  non-uniform  flow  distribution, 
with  boiling  in  some  channels  and  single-phase  flow  in  others,  leading 
to  transverse  variations  in  temperature  that  were  less  severe  when  the 
hotspot  was  at  the  downstream  position.  The  non-uniformity  was 
attributed  to  statistical  variations  in  the  small  numbers  of  stable 
bubble  nucleation  sites  per  channel.  In  addition,  the  oscillations  were 
superimposed  on  the  non-uniform  flow  distributions.  At  low  hotspot 
heat  flux,  there  were  high  amplitude  oscillations  at  a  dominant  low 
frequency  (HALF).  With  increasing  heat  flux,  this  changed  progres¬ 


sively  to  lower  amplitude  oscillations  at  higher  frequency  (LAHF).  In 
the  HALF  regime,  synchronized  fluctuations  in  the  wall  temperature 
occurred  in  channels  with  and  without  boiling,  indicating  that  the 
oscillations  were  not  driven  by  the  growth  of  confined  bubbles  in 
individual  channels. 

Wang  et  al.  [66]  investigated  experimentally  the  effects  of  inlet/ 
outlet  configurations  on  flow  boiling  instability  in  parallel  microchan¬ 
nels  have  a  30  mm  length  and  a  hydraulic  diameter  of  186  pm.  Three 
types  of  inlet/outlet  configurations  were  obtained.  In  Type-A  connec¬ 
tion,  inlet  and  outlet  conduits  were  perpendicular  to  the  parallel  MCs 
in  the  test  section.  In  Type-B  connection,  fluid  could  flow  into  and  exit 
from  the  parallel  MCs  without  any  restriction,  and  in  Type-C  connec¬ 
tion,  fluid  entered  each  of  the  MCs  through  an  inlet  restriction.  The 
results  indicate  that  the  amplitudes  of  temperature  and  pressure 
oscillations  in  Type-B  connection  was  much  smaller  than  those  in 
Type-A  connection  under  the  same  heat  flux  and  mass  flux  conditions. 
Nearly  steady  flow  boiling  exists  in  the  parallel  MCs  with  Type-C 
connection  under  the  experimental  conditions  as  shown  in  Fig.  6. 
Therefore,  this  configuration  was  recommended  for  high-heat  flux  MC 
applications.  The  experimental  data  of  pressure  drop  was  higher  and 
heat  transfer  coefficients  was  lower  for  boiling  flow  at  high  vapor 
quality  in  MCs  than  those  predicted  from  correlation  equations  for 
boiling  flow  in  MCs,  due  to  local  dryout. 

Hernando  et  al.  [67]  investigated  experimentally  fluid  flow 
and  heat  transfer  in  a  single-phase  liquid  flow  micro-heat 
exchanger,  characterized  by  MCs  of  100  x  100  pm  and 
200  x  200  pm  square  cross  section.  The  deionized  water  was 
used  as  a  working  fluid.  The  results  show  good  agreement  with 
the  general  theory  and  no  special  effect  related  to  the  small 
dimension  of  the  channels  was  observed.  In  addition,  no  heat 
transfer  increment  appeared  as  a  result  of  the  small  channel 
dimensions.  For  very  low  Reynolds  numbers,  the  results  show 
discrepancies  with  the  available  models  and  correlations,  which 
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(i)  Arrangement  of  parallel  microchannels  with 
restrictive  inlet 


(ii)  Sketch  of  the  microchannel  with  a 
restrictive  inlet 


Fig.  6.  Parallel  microchannels  with  three  different  inlet/outlet  connections  [66].  (a)  Type-A  connections:  flow  entering  and  exiting  from  paral  microchannels  with 
restrications  because  inlet/outlet  conduit  perpendicular  to  microchannels,  (b)  Type-B  connection:  flow  entering  to  and  exiting  from  microchannels  freely  withuot  restriction 
and  (c)  Type-C  connection:  flow  entering  with  restriction  and  exiting  without  restriction  in  microchannels. 
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Fig.  7.  geometrical  configuration  [68], 


overestimate  the  convection  coefficients  which  related  with 
the  adequacy  of  the  boundary  conditions.  On  the  other  hand, 
the  plate  conduction  thermal  resistance  was  a  major  restriction 
for  micro-heat  exchanger  performance.  The  plate  thickness 
and  plate  material  were  critical  in  the  design  of  micro-heat 
exchangers. 


3.2.  Numerical  studies  on  MC 

Yang  and  Lai  [68]  investigated  numerically  the  effects  of  forced 
convection  flow  of  Al203-water  nanofluids  in  the  radial  flow 
cooling  system  using  a  single-phase  approach.  The  geometrical 
configuration  of  the  confined  radial  flow  was  used  shown  in  Fig.  7. 
The  two  plates  separated  by  a  gap  were  parallel  and  stationary. 
The  lower  disk  (with  nozzle)  acts  as  the  insulated  confinement 
plate,  while  the  upper  one  (wafer)  was  subjected  to  a  constant 
heat  flux.  The  numerical  simulations  were  performed  based  on  the 
following  conditions:  volume  concentration  of  47  nm  A1203  nano¬ 
particle  suspended  within  water  is  from  0%  to  10%,  flow  Reynolds 
number  range  from  300  to  900  with  the  heat  flux  2438  W/m2  and 
3900  W/m2.  The  inlet  fluid  temperature  T0=296K.  The  results 
indicate  that  the  heat  transfer  coefficient  increased  with  the 
increases  of  Reynolds  number  and  particle  volume  fraction. 

Wu  et  al.  [69]  proposed  a  new  correlation  for  the  saturated- 
flow  boiling  critical  heat  flux  (CHF)  in  microchannel  for  both  multi 
and  single  channel  configurations.  The  saturated  CHF  correlation 
was  proposed  by  using  the  boiling  number,  length  to  diameter 
ratio  and  the  exit  quality.  The  results  reveal  that  the  boiling 
number  at  CHF  decreased  greatly  with  a  small  increase  in  length 
to  diameter  ratio  Lh  /dhe  when  this  ratio  <  50.  The  transition  was 
appeared  in  the  range  of  50  <  Lh  !dhe  <  150.  This  correlation  can 
predict  the  overall  MC  database  accurately.  This  method  almost 
predicts  97%  of  the  non-aqueous  data  and  94%  of  the  water  data 
within  the  +  30%  error  band. 

Shevade  and  Rahman  [70]  investigated  numerically  the  con¬ 
vective  heat  transfer  in  rectangular  and  square  cross  sections  MCs 
for  volumetric  heat  generation.  Gadolinium  was  used  as  the 
substrate  material  and  water  as  the  working  fluid.  The  results 
reveal  that  the  Nusselt  number  decreased  along  the  x  or  y-axis 
from  the  mid-section  of  the  channel  to  the  corner  due  to  smaller 
fluid  velocity  which  caused  by  larger  frictional  resistance  at  the 
corner.  In  rectangular  channels,  the  heat  transfer  coefficient  values 
are  higher  on  the  longer  side  than  on  the  shorter  side.  At  a  smaller 
flow  rate,  the  local  Nusselt  number  decreased  and  the  outlet 
temperature  increased  as  the  low  velocity  fluid  remains  in  contact 
with  the  solid  for  a  longer  period  of  time.  The  pressure  drop  in  the 
channel  increased  as  the  Reynolds  number  increased.  Nusselt 
number  increased  as  the  channel  dimensions  increased.  For  the 
same  channel,  Nusselt  number  increased  with  Reynolds  number. 


Gamrat  et  al.  [71]  investigated  numerically  the  roughness 
effects  on  heat  transfer  in  thermally  fully-developed  laminar  flows 
through  microchannels.  The  models  considered  parallelepipedic 
elements  periodically  distributed  on  the  plane  walls  of  a  micro- 
channel  of  infinite  span.  Three-dimensional  numerical  simulations 
were  first  conducted  in  a  computational  domain  extending  over 
one  wavelength  in  streamwise  and  spanwise  directions.  The 
numerical  results  reveal  that  the  Poiseuille  number  P„  and  the 
Nusselt  number  (Nu)  increased  with  the  relative  roughness  and 
nearly  independent  of  Reynolds  number  in  the  laminar  regime 
(Re  <  2000).  In  addition,  the  roughness  increase  the  friction  factor 
more  than  the  heat  transfer  coefficient  for  a  relative  roughness 
height  expected  in  the  fabrication  of  MCs. 

Mlcak  et  al.  [72]  investigated  numerically  the  heat  transfer  and 
laminar  fluid  flow  in  an  array  of  parallel  microchannels  etched  on  a 
silicon  substrate  with  water  as  the  circulating  fluid.  The  hydraulic 
diameter  used  was  85.6  pm  and  aspect  rations  ranged  from  0.10  to 
1.0.  Reynolds  number  used  was  ranged  from  50  to  400.  The  results 
reveal  that  the  linearly  increasing  friction  coefficient  that  was 
observed  between  Reynolds  numbers  50  and  400  as  a  result  of 
increased  hydraulic  entrance  length.  Friction  coefficient  values  were 
found  to  monotonically  increase  for  aspect  ratios  approaching  0.10. 
The  values  of  inlet  and  outlet  thermal  resistance  decreased  for 
aspect  ratios  approaching  0.10,  for  the  same  hydraulic  diameter  the 
inlet  and  outlet  thermal  resistance  and  friction  coefficient  values 
changed  very  little  for  aspect  ratios  larger  than  0.50. 

Chen  et  al.  [73]  investigated  numerically  the  heat  transfer  and 
fluid  flow  in  noncircular  microchannel  heat  sink.  The  cooling  fluid  in 
the  channel  used  was  deionized  water.  The  geometries  for  the 
triangular,  rectangular  and  trapezoidal  channel  were  used,  129.6  pm, 
40  pm  and  158  pm,  respectively.  The  results  reveal  that,  the  circum¬ 
ferential  average  Nusselt  number  has  a  much  higher  value  at  the  inlet 
region  but  it  quickly  approaches  the  constant  fully  developed  value. 
The  liquid  and  solid  temperature  increased  along  the  flow  direction. 
In  addition,  the  Poiseuille  number  in  MCs  with  a  certain  cross  section 
remained  almost  constant  and  independent  with  the  Reynolds 
number  for  fully  developed  laminar  flow  which  was  similar  to  fluid 
flow  in  MCs.  Triangular  MC  heat  sink  has  the  highest  thermal 
efficiency  among  the  three  heat  sinks,  the  trapezoidal  MC  was  inferior 
to  triangular  but  superior  to  rectangular. 

Zhu  and  Liao  [74]  investigated  theoretically  the  laminar  forced 
heat  transfer  for  a  gas  flow  through  a  microchannel  of  arbitrary 
cross  section  with  axial  constant  heat  flux  and  circumferential 
varied  wall  temperature  boundary  condition  in  the  slip  flow  and 
temperature-jump  regime.  The  influences  of  Knudsen  number,  the 
channel  aspect  ratio,  and  the  thermal  boundary  condition  on  heat 
transfer  behavior  were  estimated.  The  results  indicate  that  for  a 
rectangular  and  triangular  microchannel,  the  average  Nusselt 
numbers  were  smaller  than  those  of  a  macrochannel  with  the 
same  thermal  boundary  conditions  and  it  decreases  with  increas¬ 
ing  Knudsen  number.  Knudsen  number  has  more  significant 
influences  on  the  average  Nusselt  number  of  the  rectangular 
microchannel  with  a  smaller  channel  aspect  ratio.  The  varying 
trends  of  the  average  Nusselt  number  with  the  channel  aspect 
ratio  were  similar  at  various  Kn  numbers.  The  channel  aspect  ratio 
has  few  effects  on  the  average  Nusselt  number  ratio  at  a  fixed 
Knudsen  number  within  a  moderate  range  of  channel  aspect  ratio. 

Ji  et  al.  [75]  investigated  numerically  the  wall  roughness  on 
gaseous  flow  and  heat  transfer  in  a  microchannel  with  length 
range  from  1  pm  to  100  pm.  The  effect  was  examined  for  gas  flows 
under  inlet  Mach  number  ranged  from  0.0055  to  0.202.  The 
Reynolds  number  used  was  ranged  from  0.001  to  100.  The  results 
indicated  that  the  rough  elements  restrict  the  upstream  gas  flow 
and  cause  more  pressure  drop  and  increase  the  Poiseuille  number 
due  to  the  obstruction  effect  and  formation  of  recirculation 
regions.  Rarefied  flow  was  sensitive  to  the  shape  of  roughness 
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elements.  The  average  Poiseuille  number  increased  not  only  as  the 
roughness  height  increases  but  also  as  the  spacing  distance 
between  the  roughness  elements  was  decreased.  However,  this 
effect  was  much  stronger  on  the  gas  flow  with  a  lower  Knudsen 
number.  In  addition,  for  a  highly  compressible  flow,  the  average 
Nusselt  number  does  not  change  much  or  increases  slightly.  The 
existence  of  wall  roughness  enhances  the  heat  transfer. 

Sui  et  al.  [76]  investigated  numerically  laminar  liquid  water 
flow  and  heat  transfer  in  three-dimensional  wavy  MCs  with 
rectangular  cross  section  under  constant  wall  heat  flux.  In  the 
conjugate  simulation,  a  uniform  heat  flux  was  applied  at  the 
bottom  wall  of  the  silicon  substrate  and  on  the  upper  surface 
adiabatic  boundary  condition  was  employed.  The  results  reveal 
that  when  liquid  coolant  flows  through  the  wavy  MCs,  second¬ 
ary  flow  (Dean  vortices)  was  generated.  The  quantity  and  the 
location  of  the  vortices  may  change  along  the  flow  direction, 
leading  to  chaotic  advection,  which  can  greatly  enhance  the 
convective  fluid  mixing,  and  thus  the  heat  transfer  performance 
of  the  wavy  MCs  was  much  better  than  that  of  straight  MCs  with 
the  same  cross  section.  In  addition,  the  pressure  drop  penalty 
of  the  wavy  MCs  was  much  smaller  than  the  heat  transfer 
enhancement. 

Muhammad  et  al.  [77]  investigated  numerically  the  convective 
heat  transfer  in  a  composite  trapezoidal  microchannel  with  length 
23  mm  during  magnetic  heating.  The  channel  was  formed  by 
etching  a  silicon  100  pm  wafer  and  bonding  it  with  a  slab  of 
gadolinium  with  height  varied  between  1.5  mm  and  5  mm.  Rey¬ 
nolds  number  used  was  ranged  from  1000  to  3000.  The  roughness 
height  used  was  0.35  nm  for  silicon  MC  surface  and  1  nm 
for  gadolinium.  Magnetic  field  used  was  ranged  from  2.5  T  to 
10  T.  The  results  reveal  that  for  a  smaller  flow  rate,  the  outlet 
temperature  increased  as  a  lower  velocity  fluid  remains  in 
contact  with  the  solid  for  a  longer  period  of  time.  The  peripheral 
average  heat  transfer  coefficient  and  Nusselt  number  decreased 
along  the  length  of  the  channel  due  to  the  development  of 
thermal  boundary  layer.  Increasing  the  channel  spacing  slightly 
increased  the  local  Nusselt  number  all  along  the  length  of  the 
channel.  The  Nusselt  number  increased  with  decrease  of  channel 
hydraulic  diameter  at  any  given  Reynolds  number  due  to  higher 
fluid  velocity.  Compared  to  helium  and  FC-72,  water  provides 
better  thermal  performance  with  higher  values  of  Nusselt 
number. 

Liu  et  al.  [78]  investigated  numerically  fluid  flow  and  heat 
transfer  in  MC  cooling  passages.  Two-dimensional  simulation  was 
performed  for  low  Reynolds  number  flow  of  liquid  water  with 
diameter  D=100  pm  of  single  channel  subjected  to  localized  heat 
flux  boundary  conditions.  The  results  reveal  that  the  heat  transfer 
enhancement  due  to  viscosity-variation  was  pronounced,  though 
the  axial  conduction  introduced  by  thermal-conductivity-variation 
was  insignificant  unless  for  the  cases  with  very  low  Reynolds 
numbers. 

Renksizbulut  et  al.  [79]  investigated  numerically  rarefied  gas 
flow  and  heat  transfer  in  the  entrance  region  of  rectangular  MC  in 
the  slip-flow  regime.  The  Reynolds  number  used  was  ranged  from 
0.1  to  10,  channel  aspect  ratio  ranged  from  0  to  1  and  the  Knudsen 
number  was  Kn  <  0.1.  The  results  indicate  that  the  slip-flow  regime 
was  associated  with  very  large  reductions  in  wall  friction  and  heat 
transfer  in  the  entrance  region.  Fully  developed  values  of  friction 
coefficient  and  heat  transfer  rates  were  more  influenced  with 
rarefaction  effects  in  a  flat  channel  than  a  square  duct  due  to  the 
major  influence  of  channel  corners  where  rarefaction  effects 
largely  vanish. 

Niazmand  et  al.  [80]  investigated  numerically  the  velocity  and 
temperature  fields  in  the  slip-flow  regime  in  trapezoidal  MCs  with 
constant  wall  temperatures.  The  range  of  channel  aspect  ratio  used 
was  from  0.25  to  2,  side  angles  from  30°  to  90  and  the  Reynolds 


number  ranged  from  0.1  to  10.  The  prandtl  Number  used  was 
Pr=  1.  The  results  reveal  that  the  major  reductions  in  the  friction 
and  heat  transfer  coefficients  were  observed  in  the  entrance 
region  due  to  large  amounts  of  velocity-slip  and  temperature- 
jump.  The  friction  coefficient  decreased  strongly  both  with 
increasing  the  Knudsen  number  and  aspect  ratio  but  has  a  weaker 
dependence  on  the  side  angle. 

Khadem  et  al.  [81]  investigated  numerically  fluid  flow  and  heat 
transfer  characteristics  in  MCs  at  slip  flow  regime  with  considera¬ 
tion  of  slip  and  temperature-jump.  The  results  reveal  that  the 
rough  elements  restrict  the  upstream  gas  flow  and  cause  more 
pressure  drop  and  increases  in  the  Poiseuille  number  due  to  the 
obstruction  effect.  Roughness  has  more  significant  effect  on  higher 
Knudsen  number  flows  with  higher  relative  roughness.  The  chan¬ 
ging  in  roughness  shape  leads  to  extremely  different  flow  patterns 
except  for  small  relative  roughness. 

Zade  et  al.  [82]  investigated  numerically  the  effects  of 
variable  physical  properties  on  the  flow  and  heat  transfer 
characteristics  of  simultaneously  developing  slip-flow  in  rec¬ 
tangular  MCs  with  constant  wall  temperature.  The  Reynolds 
number  used  was  set  equal  to  0.1.  The  inlet  gas  temperature  in 
all  simulations  was  specified  to  be  T,n  =  300K,  and  the  wall 
temperature  was  Tw= 350  K.  The  results  indicate  that  in  the 
entrance  region,  Nusselt  number  was  only  affected  by  the 
change  in  the  Knudsen  number  and  was  higher  in  variable- 
density  slip-flow  cases.  The  Nusselt  number  was  generally 
higher  in  the  variable-property  simulations.  In  the  fully- 
developed  region,  since  the  effect  of  density  change  diminishes, 
friction  coefficient  and  Nusselt  number  predictions  of  the 
variable-property  simulations  were  higher  in  comparison  to 
those  of  the  constant-property  cases. 

Shokouhmand  and  Bigham  [83]  investigated  numerically  the 
effects  of  Knudsen  number  and  geometry  on  thermal  and  hydro- 
dynamic  characteristics  of  flow  in  the  wavy  microchannel  at 
constant  Reynolds  number.  The  results  reveal  that  the  Nusselt 
number  and  Cp  Re  decreased  with  Knudsen  number.  In  the 
divergent  parts,  Nusselt  number  experiences  a  rapid  decrease. 
The  very  high  heat  transfer  rate  and  Cp  Re  at  the  entrance  declines 
rapidly  as  the  thermally-hydraulically  developing  flow  approaches 
fully  developed  flow.  In  addition,  the  effect  of  viscous  dissipation 
has  a  considerable  effect  in  MCs  and  this  effect  can  be  more 
significant  by  increasing  Knudsen  number. 

Mohammed  et  al.  [84,85]  investigated  numerically  the  influ¬ 
ence  of  channel  shapes  (zigzag,  curvy,  wavy)  on  the  thermal  and 
hydraulic  performance  of  microchannel  heat  sink  (MCHS).  The 
Reynolds  number  used  was  ranged  from  100  to  1000  for  all 
channel  flow  shapes.  Water  was  used  as  a  working  fluid.  The  inlet 
temperature  of  water  used  was  293  K  and  the  heat  flux  at  the  top 
plate  of  MCHS  was  100  W/m2.  The  results  reveal  that,  the  tem¬ 
perature  and  the  heat  transfer  coefficient  of  the  zigzag  MCHS  was 
the  least  and  greatest,  respectively,  comparing  to  other  channel 
shapes.  The  second  best  thermal  performance  was  provided  by 
curvy  channels  after  the  wavy  channels.  The  zigzag  MCHS  was  the 
best  channel  for  the  best  thermal  performance,  while  the  step 
MCHS  was  the  best  channel  for  the  hydraulic  performance  with 
moderate  degradation  of  heat  transfer  compared  to  conventional 
straight  MCHS 

Kosar  [86]  investigated  numerically  the  effect  of  substrate 
thickness  and  material  on  heat  transfer  in  microchannels  heat 
sink.  The  size  of  the  MCs  used  (200  pm  x  200  pm,  5  cm  long). 
The  substrate  thickness  was  ranged  from  100  pm  to  1000  pm 
and  the  MEMS  materials  used  (polyimide,  Silica  Glass,  Quartz, 
steel,  silicon  and  copper).  Temperature,  pressure  and  velocity 
distributions  across  a  MC  heat  sink  were  found  under  a  broad 
range  of  Reynolds  number  and  heat  transfer  using  water  as  a 
working  fluid.  The  results  reveal  that  a  very  good  agreement 
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between  the  model  results  and  the  theoretical  prediction 
corresponding  to  constant  heat  flux  boundary  condition  for 
Reynolds  number  Re  <  976  was  observed,  beyond  Re  =  1000 
Nusselt  number  results  were  greater  than  the  theoretical  value 
of  Nusselt  number  for  developed  flows  due  to  developing 
effects.  The  Nusselt  number  was  lower  for  low  thermal  con¬ 
ductivity  materials  than  the  Nusselt  number  obtained  from  high 
conductivity  materials.  The  weaker  dependence  of  substrate 
thickness  on  Nusselt  number  compared  to  the  effect  of  thermal 
conductivity  was  because  of  the  existing  constraint  in  the  ratio 
of  substrate  thickness  to  the  MC  size. 

Mohammed  et  al.  [87,88]  investigated  numerically  the  influ¬ 
ences  of  nanofluids  on  parallel  flow  square  and  rectangular 
microchannel  heat  exchanger  (MCHE)  performance.  Four  differ¬ 
ent  nanofluids  (aluminum  oxide  (A1203),  silicon  dioxide  (Si02), 
silver  (Ag)  and  titanium  dioxide  (Ti02)  with  three  nanoparticle 
volume  fraction  of  2%,  5%  and  10%  were  used  with  water  as 
base  fluid.  Reynolds  number  used  was  ranged  from  100  to  800. 
The  results  reveal  that  as  the  nanoparticle  volume  fraction 
increased,  the  bulk  average  temperature  of  the  cold  fluid 
increased  and  the  heat  transfer  rate  decreased.  On  the  other 
hand,  as  the  Reynolds  number  increased,  the  average  bulk 
temperature  of  the  cold  fluid  decreased  and  the  heat  transfer 
rate  increased.  In  addition,  aluminum  oxide  nanofluid  has  the 
highest  heat  transfer  coefficient  among  all  the  nanofluids 
tested.  The  presence  of  nanoparticles  increase  the  pressure 
drop  slightly  with  silver  nanofluid  has  the  lowest  pressure  drop, 
while  the  highest  difference  of  pressure  drop  from  water  was 
lower  than  5%. 

Revillin  and  Thome  [89]  reported  a  theoretical  model  for  the 
prediction  of  the  critical  heat  flux  (CHF)  of  refrigerants  flowing  in 
heated  microchannel.  The  prediction  was  based  on  the  local 
dryout  of  the  liquid  film  in  annular  flow  occurring  when  the  film 
thickness  becomes  equal  to  the  interfacial  wave  height  during 
evaporation.  The  model  was  based  on  the  conservation  of  mass, 
momentum  and  energy,  the  Laplace-Young  equation  and  a  semi- 
empirical  expression  for  the  height  of  interfacial  waves.  The  results 
reveal  that  more  than  96%  of  the  data  was  predicted  within  a 
+  20%  error  band  and  a  mean  absolute  error  of  8%.  It  is  possible  to 
predict  CHF  data  from  a  third  laboratory  for  water  and  R-113 
flowing  in  rectangular  (using  the  width  of  the  channel  as  the 
characteristic  dimension)  and  circular  MCHS  with  multiple 
channels. 

Revellin  et  al.  [90]  reported  status  of  prediction  methods  for 
critical  heat  flux  in  mini  and  microchannel.  2996  data  point 
collected  from  19  different  laboratories  includes  nine  different 
fluids  (R-134a,  R-245fa,  R-236fa,  R-123,  R-32,  R-113,  nitrogen,  C02 
and  water)  for  a  wide  range  of  experimental  conditions  were 
compared  to  6  different  correlations  and  1  theoretically  based 
model  to  provide  the  best  CHF  prediction  method.  It  predicts 
almost  86%  of  the  CHF  data  for  non-aqueous  fluids  within  a  30% 
error  band.  This  method  predicts  almost  83%  of  the  CHF  data  for 
water  within  a  30%  error  band. 

Rosa  et  al.  [91  ]  reviewed  an  experimental  and  numerical  result 
from  different  references  about  the  importance  of  scaling  effects 
(entrance  effects,  conjugate  heat  transfer,  viscous  heating,  electric 
double  layer  (EDL)  effects,  temperature  dependent  properties, 
surface  roughness,  rarefaction  and  compressibility  Effects)  in  the 
single-phase  heat  transfer  in  microchannels.  Their  review  con¬ 
cludes  with  a  concise  set  of  recommendations  for  purposes  of 
performance  and  design.  For  single  channels,  available  correlations 
for  macro-channels  can  also  give  reliable  predictions  at  the 
micro-scale,  but  only  if  all  the  scaling  effects  can  be  considered 
negligible.  Otherwise,  when  scaling  effects  cannot  be  neglected  or 
for  the  case  of  heat  exchangers  with  parallel  channels,  suitable 
numerical  simulations  may  be  the  sole  alternative  to  carefully 


designed  experiments  to  evaluate  the  heat  transfer  rates.  A 
summary  for  the  numerical  and  experimental  studies  in  MC  is 
shown  in  Table  2. 

4.  Nanofluids 

4.1.  Preparation  of  nanofulids 

The  powder  form  nanoparticles  which  disperse  in  host  liquids 
are  called  nanofluids.  Nanofluids  can  be  produced  by  two  techni¬ 
ques;  the  two-step  (double-step)  method,  and  one-step  (single-step) 
method.  These  methods  have  been  utilized  using  different  types  of 
chemical  and  physical  techniques  to  make  sure  that  the  solid-liquid 
mixture  is  stable  to  avoid  agglomeration,  additional  flow  resistance, 
possible  erosion  and  clogging,  poor  thermal  conductivity,  and  poor 
heat  transfer.  The  two-step  method  is  done  by  producing  the 
nanoparticles  powder  initially  as  introduced  in  the  previous  section, 
and  then  disperses  them  into  a  host  liquid.  However,  in  one-step 
method  the  nanoparticles  are  simultaneously  made  and  directly 
dispersed  into  the  base  fluid  [92].  It  is  noticed  in  the  literature  that 
nanofluids  with  oxide  nanoparticles  and  carbon  nanotubes  are 
produced  well  by  the  two-step  method,  while  it  is  not  suitable  for 
nanofluids  with  metallic  nanoparticles. 

Chopkar  et  al.  [93]  presented  a  simple  method  of  preparing 
nanofluids  comprising  a  small  amount  (  <  1  vol%)  of  nanocrystal¬ 
line  metallic  particles  (Al70Cu30)  as  a  stable  colloidal  dispersion  in 
ethylene  glycol.  The  size/microstructure  of  the  nanoparticles 
was  characterized  by  X-ray  diffraction  and  transmission  electron 
microscopy,  and  the  thermal  conductivity  of  the  nanofluid  was 
measured  by  using  a  modified  thermal  comparator.  The  results 
showed  that  the  thermal  conductivity  measurements,  using  a 
purpose-built  thermal  comparator  device,  revealed  to  that  the 
conductivity  of  this  nanofluid  is  significantly  greater  (1.2-2  times) 
than  that  of  the  base  fluid.  Moreover,  the  increase  in  conductivity 
is  a  function  of  volume  fraction,  size  and  thermal  property  of  the 
solid  suspension.  The  results  also  demonstrated  that  the  heat 
removal  rate  or  efficiency  of  the  nanofluid  is  better  than  that  of  the 
base  fluid,  especially  at  lower  temperatures.  The  results  revealed 
to  that  the  nanofluid  could  be  useful  in  microthermal  heat  transfer 
applications  including  in  MEMS  and  automobiles. 

4.2.  Experimental  studies  on  nanofluids 
4.2.1.  Thermal  conductivity 

It  is  reported  that  many  experimental  studies  have  been  done  on 
thermal  conductivity  of  nanofluids  because  it  is  a  primary  assess¬ 
ment  of  the  heat  transfer  performance  of  nanofluids.  The  physical 
mechanism  for  thermal  conductivity  enhancement  of  nanofluids  is 
not  well  understood  up  to  date.  Li  et  al.  [94]  presented  investigation 
of  the  effective  thermal  conductivity  of  Al203/water  nanofluids  with 
normal  diameters  of  47  nm  at  different  volume  fractions  (0.5%,  2%, 
4%,  and  6%).  In  this  investigation  two  different  methodologies  were 
utilized:  a  transient  hot-wire  method  and  a  steady-state  cut-bar 
method.  The  comparison  of  the  measured  data  obtained  by  using 
two  different  experimental  systems  at  room  temperature  was 
conducted  and  the  experimental  data  at  higher  temperatures  were 
obtained  with  steady-state  cut-bar  method  and  compared  with 
previously  reported  data  obtained  using  a  transient  hot-wire 
method.  The  results  showed  that  at  room  temperature,  both  the 
steady-state  cut-bar  and  transient  hot-wire  methods  result  in 
closely  identical  values  for  the  effective  thermal  conductivity  of 
the  nanofluids  tested.  Furthermore,  at  higher  temperatures,  the 
onset  of  natural  convection  results  in  larger  measured  effective 
thermal  conductivities  for  the  hot-wire  method  than  those  obtained 
using  the  steady-state  cut-bar  method.  The  results  also  confirmed 


Table  2 

Summary  of  experimental  and  numerical  studies  for  MC. 


Researchers  Diameter/flow  nature  Reynolds  Fluid/gas 

number 


Findings 


Jung  et  al.  [29] 

50-100  pm/laminar 

5-300 

Water  Water+EG/Al203 

Cheng  Wang 

25  mm/laminar 

- 

NaCl  and  diluted  glycerol 

et  al.  [30] 

Agarwal  et  al. 

0.424-0.939  mm/Laminar 

- 

R134a 

[31] 

Roland  et  al.  [32] 

700-*200  pm/laminar  and 
turbulent 

200- 

8000 

Demineralized  water 

Wu  and  Cheng 

L/D= 285.41 -451.40/laminar 

0-100 

Water 

[33] 

Park  and  Punch 

106-307  pm/laminar 

69-800 

Deionized  water 

[34] 

Quan  et  al.  [35] 

127-173  pm/laminar 

240 

Water 

Lee  et  al.  [36] 

194-534  pm 

300- 

3500 

Deionized  water 

Harirchian  and 
Garimella  [37] 

Width  100-5850  pm  with  depth 
400  pm/laminar 

- 

Fluorinert  FC-77 

Ergu  et  al.  [38] 

Width=3.70  mm, 
height =0.107  mm, 
length =35  mm/laminar 

100-845 

Distilled  water 

Owhaib  et  al. 

0.862-1.7  mm/laminar 

- 

R134a 

[39] 

Chiu  et  al.  [40] 

W/H=  1.67-14.29/laminar 

50-1000 

Water 

Shen  et  al.  [41  ] 

Width =300  pm,  Height=800  pm/ 
laminar 

162-1257 

Deinonized  water 

Barber  et  al.  [42] 

727  pm/laminar  to  turbulent 

- 

Not  given 

Seng  Lee  and 
Garimella  [43] 

Width  =  102-997  pm, 
depth=400  pm 

- 

Deionized  water 

Celata  et  al.  [44] 

480  pm/laminar 

- 

FC-72 

Megahed  [45] 

248  pm/laminar 

- 

FC-72 

Diaz  and 

Schmidt  [46] 

Height=0.3  mm,  width  =  12.7  mm/ 
laminar 

- 

Water/ethanol 

Chen  et  al.  [47] 

100-250  pm/laminar 

200-1700 

Water 

Betz  and  Attinger 

2  mm/laminar 

160-1580 

Water 

[48] 

Wojtan  et  al.  [49] 

0.5  and  0.8  mm/laminar 

- 

R-134a  and  R-245fa 

Asthana  et  al. 

100  pm/laminar 

- 

Water/oil 

[50] 

Wang  and  Cheng 

155  pm/laminar 

- 

Deionized  water 

[51] 

Wu  et  al.  [52] 

90.6  pm/laminar 

300-900 

Water 

Steinke  and 
Kandlikar  [53] 

Width =200  pm,  depth =250  pm 
and  a  length=10  mm 

14-789 

Water 

Wu  et  al.  [54] 

Width =300  pm,  depth  =  800  pm/ 
laminar 

- 

Water/methanol 

Zhang  et  al.  [55] 

Width =61  pm,  height =272  pm, 
length =15  mm/laminar 

- 

Deionized  water 

Kohl  et  al.  [56] 

25-100  pm/laminar  and  turbulent 

4.9- 

18,814 

Water  and  air 

Rosengarten 

100  mm/laminar 

0-70 

Distilled  water 

et  al.  [57] 

Tang  et  al.  [58] 

52-300  pm/laminar 

3-6300 

Nitrogen  and  helium 

Huh  et  al.  [59] 

103.5  pm/laminar 

19.2-38.3 

Deionized  water 

The  convective  heat  transfer  coefficient  of  the  A1203  was  increased  up  to  32%  compared  to  the  distilled  water 
Demonstrated  a  new  method  to  solve  the  unmatched  viscosity  problem  of  the  two-fluid  flow  in  MC 

Developed  a  correlation  for  heat  transfer  during  condensation  of  refrigerant  R134a  in  horizontal  non-circular  MC 

A  reduction  in  the  Nusselt  number  for  laminar  flow  at  the  microscale  is  occurred 

Nusselt  number  increased  almost  linearly  with  the  rise  of  Reynolds  number 

The  conventional  theory  for  fully-developed  flow  was  applicable  within  the  range  of  the  values  of  parameters  they  use 

Heat  transfer  coefficient  in  MC  increased  with  mass  flux  and  quality  and  decreased  with  the  hydraulic  diameter 
Heat  transfer  coefficient  increased  with  the  decrease  of  the  channel  size  depending  on  the  flow  rate 

A  strong  dependence  of  pressure  drop  on  both  channel  size  and  mass  flux  was  observed 

The  classical  laminar  flow  equations  can  be  used  for  the  calculation  of  the  friction  factor  in  MCs 


Heat  transfer  coefficient  was  higher  for  smaller  diameter  tubes.  An  increase  in  the  system  pressure  improves  the  heat  transfer 
performance 

Heat  transfer  enhancement  was  more  significant  on  MCs  with  high  aspect  ratio 

Nusselt  numbers  were  significantly  lower  than  the  conventional  theory  predictions  for  two  main  factors,  the  cross-sectional 
aspect  ratio  and  the  surface  roughness 

Flow,  pressure  and  temperature  conditions  leading  to  nucleate  boiling,  flow  instabilities  and  transition  regimes  during  flow 
boiling  in  a  MC 

The  saturated  heat  transfer  coefficient  becomes  largely  insensitive  to  heat  flux 

Heat  transfer  coefficient  was  rather  independent  of  heat  flux 

The  instability  increased  with  decreasing  the  mass  flux 

Heat  transfer  coefficient  increased  with  increasing  quality  at  low  heat  fluxes 

The  condensation  heat  transfer  was  enhanced  by  a  reduction  in  the  channel  scale 
The  segmented  flow  increases  the  Nusselt  number  of  laminar  flow  by  more  than  100% 

A  strong  dependence  of  CHF  on  mass  velocity,  heated  length  and  MC  diameter 

Increased  heat  removal  increased  the  pressure  drop,  which  is  strongly  affected  by  the  fluid  temperature 

The  stable  flow  boiling  occurred  in  a  single  MC  depending  on  the  vapor  quality 

Cross-sectional  shape  of  the  MC  plays  a  significant  role  on  the  instability  of  condensation  flow 
The  conventional  Stokes  and  Poiseuille  flow  theories  apply  for  single-phase  liquid  flow  in  MC  flows 

The  wall  temperatures  generally  match  the  oscillation  trends  of  the  inlet  pressure  and  the  outlet  temperature,  but  a  time  delay 
occurs  due  to  the  thermal  inertia  of  the  copper  block 

The  predictions  from  the  flow  oscillation  model  were  agreed  well  with  the  experimental  pressure-drop  observations  across  a 
flow  meter  and  a  MCHS 

The  results  suggest  that  friction  factors  for  MCs  can  be  accurately  determined  from  data  for  standard  large  channel 

Slip  flow  was  an  extremely  important  for  micro  and  in  particular  nano  heat  exchanger  designs  as  heat  transfer  rates  may 

decrease  more  than  10%  at  high  shear  rates 

Surface  roughness  in  MCs  affected  on  the  friction  factor 

As  the  supplied  heat  flux  and  vapor  quality  increased,  the  period  of  the  unstable  fluctuations  lengthened  and  the  amplitude 
also  increased 
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Table  2  ( continued ) 


Researchers 

Diameter/flow  nature 

Reynolds 

number 

Fluid/gas 

Findings 

Mokrani  et  al. 

[60] 

1  and  100  pm/laminar  and 
turbulent 

100- 

5000 

Water 

Conventional  laws  and  correlations  were  directly  applicable  to  the  MCs  of  heights  between  500  pm  and  50  pm 

Bertsch  et  al.  [61] 

1.09  and  0.54  mm 

- 

R-134a  and  R-245fa 

Heat  transfer  coefficient  for  R-245fa  in  comparison  with  R-134a  in  single-phase  flow  was  higher 

Ngo  et  al.  [62] 

1.09  mm/laminar 

250- 

2500 

co2 

Pressure-drop  factor  of  the  MCHE  with  S-shaped  fins  was  4-5  times  less  than  that  of  MCHE  with  zigzag  fins 

Wei  and  Joshi 

[63] 

Fang  et  al.  [64] 

53-112  pm/laminar 

- 

Water 

Careful  design  and  control  of  the  micro-fabrication  process  were  crucial  to  the  performance  of  the  micro-scale  devices 

100-300  pm/laminar 

_ 

Degassed  water 

Smaller  channel  was  more  favorable  due  to  its  high  mean  heat  flux 

Bogojevic  et  al. 

[65] 

Wang  et  al.  [66] 

194  |im/laminar 

82 

Water 

The  onset  of  boiling  caused  metastable  non-uniform  flow  distribution 

186  pm/laminar 

- 

Deionized  water 

Pressure  drop  was  higher  and  heat  transfer  coefficients  was  lower  for  boiling  flow  at  high  vapor  quality  in  MCs  than  those 
predicted  from  correlation  equations  for  boiling  flow  in  MCs 

Hernando  et  al. 

[67] 

100-200  pm/laminar 

10-5000 

Deionized  water 

Plate  conduction  thermal  resistance  was  a  major  restriction  for  micro-heat  exchanger  performance.  The  plate  thickness  and 
plate  material  were  critical  in  the  design  of  micro-heat  exchangers 

Yang  and  Lai  [68] 

15  mm/laminar 

300-900 

A1203/ water 

Heat  transfer  coefficient  increased  with  the  increases  of  Reynolds  number  and  particle  volume  fraction 

Wu  et  al.  [69] 

L/D= 50-150/laminar 

- 

R134a,  R123,  R245fa, 

R236fa,  nitrogen  and  water 

This  proposed  method  almost  predicts  97%  of  the  non-aqueous  data  and  94%  of  the  water  data  within  the  +  30%  error  band 

Shevade  and 
Rahman  [70] 

0.06-0.2  cm/turbulent 

1600- 

3000 

Water 

The  pressure  drop  in  the  channel  increased  as  the  Reynolds  number  increases.  Nusselt  number  increased  as  the  channel 
dimensions  increase 

Mlcak  et  al.  [72] 

85.6  pm/laminar 

50-400 

Water 

Friction  coefficient  linearly  increased  as  the  hydraulic  entrance  length  increase 

Chen  et  al.  [73] 

40-158  pm/laminar 

30-500 

Deionized  water 

Triangular  MC  heat  sink  has  the  highest  thermal  efficiency  among  the  three  heat  sinks,  the  trapezoidal  MC  was  inferior  to 
triangular  but  superior  to  rectangular 

Zhu  and  Liao  [74] 

-/laminar 

- 

Gas 

For  a  rectangular  and  triangular  microchannel,  the  average  Nusselt  numbers  were  smaller  than  those  of  a  macrochannel  with 
the  same  thermal  boundary  conditions  and  it  decreases  with  increasing  Knudsen  number. 

Ji  et  al.  [75] 

Length  1-100  |im/laminar 

0.001- 

100 

Gas 

Average  Poiseuille  number  increased  not  only  as  the  roughness  height  increases  but  also  as  the  spacing  distance  between  the 
roughness  elements  was  decreased 

Sui  et  al.  [76] 

1.5  x  10~4  m 

100-800 

Water 

The  quantity  and  the  location  of  the  vortices  may  change  along  the  flow  direction,  leading  to  chaotic  advection,  which  can 
greatly  enhance  the  convective  fluid  mixing,  and  thus  the  heat  transfer  performance  of  the  wavy  MCs  was  much  better  than 
that  of  straight  MCs  with  the  same  cross  section 

Muhammad 
et  al.  [77] 

Height=300  [im.  Depth  =  100- 
250  pm 

1000  to 

3000 

FC-72,  helium  and  water 

Compared  to  helium  and  FC-72,  water  provides  better  thermal  performance  with  higher  values  of  Nusselt  number 

Liu  et  al.  [78] 

100  pm/laminar 

10-200 

Water 

Heat  transfer  enhancement  due  to  viscosity-variation  was  pronounced,  though  the  axial  conduction  introduced  by  thermal- 
conductivity-variation  was  insignificant  unless  for  the  cases  with  very  low  Reynolds  numbers 

Renksizbulut 

Width/height =0-1  /laminar 

0.1-10 

Gas 

Slip-flow  regime  was  associated  with  very  large  reductions  in  wall  friction  and  heat  transfer  in  the  entrance  region 

et  al.  [79] 

Niazmand  et  al. 
[80] 

Width/height= 0.25-2/laminar 

0.1-10 

Gas 

Major  reductions  in  the  friction  and  heat  transfer  coefficients  were  observed  in  the  entrance  region  due  to  large  amounts  of 
velocity-slip  and  temperature-jump 

Khadem  et  al. 

[81] 

Zade  et  al.  [82] 

L/H =20  laminar 

- 

Gas 

Roughness  has  more  significant  effect  on  higher  Knudsen  number  flows  with  higher  relative  roughness. 

10  pm/laminar 

0.1 

Gas 

In  the  entrance  region,  Nusselt  number  was  only  affected  by  the  change  in  the  Knudsen  number  and  was  higher  in  variable- 
density  slip-flow  cases. 

Shokouhmand 

- 

2 

Gas 

Nusselt  number  and  Q.  Re  decreased  with  increaseing  Knudsen  number 

and  Bigham 

[83] 

Mohammed 
et  al.  [84,85] 

339.15  pm/laminar 

100-1000 

Water 

Temperature  and  the  heat  transfer  coefficient  of  the  zigzag  MCHS  was  the  least  and  greatest  comparing  to  other  channel 
shapes.  The  second  best  thermal  performance  was  provided  by  curvy  channels  after  the  wavy  channels 

Kosar [86] 

Height=200  pm  Width  =  200  pm/ 
Laminar 

50-2000 

Water 

The  Nusselt  number  was  lower  for  low  thermal  conductivity  and  vice  versa 

Mohammed 

et  al.  [87,88] 

100  to 

800 

Al203/water 

Si02/water 

Ag/  water 

Ti02/water 

A1203  nanofluid  has  the  highest  heat  transfer  coefficient 
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that  the  observed  enhancement  of  the  thermal  conductivity  of 
Al203/water  nanofluids  exists  and  is  independent  of  the  measure¬ 
ment  technique. 

Zhang  et  al.  [95]  reported  measurements  of  the  effective  thermal 
conductivity  and  thermal  diffusivity  of  Al203/water,  Zr02/water, 
Ti02/water,  and  CuO/water  nanofluids  for  various  volume  fractions 
and  temperature  by  using  the  transient  short-hot-wire  technique. 
The  average  diameters  of  Al203t  Zr02,  Ti02,  and  CuO  particles  were 
20  nm,  20  nm,  40  nm  and  33  nm,  respectively.  The  results  demon¬ 
strated  that  the  effective  thermal  diffusivity  and  thermal  conduc¬ 
tivity  increase  with  an  increase  in  the  particle  concentration  and  in 
the  thermal  conductivity  of  nanoparticles.  However,  the  effective 
thermal  conductivity  did  not  show  any  anomalous  enhancement. 

Murshed  at  al.  [96]  presented  investigation  of  enhancing  the 
thermal  conductivity  of  Ti02-water  based  nanofluids.  The  nano¬ 
fluids  were  prepared  by  dispersing  Ti02  nanoparticles  in  rod- 
shapes  of  <p= 10  nm  x  40  nm  (diameter  by  length)  and  in  spherical 
shapes  of  $=  15  nm  in  deionized  water.  The  thermal  conductivities 
of  these  nanofluids  were  measured  by  using  the  transient  hot-wire 
apparatus  with  an  integrated  correlation  model.  The  viscosity  and 
pH  value  of  the  nanofluids  were  also  characterized.  The  experi¬ 
mental  results  showed  that  the  thermal  conductivity  increases 
with  an  increase  of  particle  volume  fraction  and  the  nanofluids 
containing  a  small  amount  of  nanoparticles  have  much  higher 
thermal  conductivities  than  normal  base  fluids.  The  particle  size 
and  shape  also  have  influence  of  the  thermal  conductivity  enhance¬ 
ment  of  nanofluids.  The  results  of  Ti02  particles  of  ^10  nm  x  40  nm 
and  ^15  nm  dimensions  with  maximum  5%  volume  fraction 
showed  that  the  enhancement  was  nearly  33%  and  close  to  30%, 
respectively  over  the  base  fluid. 

Murshed  et  al.  [97]  conducted  an  experimental  and  theoretical 
analyzing  on  the  effective  thermal  conductivity  and  viscosity  of 
nanofluids.  Two  static  mechanisms-based  models  were  presented 
to  predict  the  enhanced  thermal  conductivity  of  nanofluids  having 
spherical  and  cylindrical  nanoparticles.  The  transient  hot-wire 
technique  was  utilized  to  measure  the  thermal  conductivity  of 
nanofluids  at  different  temperatures  ranging  from  20  °C  to  60  C. 
The  results  of  this  investigation  showed  that  the  thermal  con¬ 
ductivity  and  viscosity  of  nanofluids  were  found  to  be  substan¬ 
tially  higher  than  the  values  of  the  base  fluids.  Both  the  thermal 
conductivity  and  viscosity  of  nanofluids  increase  with  increases 
the  nanoparticle  volume  fraction.  The  results  also  presented  that 
the  thermal  conductivity  of  nanofluids  observed  to  be  strongly 
dependent  on  temperature  and  increases  linearly  with  the  increases 
of  temperature.  The  results  of  the  proposed  models  which  consider 
particle  size,  interfacial  layer,  and  volume  fraction  showed  reason¬ 
ably  good  agreement  with  the  experimental  results  and  give  better 
predictions  for  the  effective  thermal  conductivity  of  nanofluids 
compared  to  existing  classical  models. 

Oh  et  al.  [98]  presented  experimental  data  of  the  thermal 
conductivity  enhancement  in  A1203  nanofluids  with  D1  water  and 
ethylene  glycol  as  base  fluids  by  using  the  modified  3a>  method. 
The  results  revealed  that  the  thermal  conductivity  of  the  nanofluid 
continuously  increase  with  time.  The  effect  of  gravitational  on  the 
nanoparticle  aggregates  is  believed  to  be  the  reason  behind  this 
observation.  The  increscent  of  thermal  conductivity  due  to  sedi¬ 
mentation  effect  was  also  measured  quantitatively.  In  comparison 
with  the  conventional  thermal  conductivity  measurement  meth¬ 
ods,  the  proposed  3 co  measurement  technique  offered  several 
advantages  and  it  can  be  used  with  only  a  small  sample  volume 
to  show  whether  homogeneous  mixing  is  achieved  for  the 
stationary  nanofluid. 

Petal  et  al.  [99]  presented  comprehensive  experimental  inves¬ 
tigation  of  thermal  conductivity  of  nanofluids  with  variation  in 
nanoparticle  material,  base  liquid,  particle  size,  particle  volume 
fraction  and  suspension  temperature.  Transient  hot  wire  (THW) 


equipment  as  well  as  temperature  oscillation  equipment  were 
developed  for  the  measurement  of  thermal  conductivity  of  liquids. 
The  results  showed  that,  in  general,  thermal  conductivity  values  of 
all  the  nanofluids  are  higher  than  that  of  the  equivalent  macro¬ 
particle  suspensions.  A  very  significant  effect  of  particle  size  and 
suspension  temperature  was  observed  for  the  thermal  conductiv¬ 
ity  of  nanoparticle  suspensions.  The  results  also  revealed  to  that 
the  metallic  nanofluids  were  found  to  give  higher  enhancements 
than  that  of  oxide  nanofluids.  The  investigation  results  also 
presented  that  the  increasing  in  the  temperature  significantly 
increase  the  thermal  conductivity  of  nanofluid.  It  was  also  found 
that  the  thermal  conductivity  of  nanoparticle  suspensions  is 
relatively  higher  at  lower  volume  fractions,  thereby  giving  a 
non-linear  dependence  on  the  particle  volume  fraction.  However, 
the  enhancement  in  thermal  conductivity  of  suspensions  increases 
with  the  decrease  in  the  base  fluid  conductivity.  This  is  due  to  the 
fact  that  the  relative  contribution  of  particles  in  the  total  heat 
conduction  increases  with  reduced  contribution  from  the  lower 
conductivity  base  fluids. 

Xie  et  al.  [100]  investigated  the  thermal  conductivity  behavior 
of  SiC  dilute  nanoparticle-fluid  mixtures.  The  transient  hot-wire 
method  was  utilized  to  measure  the  thermal  conductivities  of  four 
series  suspensions.  The  experimental  results  showed  that  the 
thermal  conductivities  of  SiC  suspensions,  containing  a  small 
amount  of  solid  particles,  are  significantly  higher  than  those  of 
the  base  liquids.  For  low  volume  fraction  range,  the  thermal 
conductivity  enhancing  ratios  increasing  almost  linearly  with  the 
increase  of  the  volume  fraction  of  solid  particles.  The  results  of 
different  systems  containing  the  same  kind  of  particles  showed 
that  the  thermal  conductivity  ratios  were  independent  of  the  base 
fluids,  but  the  absolute  value  of  the  thermal  conductivity  of  the 
suspension  is  proportional  to  that  of  the  base  liquids. 

Sharma  et  al.  [101]  presented  nanofluids  synthesized  using 
silver  nitrate  (precursor),  ethylene  glycol  (reducing  agent),  and 
polyacrylamide-co-acrylicacid  (dispersion  stabilizer)  in  order  to 
enhance  the  thermal  conductivity.  The  different  concentrations  of 
silver  nanofluid  (1000-10,000  ppm)  were  synthesized.  The  thermal 
conductivity  as  well  as  stability  of  these  nanofluids  were  deter¬ 
mined  with  a  transient  hot-wire  apparatus.  Typically,  10,000  ppm 
silver  nanofluid  exhibited  rapid  increase  in  the  particle  size  with  the 
passage  of  time.  The  results  presented  that  the  thermal  conductivity 
of  silver  nanofluids  increased  with  increases  the  amount  of  silver 
particles  in  nanofluid.  The  results  showed  that  after  30  days  of 
preparation,  the  thermal  conductivity  of  1000  ppm  and  5000  ppm 
silver  nanofluids  slightly  decreased  from  10%,  and  15%  to  9%,  and 
14%,  respectively.  While,  the  thermal  conductivity  of  10,000  ppm 
nanofluid  decreased  after  30  days  of  preparation  from  18%  to  14%.  It 
was  observed  also  that  the  silver  particles  were  aggregated  in  early 
stage  of  preparation  (upto  15  days),  which  leads  to  increase  in  the 
size  of  silver  particles.  However,  the  stability  of  silver  nanofluid  was 
strongly  affected  by  the  characteristics  of  the  suspended  particle 
and  base  fluid  such  as  the  particle  morphology,  the  chemical 
structure  of  the  particles  and  base  fluid. 

Wong  and  Kurma  [102]  presented  study  of  three  transport 
properties,  namely  thermal  conductivity,  electrical  conductivity  and 
viscosity,  for  Alumina  nanoparticles  with  a  mean  diameter  of  36  nm 
dispersed  in  water.  Experiments  were  performed  both  as  a  function 
of  volumetric  concentration  (3-8%)  and  temperature  (2-50  °C).  The 
results  illustrated  that  the  thermal  conductivity  of  alumina  nanofluid 
increased  with  both  increase  in  temperature  and  concentration. 
A  maximum  thermal  conductivity  of  0.7351  W/m  I<  1  was  observed 
for  an  8.47%  volume  concentration  of  alumina  nanoparticles  at 
46.6  C.  The  effective  thermal  conductivity  at  this  concentration 
and  temperature  was  1.1501,  which  lead  to  an  increase  in  thermal 
conductivity  by  22%  when  compared  to  water  at  room  temperature. 
In  order  to  determine  the  kinematic  viscosity  of  alumina  nanofluid,  a 
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standard  kinematic  viscometer  with  constant  temperature  bath  was 
used.  The  results  revealed  to  an  increase  of  35.5%  in  the  kinematic 
viscosity  for  an  8.47%  volumetric  concentration  of  alumina  nanopar¬ 
ticles  in  water.  The  maximum  kinematic  viscosity  of  alumina  nano¬ 
fluid  was  2.901  42  mm2/  s,  obtained  at  0  °C  for  an  8.47%  volumetric 
concentration  of  alumina  nanoparticles. 

Turgut  et  al.  [103]  presented  investigation  of  the  thermal 
conductivity  and  viscosity  of  Ti02  nanoparticles  in  deionized  water 
for  volume  fraction  from  0.2%  to  3%  of  particles.  The  mean 
diameter  of  Ti02  nanoparticles  was  21  nm.  While,  the  measure¬ 
ment  of  nanofluids  thermal  conductivity  was  done  by  using  the  3m 
method.  The  effective  thermal  conductivity  of  Ti02  nanoparticles 
in  deionized  water  was  measured  at  temperatures  of  13  °C,  23  °C, 
40  “C,  and  55  °C.  The  experimental  results  demonstrated  that  the 
thermal  conductivity  increases  with  an  increase  of  particle  volume 
fraction,  and  the  enhancement  was  observed  to  be  7.4%  over  the 
base  fluid  for  a  nanofluid  with  3%  volume  fraction  of  Ti02 
nanoparticles  at  13  °C.  The  results  also  showed  that  there  is  no 
dependence  related  to  temperature;  the  thermal  conductivity 
increased  by  the  same  order  of  magnitude  as  the  base  fluid  which 
is  water.  The  effective  viscosities  of  Ti02-water  nanofluids  with 
concentrations  from  0.2  vol%  to  3.0  vol%  were  measured.  The 
results  revealed  that  for  low  volume  concentration  of  nanoparti¬ 
cles  for  0.2  vol%.  fraction,  the  viscosity  values  follow  quite  well  the 
viscosity  values  of  pure  water  with  a  decrease  in  viscosity  with 
increasing  temperature.  While,  for  higher  additions  of  Ti02  nano¬ 
particles,  the  increase  in  the  nanofluid  viscosity  was  larger  than 
the  enhancement  in  the  thermal  conductivity. 

Chandrasekar  and  Suresh  [104]  presented  development  of  lower 
and  upper  limits  for  thermal  conductivity  of  nanofluids.  The  upper 
limit  is  estimated  by  coupling  heat  transfer  mechanisms  like 
particle  shape,  nanolayer  and  Brownian  motion.  While  the  lower 
limit  is  based  on  Maxwell's  equation.  Experimental  data  from  a 
range  of  independent  published  sources  is  used  for  validation  of  the 
developed  limits.  The  comparison  with  the  published  experimental 
data  indicated  that,  the  experimental  data  considered  lie  between 
the  newly  developed  limits.  Comparison  also  showed  that  the  present 
limits  are  more  rigorous  in  placing  a  narrow  lower  and  upper  limit.  As 
most  of  the  experimental  data  lies  within  the  newly  developed  limits, 
it  can  be  concluded  that  particle  shape,  nanolayer  and  Brownian 
motion  are  significant  in  enhancing  the  thermal  conductivity  of 
nanofluids. 

Chen  et  al.  [105]  proposed  a  methodology  for  predicting  the 
effective  thermal  conductivity  of  dilute  suspensions  of  nanoparti¬ 
cles  (nanofluids)  based  on  rheology.  The  methodology  utilizes  the 
rheological  data  to  infer  microstructures  of  nanoparticles  quanti¬ 
tatively,  which  is  then  included  into  the  conventional  Hamilton- 
Crosser  equation  to  predict  the  effective  thermal  conductivity  of 
nanofluids.  The  results  showed  that  the  conventional  correlations 
failed  to  predict  the  thermal  conductivity  of  nanofluids,  while,  by 
incorporating  the  microstructure  parameters  acquired  from  rheo¬ 
logical  measurements,  the  thermal  conductivity  can  be  predicted 
accurately.  The  results  also  presented  the  importance  of  the 
rheological  characterization  of  nanofluids,  which  could  provide 
structural  information  required  for  the  prediction  of  the  suspen¬ 
sion  thermal  conductivity. 

Hong  et  al.  [106]  proposed  a  mixture  of  nanoparticles  and  fluid, 
have  enormous  potential  to  improve  the  efficiency  of  heat  transfer 
fluids.  Fe  nanofluids  are  prepared  with  ethylene  glycol  and  Fe 
nanocrystalline  powder  synthesized  by  a  chemical  vapor  conden¬ 
sation  process.  The  results  revealed  to  that  the  thermal  conduc¬ 
tivity  of  a  Fe  nanofluid  is  increased  nonlinearly  up  to  18%  as  the 
volume  fraction  of  particles  is  increased  up  to  0.55  vol%.  In 
comparison  with  the  Cu  nanofluids  dispersed  with  a  little  agglom¬ 
eration,  Fe  nanofluids  showed  a  more  effective  thermal  transport 
property.  Hence,  the  highly  thermally  conductive  material  is  not 


always  the  best  candidate  for  the  suspension  in  improving  thermal 
transport  property  of  fluids. 

Liu  et  al.  [107]  presented  study  of  enhancement  of  the  thermal 
conductivity  of  water  in  the  presence  of  copper  (Cu)  using  the 
chemical  reduction  method.  Copper  nanoparticles  showed  the  char¬ 
acteristics  of  small  particle  sizes  and  uniform  size  distribution.  The 
typical  copper  nanoparticles  diameter  were  around  50-100  nm  with 
square  and  spherical  shapes.  The  volume  fractions  of  Cu  nanoparticle 
suspensions  in  water  liquid  were  in  the  range  from  0.05  vol%  to 
0.2  vol%.  The  results  showed  that  without  the  addition  of  dispersant 
and  surfactant,  the  thermal  conductivity  of  the  produced  nanofluids 
reveals  a  time-dependent  characteristic.  The  results  revealed  to  that, 
Cu  water  nanofluids  with  low  concentration  of  nanoparticles  have 
noticeably  higher  thermal  conductivities  than  the  water  base  fluid 
without  Cu.  Furthermore,  for  Cu  nanoparticles  at  a  volume  fraction  of 
0.1  vol%,  thermal  conductivity  was  enhanced  by  up  to  23.8%. 

Yu  et  al.  [108]  presented  investigation  of  the  thermal  conductivity 
and  viscosity  of  Ethylene  glycol  (EG)  based  nanofluids  containing  ZnO 
nanoparticles.  The  results  showed  that  the  thermal  conductivity  of 
ZnO-EG  nanofluids  is  independent  of  setting  time  from  20  min  to 
360  min.  However,  the  absolute  thermal  conductivity  increases  with 
the  temperature  range  from  10  °C  to  60  °C,  while  the  enhanced  ratios 
were  almost  constant.  The  thermal  conductivity  of  ZnO-EG  nanofluids 
was  strongly  depends  on  the  particle  concentration  and  it  increases 
nonlinearly  with  the  volume  fraction  of  nanoparticles.  The  enhanced 
in  thermal  conductivity  of  ZnO-EG  nanofluid  was  26.5  vol%  of 
5.0  vol%.  The  rheological  behaviors  of  the  nanofluids  showed  that 
ZnO-EG  nanofluids  demonstrated  Newtonian  behaviors  with  low 
volume  concentrations.  While,  the  shear-shinning  behavior  observed 
for  higher  volume  concentrations  nanofluids  and  showing  non- 
Newtonain  behavior  due  to  the  effective  volume  fraction  of  aggregates 
is  much  higher  than  the  actual  solid  volume  fraction. 

Xie  et  al.  [109]  presented  a  technical  route  for  preparing  homo¬ 
geneous  and  stable  nanofluids  including  diamond  nanoparticles 
(DNPs)  with  the  mixture  of  distilled  water  (DW)  and  ethylene  glycol 
(EG)  (DW-EG)  as  the  base  fluid.  The  results  showed  that  the  thermal 
conductivity  enhancement  of  the  DNP  nanofluids  increases  with 
increase  the  volume  fraction  of  DNPs.  The  thermal  conductivity 
enhancement  was  more  than  18.0%  for  a  nanofluid  at  a  DNP  volume 
fraction  of  0.02.  The  convective  heat  transfer  coefficient  increases 
with  increasing  volume  fraction  of  DNPs,  first,  and  then  it  decreases 
with  a  further  increase  in  the  volume  fraction  of  DNPs. 

Jimenez-Perez  et  al.  [110]  proposed  a  study  of  using  thermal 
lens  spectrometry  (TLS)  and  photopyroelectric  (PPE)  techniques  to 
obtain  the  thermal  diffusivity  and  effusivity  of  three  different 
solvents  samples  ethanol,  water,  and  ethylene  glycol  (EG)  with 
gold  nanoparticles.  The  thermal  effusivity  of  these  samples  was 
determined  by  the  PPE  technique  in  a  front  detection  configura¬ 
tion.  The  experimental  results  revealed  that  the  thermal  diffusivity 
values  of  different  solvents  (water,  ethanol,  and  ethylene  glycol) 
were  the  presence  of  gold  nanoparticles  increases  when  compared 
with  the  pure  solvents. 


4.2.2.  Effects  of  thermal  conductivity  on  heat  transfer 

Many  experimental  studies  have  been  done  by  different 
researchers  for  examining  and  improving  the  thermal  conductivity 
enhancement  of  certain  types  of  nanofluids  using  different  meth¬ 
ods  and  parameters.  They  calculated  the  thermal  conductivity 
enhancement  by  using  the  ratio  of  the  thermal  conductivity  of  the 
nanofluid  over  its  base  fluid  thermal  conductivity.  It  is  noticed  that 
most  of  the  studies  did  not  concerned  or  noticed  the  particle  size 
distribution.  In  addition,  the  effect  of  the  particles  agglomeration 
in  the  base  fluid  is  also  unknown.  It  is  declared  that  the  thermal 
conductivity  of  nanofluid  is  sensitive  to  certain  parameters  which 
are  responsible  for  constant  nanofluids  such  as  particle  size, 
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particle  material,  particle  shape,  base  fluid  material,  temperature, 
additives  and  surfactants,  and  acidity.  These  parameters  will  be 
introduced  in  the  following  sections  separately  base  on  the  previous 
experimental  investigations  and  findings. 

4.2.2.).  Effects  of  particle  size  on  heat  transfer.  The  size  of  particle 
plays  a  significant  role  in  thermal  conductivity  and  heat  transfer 
enhancement  in  base  fluids.  As  introduced  before  that  different 
suspension  sizes  of  millimeter  and  micrometer  particles  in  host 
liquid  have  been  used  to  increase  the  thermal  conductivity,  but 
the  particles  agglomerate  quickly  and  settled  out  of  the  liquid. 
Nanoparticles  are  formed  to  decrease  the  size  of  the  particles  so 
that  the  time  of  sedimentation  decreases  to  reach  in  some  cases  to 
more  than  few  weeks  or  months. 

Petal  et  al.  [Ill]  presented  comprehensive  experimental  inves¬ 
tigation  of  thermal  conductivity  of  nanofluids  with  variation  in 
nanoparticle  material,  base  liquid,  particle  size,  particle  volume 
fraction  and  suspension  temperature.  Transient  hot  wire  (THW) 
equipment  as  well  as  temperature  oscillation  equipment  were 
developed  for  the  measurement  of  thermal  conductivity  of  liquids. 
The  results  showed  that,  in  general,  thermal  conductivity  values  of 
all  the  nanofluids  are  higher  than  that  of  the  equivalent  macro¬ 
particle  suspensions.  A  very  significant  effect  of  particle  size  and 
suspension  temperature  was  observed  for  the  thermal  conductiv¬ 
ity  of  nanoparticle  suspensions.  The  results  also  revealed  to  that 
the  metallic  nanofluids  were  found  to  give  higher  enhancements 
than  that  of  oxide  nanofluids.  The  particle  size  was  found  to  have  a 
tremendous  effect  on  the  thermal  conductivity  of  nanofluids  with 
enhancement  in  the  thermal  conductivity  increasing  with  the 
reduction  in  the  particle  size.  The  effect  of  particle  size  may  be 
attributed  mainly  to  two  reasons  which  are  the  high  specific 
surface  area  of  the  nanoparticles  and  Brownian  motion. 

Timofeeva  et  al.  [112]  investigated  the  effect  of  average  particle 
sizes  on  basic  macroscopic  properties  and  heat  transfer  performance 
of  water-based  silicon  carbide  suspensions.  The  average  particle  sizes, 
calculated  from  the  specific  surface  area  of  nanoparticles,  were 
varied  from  16  nm  to  90  nm.  The  results  revealed  that  the  nanofluids 
of  same  material  and  volume  concentration  with  larger  particles 
provide  higher  thermal  conductivity  and  lower  viscosity  increases 
than  those  with  smaller  particles  due  to  the  smaller  solid/liquid 
interfacial  area  of  larger  particles. 

Anoop  et  al.  [113]  presented  an  experimental  investigation  on 
the  convective  heat  transfer  characteristics  in  the  developing 
region  of  tube  flow  with  constant  heat  flux  carried  out  with 
alumina-water  nanofluids.  In  this  study  the  effect  of  particle  size 
on  convective  heat  transfer  in  laminar  developing  region  was 
evaluated.  Two  particle  sizes  were  used,  one  with  average  particle 
size  off  45  nm  and  the  other  with  150  nm.  The  results  revealed 
that  both  nanofluids  showed  higher  heat  transfer  characteristics 
than  the  base  fluid.  Nanofluids  with  of  45  nm  particles  showed 
higher  heat  transfer  coefficient  than  that  with  150  nm  particles. 
The  results  also  showed  with  increase  in  flow  rate  and  particle 
concentration  the  average  heat  transfer  coefficient  value  was 
increased. 

Beck  et  al.  [114]  presented  a  new  data  for  the  thermal  conduc¬ 
tivity  enhancement  in  seven  nanofluids  containing  8-282  nm  dia¬ 
meter  alumina  nanoparticles  in  water  or  ethylene  glycol.  The  results 
showed  that  the  thermal  conductivity  enhancement  in  these  nano¬ 
fluids  decreases  as  the  particle  size  decreases  below  about  50  nm. 
This  decrease  in  enhancement  can  be  attributed  to  a  decrease  in  the 
thermal  conductivity  of  the  nanoparticles  themselves  as  the  particle 
size  becomes  small  enough  to  be  affected  by  increased  phonon 
scattering. 

Chopkar  et  al.  [115]  presented  experimental  investigation  of  the 
impact  of  Al2Cu  and  Ag2Al  nanoparticle  size  and  volume  fraction  on 


the  effective  thermal  conductivity  of  water  and  ethylene  glycol  based 
nanofluid  prepared  by  a  two-stage  process.  The  thermal  conductivity 
ratio  of  nanofluid,  was  measured  by  using  an  indigenously  developed 
thermal  comparator  device.  The  results  revealed  that  the  thermal 
conductivity  ratio,  relative  to  that  of  base  fluid,  increases  nonlinearly 
with  the  increase  in  volume  fraction  and  the  decrease  in  the  size/ 
diameter  of  nanoparticles.  This  increase  was  up  to  100  pet  with  only 
1.5  vol%/pct  nanoparticles  of  30-  to  40-nm  average  diameter.  More¬ 
over,  the  analytical  model  denoted  that  the  interfacial  layer  signifi¬ 
cantly  influences  the  effective  thermal  conductivity  ratio  of  nanofluid 
for  the  comparable  amount  of  nanoparticles. 

Ding  and  Wen  [116]  presented  investigation  of  migration 
particle  migration  in  pressure-driven  laminar  pipe  flows  of  rela¬ 
tively  dilute  suspensions  of  nanoparticles  (nanofluids).  A  theore¬ 
tical  model  is  formulated  to  predict  particle  concentration,  and 
velocity  field  of  nanofluids  in  the  transverse  plane  of  the  pipe. 
The  results  showed  that  the  particle  concentration  in  the  wall 
region  can  be  much  lower  than  that  in  the  central  core  region. 
That  indicates  a  highly  nonuniform  thermal  conductivity  profile 
across  the  transverse  plane  of  the  pipe,  and  thus  has  a  significant 
implication  to  heat  transfer  intensification  using  nanofluids. 

4.2.22.  Effects  of  acidity  (PH).  Limited  studies  have  been  published 
on  the  effect  of  fluid  acidity  on  the  thermal  conductivity  enhance¬ 
ment  of  nanofluids,  but  the  general  trend  is  that  acidity  increases 
the  thermal  conductivity  enhancement.  Xian-Ju  and  Xian-Fang  [117] 
investigated  the  viscosity  and  the  thermal  conductivity  of  Cu  and 
A1203  nanoparticles  dispersed  in  water  under  different  pH  values. 
The  experimental  results  showed  that  the  viscosity  is  closely  related 
to  the  enhancement  of  thermal  conductivity.  Furthermore,  there 
exists  an  optimal  pH  value  can  results  in  the  highest  thermal 
conductivity  and  the  lowest  viscosity  of  nanofluids. 

Wang  et  al.  [118]  investigated  the  enhancement  of  heat  transfer 
of  Al203  and  Cu  nanoparticles  dispersed  in  water  as  a  base  fluid 
under  different  pH  values  and  different  sodium  dodecylbenzene 
sulfonate  (SDBS)  dispersant  concentrations.  The  experimental 
results  showed  that  the  stability  of  nanofluids  has  a  good  corre¬ 
sponding  relation  with  thermal  conductivity  and  the  better  dis¬ 
persion  behavior,  the  higher  thermal  conductivity. 

4.2.23.  Effects  of  particles  material.  There  are  many  types  of 
materials  have  been  utilized  to  produce  nanofluids.  It  was  found 
that  metallic  and  allotrops  of  carbon  nanoparticles  (such  as  diamond) 
increase  the  thermal  conductivity  of  a  host  fluid  more  than  other  types 
of  nanoparticles  having  the  same  volume  concentration.  Nguyen  et  al. 
[119]  investigated  experimentally  the  behavior  and  heat  transfer 
enhancement  of  a  particular  nanofluid,  A1203  nanoparticle-water 
mixture,  flowing  inside  a  closed  system  that  is  destined  for  cooling 
of  microprocessors  or  other  electronic  components.  The  heat  transfer 
coefficient  of  a  particular  nanofluid  with  6.8%  particle  volume 
concentration,  was  found  to  increase  as  much  as  40%  compared  to 
that  of  the  base  fluid,  and  the  results  clearly  revealed  that  the  heat 
transfer  coefficient  was  considerably  increase  with  an  increase  of 
particle  concentration.  However,  an  increase  of  particle  concentration 
produced  a  clear  decrease  of  the  heated  component  temperature. 

4.22.4.  Effects  of  temperature.  The  nature  of  the  nanofluid  thermal 
conductivity  is  sensitive  to  temperature.  Many  studies  have  been 
done  for  the  effect  of  temperature  on  thermal  conductivity.  Minsta 
et  al.  [120]  presented  experimental  determination  of  thermal 
conductivity  values  for  three  types  of  water-based  nanofluids  for 
various  temperatures  ranging  between  20  °C  and  40  °C.  The  consi¬ 
dered  nanofluids  were  composed  of  29  nm  CuO  particles  as  well  as 
36  nm  and  47  nm  A1203  particles.  For  each  nanofluid  considered, 
three  volume  fractions  of  particles  were  investigated:  3%,  6%  and  9% 
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between  20  °C  and  40  °C.  The  thermal  conductivity  measurement 
technique  utilized  was  based  on  the  transient  hot  wire  method. 
The  results  clearly  demonstrated  that  the  effective  thermal 
conductivity  of  nanofluids  increases  with  temperature.  In  the 
temperature  range  of  20  "C  and  40  °C,  the  average  an  increase  in 
thermal  conductivity  was  approximately  16%  for  each  type  nanofluid. 

Yu-Hua  et  al.  [121]  studied  the  temperature  effects  on  the  thermal 
conductivity  of  nanofluids.  The  mechanism  of  thermal  conductivity  of 
nanofluids  was  analyzed  and  calculated,  including  particle  agglom¬ 
eration,  Brownian  motion  effects  and  viscosity,  together  influenced  by 
temperature.  The  results  showed  that  the  higher  temperature  lead  to 
decrease  the  of  agglomeration  of  nanoparticles  through  the  reduction 
of  the  particle  surface  area,  and  makes  the  Brownian  motion  more 
intensive  due  to  the  decreases  of  viscosity. 

Nguyen  et  al.  [122]  investigated  experimentally  the  influence  of 
both  the  temperature  and  the  particle  size  on  the  dynamic  viscosities 
of  two  particular  water-based  nanofluids,  namely  water-Al203  and 
water-CuO  mixtures.  The  results  revealed  that  for  particle  volume 
fractions  less  than  4%,  viscosities  of  36  nm  and  47  nm  particle-size 
alumina-water  nanofluids  are  approximately  identical.  However,  the 
viscosities  of  47  nm  particle-size  and  higher  particle  fractions  were 
clearly  higher  than  those  of  36  nm  size.  The  results  showed  that  for 
all  the  nanofluids  tested,  their  viscosities  were  strongly  dependent  on 
both  particle  volume  fraction  and  temperature. 

Hosseini  and  Ahmadi  [123]  presented  study  of  the  structural 
changes  of  nanofluid  as  increasing  the  temperature  which  is  very 
important  for  engineering  designs.  The  results  revealed  that  the 
nanoparticles  increase  the  temperature  of  nanofluid  when  it  is 
prepared  in  ultrasonic  mixer  than  other  impure  fluids.  The  main 
reason  of  increasing  the  temperature  of  operating  suspension  than 
solution  mainly  due  to  sending  out  the  waves  from  ultrasonic  prop, 
collapsing  between  particles,  and  touching  fraction  between  particles. 

42.2.5.  Effects  of  base  fluids  material.  Many  host  liquids  have  been 
used  to  produce  nanofluids  such  as  aqueous  and  organic  liquids 
(ethylene  glycol,  and  oils)  to  enhance  thermal  conductivity.  Xie 
et  al.  [124]  analyzed  the  effect  of  the  base  fluid  on  the  thermal 
conductivity  of  nanofluids.  Nanofluids  with  A1203  nanoparticles 
were  prepared  by  using  different  base  fluids;  deionized  water, 
ethylene  glycol,  glycerol,  and  pump  oil.  In  addition,  glycerol-water 
and  ethylene  glycol-water  mixtures  with  different  volume 
fractions  were  also  used  as  base  fluids  and  the  variation  of  the 
thermal  conductivity  ratio  with  thermal  conductivity  of  the  base 
fluid  mixture  was  examined.  The  results  revealed  that,  the  thermal 
conductivity  ratio  decreased  with  increasing  thermal  conductivity 
of  the  base  fluid. 

Li  and  Peterson  [125]  presented  investigation  of  analyzing, 
modeling  for  the  mixing  effect  of  the  base  fluid  in  the  immediate 
vicinity  of  the  nanoparticles  caused  by  the  Brownian.  The  results 
indicated  that  the  micro  convection/mixing  induced  by  the  Brow¬ 
nian  motion  of  the  nanoparticles  could  significantly  affect  the 
macro  heat  transfer  capability  of  the  nanofluids. 

Xie  et  al.  [126]  investigated  the  thermal  conductivities  of  various 
suspensions  containing  A1203  nanoparticles  with  specific  surface 
areas  prepared  in  a  range  of  5-124  m2  / g.  The  results  illustrated  that 
the  nanoparticle  suspensions,  containing  a  small  amount  of  A1203, 
have  substantially  higher  thermal  conductivity  than  the  base  fluid. 
The  enhancement  of  thermal  conductivity  increases  with  the 
volume  fraction  of  A1203.  Whereas,  the  enhanced  thermal  conduc¬ 
tivity  increases  with  an  increase  in  the  difference  between  the  pH 
value  of  aqueous  suspension  and  the  isoelectric  point  of  A1203 
particle. 

Maiga  et  al.  [127]  presented  a  study  of  forced  convection  flow 
of  water-yAl203  and  ethylene  glycol-yAl203  nanofluids  inside  a 
uniformly  heated  tube  that  is  subjected  to  a  constant  and  uniform 


heat  flux  at  the  wall.  The  results  showed  that  the  inclusion  of 
nanoparticles  has  increased  considerably  the  heat  transfer  at  the 
tube  wall  for  both  the  laminar  and  turbulent  regimes.  However, 
with  increasing  the  particle  concentration  the  improvement  of 
heat  transfer  becomes  more  pronounced.  The  results  also  demon¬ 
strated  that  the  ethylene  glycol-yAl203  mixture  gives  a  far  better 
heat  transfer  enhancement  than  the  water-yAl203  mixture. 


4.22.6.  Specific  heat.  Shin  et  al.  [128]  presented  investigation  of 
using  Silica  nanoparticles  (1%  by  weight)  dispersed  in  a  eutectic  of 
lithium  carbonate  and  potassium  carbonate  (62:38  ratio)  to  obtain 
enhancement  of  specified  heat.  A  differential  scanning  calorimeter 
instrument  was  used  to  measure  the  specific  heat  of  the  neat 
molten  salt  eutectic  and  after  addition  of  nanoparticles.  The  results 
showed  that  the  specific  heat  of  the  nanofluid  was  enhanced  by 
19-24%.  The  microscopy  results  revealed  that  the  nanofluids  were 
fairly  stable  with  the  minimal  agglomeration  of  the  nanoparticles 
after  repeated  thermocycling. 

4.3.  Viscosity 

Lee  et  al.  [  129]  experimentally  investigated  the  effective  viscos¬ 
ities  and  thermal  conductivities  of  water-based  nanofluids  contain¬ 
ing  very  low  concentrations  of  A1203  nanoparticles.  The  results 
showed  that  the  alumina  nanoparticles  can  be  best  dispersed  and 
stabilized  in  DI  water  with  little  evidence  of  aggregation  at  5  or 
more  hours  of  ultrasonic  vibration.  The  results  of  viscosity  measure¬ 
ments  showed  that  the  viscosity  of  the  Al203-water  nanofluids 
decreases  significantly  with  increasing  temperature.  The  results  also 
showed  that  in  contrast  to  viscosity,  the  measured  thermal  con¬ 
ductivities  of  the  dilute  Al203-water  nanofluids  increase  closely 
linearly  with  the  concentration. 

Kole  and  Dey  et  al.  [130]  presented  investigation  of  thermal 
conductivity  and  viscosity  of  the  nanofluids  as  a  function  of  volume 
fraction  of  A1203  nanoparticles  as  well  as  temperature  between 
10  °C  and  80  °C.  Many  of  the  suspensions  containing  A1203  nano¬ 
particles  (  <  50  nm)  in  a  car  engine  coolant  prepared  by  using  oleic 
acid  as  the  surfactant  and  are  tested  to  be  stable  for  more  than  80 
days.  The  results  indicated  that  the  volume  fraction  of  0.035  of 
A1203  nanoparticles  in  the  engine  coolant  enhanced  the  thermal 
conductivity  of  the  fluid.  This  enhancement  in  thermal  conductivity 
of  the  nanofluid  varied  linearly  with  the  volume  fraction  of  the 
nanoparticles  and  reaches  a  highest  of  11.25%  at  80  °C  for  the 
nanofluid  of  0.035  vol%  of  A1203  nanoparticles.  The  results  revealed 
that  the  viscosity  of  the  nanofluids  increases  with  alumina  volume 
fraction  and  decreases  with  the  rise  in  temperature. 

Masoumi  et  al.  [131]  introduced  a  new  equation  for  calculating 
the  nanofluid  viscosity  by  considering  the  Brownian  motion  of 
nanoparticles.  The  proposed  equation  calculates  the  nanofluid 
viscosity  as  a  function  of  the  temperature,  the  nanoparticle  volume 
fraction,  the  mean  nanoparticle  diameter,  the  nanoparticle  density 
and  the  base  fluid  physical  properties.  The  results  showed  with  the 
comparison  with  many  experimental  works  that  the  presented 
model  could  well  predict  the  effective  viscosity  of  different 
nanofluids  (CuO-H20,  A1203-H20,  CuO-EG,  Ti02-EG  and  CuO- 
EG-H20)  at  different  nanoparticle  mean  diameters,  volume  frac¬ 
tions  and  also  temperature  variations. 

Abu-Nada  and  Chamkha  [132]  presented  study  of  natural 
convection  heat  transfer  characteristics  in  a  differentially  heated 
enclosure  filled  with  a  CuO-EG-Water  nanofluid  for  different 
published  variable  thermal  conductivity  and  variable  viscosity 
models.  The  results  showed  that  for  high  Rayleigh  numbers,  the 
average  Nusselt  number  decreased  as  the  nanoparticles  volume 
fraction  increased. 
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Table  3 

The  thermal  conductivity  and  viscosity  correlations  for  nanofluids. 


Researchers 


Equation 


Remarks 


Koo  and  Kleinstreuer  [140] 


Chon  [178] 

Vajjha  and  Das  [179] 


Corcione  [180] 

Seyf-Feizbakhshi  [181] 
Nguyen  et  al.  [122] 

Masoumi  et  al.  [131] 


Abu-Nada  [173] 


Tseng  and  Lin  [175] 
Chen  et  al.  [176] 


Chen  et  al.  [177] 

Koo  and  Klennstreuer  [182] 


Vajjha  [183] 

Mahbubul  [184] 

Tseng  and  Chen  [185] 

Cheng  and  Law  [186] 

Kulkarni  et  al.  [187] 
Namburu  et  al.  [188] 

Masoud  Hosseini  et  al.  [189] 


kp  +  2  kBF + 2  (kp-kBF)<P. 
kp  +  2kBF-(kp-kBF)& 
0.013 


kBF +5  X  104^p  Cp 


(100  0)"u^ 
0.0011 


for  0  <  1% 


for  <!>  >  1% 


P 

(1000)- 
f(T,  ad)  =  (— 6. 04ad  +  0.4705)7 

+  (1722.3^—134.63) 


g=l +64.70°™°  g 


0.3690  />,  x  0.7476 

'y_  j  prO. 9955  pe^.232^ 


Kpff  = 


kp  +  2  kBF +2  (kp-  kBF)0. 
kp  +  2  kBF-(kp-kBF)0 


kBF +  5  x  W4/)0Pp  Cp  y^f(T,0) 


/(T,ad)  =  (2.8217  x  1(T20+3.917  x  10-3)+±^+ 
+  (-3.0669  x  1CT20— 3.91123  x  lO-3) 


£  =  8.4407(1000)-’  07304 

^=  1  +4.4 Re04Pr066(+)10  &003  0° 

Kf  If r  Kf 


lfr 


1 


+  =„/fcI./dv 

Pnf-Fbf  X  0  904  e01482* 

6n/-  M>/(1  +  O.O250  +  O.O1502) 
,  l>PVBd 2 

Pnf-Pbf+^j. 2CT 


1  j18KbT 


Kb 


Vb  d„\l  xpPdP" 

C  =  jtfj  [(Cfdp+C2)i/>+(C3dp+C4)] 

Pai2o3  =  0.155-'9^82  +  0.7940 +2°^'47— 0.1 9202 

8  H0  2746^863  127  3  i  6044^  +2.17544 

T  j3  T  J2 


Pnf.Pbf  X  13.47  e35-98* 

i  \  \  —2.5</)m 

.'Pm 


Pnf  -  Pbf 


Pnf  =  pbfO  +  10.60 +  (1O.60)2) 

Peff  =  Astatic  +  ^Brownian  /^static  =  Pf/(  \-$)2S 

P-Brownian  =  5x10  p(j>pf  i  / ~  5  f  (T,  0) 

y  ^Pnp^np 


ft  =  0.01 3(1  OO0)_U 
/?  =  0.001 1(1  OO0)- 


for  0<1% 
for  0  >  1% 


f  (T,  «d)  =  (-6.04ad  +  0.4705)7 + (1 722.3«d - 1 34.63) 


P-nf 

Pbf 


=  At  e(/l2 


6n/_  Pbf  X  exp  [m+a  +£(  e>b)+r  (y^T- 
pbf  X  0.4513  e0-6965* 

Pnf=  PbfO  +2.50e  +  (2.50e)2+  (2.50e)3+  (2.50e)4+  ••• ) 

<Pe  =  <p(  1+^ 


Log(/in/)=Ae-ST 

4  =  1.8375(0)2— 29.643(01  +  165.56  with  R2  =  0.9873 
6  =  4  x  lO-6(0)2— 0.001(01+0.0186  with  R2  =  0.988 

Pnf  =  Pbf  x  exp [m  +  a(+j  +£(01,1  +  1' 


1  %  <  ad  <  4%  300  <  T  <  325  K 


2pf  I<bT 

Re  =  (-++-) 

npjdp 

Where  V  is  the  volume  of  computational  domain 


Q  C2  C3  and  C4  are  constants 


A!  and  A2  depends  on  dp  and  0 


Here  A  and  B  are  functions  of  volume  percentage  0 
A,  B  are  functions  of  particle  volume  percentage  0 


Nguyen  et  al.  [133]  presented  experimental  investigation  of  the 
effect  due  to  temperature  and  particle  volume  concentration  on 
the  dynamic  viscosity  for  the  water-Al203  nanofluid.  The  viscosity 


data  were  collected  using  a  ‘piston-type'  commercial  viscometer 
for  temperatures  ranging  from  room  condition  up  to  75  °C.  Two 
different  particle  sizes  were  considered  36  nm  and  47  nm  and 
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particle  volume  fraction  varying  from  1  to  nearly  13%.  The  results 
illustrated  that  the  nanofluid  dynamic  viscosity  increases  signifi¬ 
cantly  with  particle  volume  fraction  but  it  obviously  decreases 
with  a  temperature  increase.  The  results  also  showed  that  the 
Einstein's  formula  and  some  other  ones  that  were  originated 
from  the  classical  theory  of  linear  fluid  appeared  to  be  limited  to 
nanofluids  with  low  particle  volume  fractions. 

Jia  et  al.  [134]  presented  investigation  of  the  viscosity  of  silicon 
dioxide  nanofluid  at  different  particle  sizes  and  pH  values  con¬ 
sidering  nanoparticle  aggregation.  The  experimental  and  simula¬ 
tion  results  revealed  that  the  nanoparticle  diameter  is  of  crucial 
importance  for  the  viscosity  of  a  nanofluid.  The  results  showed 
that  as  the  particle  size  is  smaller,  the  viscosity  is  much  larger  and 
is  much  more  dependent  on  the  nanofluid  volume  fraction.  The 
researcher's  correlations  for  thermal  conductivity  and  viscosity  of 
nanofluids  are  summarized  in  Table  3. 


5.  Numerical  studies  on  nanofluids 

Feng  et  al.  [135]  proposed  a  new  model  for  the  effective 
thermal  conductivities  by  taking  into  account  the  nanolayer  and 
nanoparticles  aggregation.  This  model  is  expressed  as  a  function  of 
the  thickness  of  the  nanoparticle  size,  the  nanolayer,  the  nano¬ 
particle  volume  fraction,  the  thermal  conductivities  of  suspended 
nanoparticles  and  base  fluid.  The  model  predictions  of  effective 
thermal  conductivity  showed  a  good  agreement  with  the  available 
experimental  data. 

Tillman  and  Hill  [136]  presented  a  survey  of  some  of  the  attempts 
to  model  the  enhanced  thermal  conductivity  of  such  nanofluids,  and 
address  issues  such  as  the  nanoparticles  themselves,  cluster  structure, 
the  surrounding  layer,  the  fluid  environment,  and  different  heat 
transport  processes  at  the  micro  and  nano-scales.  The  reported 
models  showed  that  the  effective  thermal  conductivity  of  nanofluids 
is  a  function  of  the  thermal  conductivities  of  base  fluid  and  nano- 
particles,  nanoparticle  size  and  shape,  aggregation,  particle  clustering 
and  fluid  temperature.  The  essential  limitations  of  conventional 
heat  conduction  models  for  solid/liquid  suspensions,  which  were 
rooted  in  macroscopic  transport  laws,  such  as  the  Fourier  law  of  heat 
conduction  diffusion  have  been  recognized.  However,  the  ballistic/ 
diffusive  conduction  has  only  briefly  been  introduced  in  solid/liquid 
system.  The  core-shell-medium,  multi-component  models  reflect  the 
nanofluid  structure  and  incorporate  the  correct  thermal  conductivity 
of  nanoparticles  which  presented  to  be  less  than  that  of  bulk  material. 

Wang  et  al.  [137]  presented  a  method  for  modeling  the 
effective  thermal  conductivity  of  nanofluid  based  on  the  effective 
medium  approximation  and  the  fractal  theory  for  the  description 
of  nanoparticle  cluster  and  its  radial  distribution.  The  surface 
adsorption  and  the  size  effect  of  nanoparticles  were  taken  into 
considerations.  The  results  revealed  to  that  the  proposed  fractal 
model  predicts  well  the  trend  for  variation  of  the  effective  thermal 
conductivity  with  dilute  suspension  of  nanoparticles,  and  fits 
successfully  with  the  experimental  data  for  50  nm  CuO  particles 
suspension  in  deionized  water  when  the  volume  concentration 
<  0.5%.  However,  the  calculated  result  also  showed  that  the 
predictive  calculation  of  effective  thermal  conductivity  is  compli¬ 
cated.  The  effective  thermal  conductivity  of  nanoparticle  suspen¬ 
sion  can  be  expressed  as: 

keff  _  (1—0)  +  30  f™(kci(r)n(r)/kci(r)  +  2kf)dr 
kf  ~  (l-0)  +  30/o”(kdn(r)/kd(r)+2/</)dr 

where  kci  is  the  effective  thermal  conductivity  of  clusters,  n(r)  is 
the  radius  distribution  function. 

Xue  and  Xu  [138]  derived  an  expression  for  the  effective 
thermal  conductivity  of  nanofluids  with  interfacial  shells.  The 
expression  is  not  only  depended  on  the  thermal  conductivity  of 


the  liquid  and  solid  and  their  relative  volume  fraction,  but  also 
depended  on  the  interfacial  properties  and  particle  size.  The 
theoretical  results  of  the  proposed  model  on  the  effective  thermal 
conductivity  of  CuO/water  and  CuO/ethylene  glycol  nanofluids 
were  in  good  agreement  with  the  experimental  data. 

9("l  v\  ke~km  v  (ke-k2)(2k2  +  k1  )-a(ki -k2)(2k2  +  ke) 

V  a)  2ke+km  a(2ke+k2)(2k2+kj)+2a(ki-k2)(k2-ke) 

(2) 

where  ke  the  effective  thermal  conductivity:  via  are  the  volume 
fraction  of  the  complex  nanoparticles;  km  and  k2  are  the  thermal 
conductivity  of  the  base  fluid  and  particle,  respectively. 

Leong  et  al.  [139]  proposed  a  model  for  predicting  the  effective 
thermal  conductivity  of  nanofluids  by  considering  the  effect  of  the 
interfacial  layer  at  the  solid  particle/liquid  interface.  It  has  been 
documented  that  the  interfacial  layer  at  the  solid  (particle)/liquid 
interface  and  particle  size  is  one  of  the  major  mechanisms  for 
enhancing  the  thermal  conductivity  of  nanofluids.  The  proposed 
model  takes  into  account  some  additional  effects  including  thickness, 
volume  fraction,  particle  size  and  thermal  conductivity  of  the  inter¬ 
facial  layer.  The  results  of  this  model  showed  good  predictions  for  the 
effective  thermal  conductivity  of  different  types  of  nanofluids.  The 
proposed  model  is  found  to  be  better  than  the  existing  models  due  to 
the  predicted  effective  thermal  conductivity  of  different  types  of 
nanofluids  were  closer  to  the  experimental  results. 

y  _(kp— kir)0ikir[2/?i—  /?3  +  l]+(kp  +  2kir)/^[0i/?3(lqr— kf)+kf] 

/)}(k„  +  2k]r)  (kp  ki,-)0i  \Jl\  -F/?3  —  1  ] 

(3) 

where  klr  and  kf  is  interfacial  layer  at  the  particle/liquid  interface; 
fh=P=\- 

Koo  and  Kleinstreuer  [140]  developed  a  new  thermal  conductiv¬ 
ity  model  for  nanofluids,  which  takes  the  effects  of  particle  volume 
fraction  and  temperature  dependence,  particle  size,  as  well  as 
properties  of  base  liquid  and  particle  phase  into  consideration  by 
considering  surrounding  liquid  traveling  with  randomly  moving 
nanoparticles.  The  nanoparticles  in  quiescent  suspension  move 
randomly  and  thereby  cany  relatively  large  volumes  of  surrounding 
liquid  with  them.  This  micro-scale  interaction  may  occur  between 
cold  and  hot  regions,  resulting  in  a  lower  local  temperature  gradient 
for  a  given  heat  flux  compared  with  the  pure  liquid  case.  Therefore, 
the  effective  thermal  conductivity,  kef{,  as  a  result  of  Brownian 
motion,  increases  to  result  in  a  lower  temperature  gradient  for  a 
given  heat  flux. 


/(7>d)  =  (—6.04«d  +  0.4705)T+(1722.3«d— 134.63)  (5) 


Yu  and  Choi  [141]  developed  a  new  model  for  the  prediction  of 
the  thermal  conductivities  of  three-phase  suspensions  of  ellipsoi¬ 
dal  particle-in-liquid  suspensions  by  extending  the  Hamilton- 
Crosser  model  to  include  the  particle-liquid  interfacial  layer.  The 
solid/liquid  interface  is  described  as  a  confocal  ellipsoid  with  a 
solid  particle.  The  model  was  utilized  to  predict  the  thermal 
conductivities  of  nanotube-in-oil  suspensions,  in  which  both  the 
nanotube-to-liquid  thermal  conductivity  ratio  and  nanotube 
eccentricity  are  extreme  values.  The  results  of  developed  model 
were  found  to  be  in  good  agreement  with  the  experimental  data. 
However,  this  new  model  is  not  able  to  predict  the  nonlinear 
behavior  of  the  nanofluid  thermal  conductivity. 

ke—k]  fg  „  kpj—k\ 

ke+(n- l)ki  -  3J=~hckPJ+(n-l)k, 
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where  ke  effective  thermal  conductivity;  k i  is  thermal  conductivity 
of  liquid  medium;  n  is  the  empirical  shape  factor;  fe,  equivalent 
volume  concentration. 

Moghadassi  et  al.  [142]  presented  a  novel  model  for  the  predic¬ 
tion  of  the  effective  thermal  conductivity  of  nanofluids  based  on 
dimensionless  groups.  The  model  expresses  the  thermal  conduc¬ 
tivity  of  a  nanofluid  as  a  function  of  the  thermal  conductivity  of 
the  solid  and  liquid,  their  particle  size,  volume  fractions  and 
interfacial  shell  properties.  The  results  of  this  model  revealed  that 
the  thermal  conductivity  changes  nonlinearly  with  nanoparticle 
loading  as  shown  in  this  equation; 


Kff 

km 


=  1  +171 


(Vpf 

(dpf 


(7) 


where  m  is  a  factor  that  depends  on  the  properties  of  the  solid 
particle  and  interfacial  shell,  while  a  and  p  are  empirical  constants. 

Nan  et  al.  [143]  proposed  a  quite  simple  formula  of  thermal 
conductivity  enhancement  in  carbon  nanotube  composites  based 
on  a  conventional  model.  The  theoretical  results  of  this  simple 
formula  predicted  much  higher  thermal  conductivity  enhance¬ 
ment  even  in  the  dilute  case  of  the  carbon  nanotubes,  due  to 
ultrahigh  thermal  conductivity  and  aspect  ratio  of  the  carbon 
nanotubes.  The  large  inconsistency  between  the  predictions  and 
current  experiments  could  be  attributed  to  interfacial  thermal 
resistance  between  the  matrix  and  nanotubes,  and  aggregation 
and  twist  of  the  nanotubes  in  the  composites.  The  results  also 
concluded  that  the  reasonable  conventional  models  were  still  valid 
for  the  carbon  nanotube  based  composites. 

Petal  et  al.  [144]  presented  a  simple  theoretical  model  for  predic¬ 
ting  the  thermal  conductivity  of  carbon  nanotube  (CNT)  nanofluids. 
The  effects  due  to  the  high  thermal  conductivity  of  CNTs  and  the 
percolation  of  heat  through  it  were  considered  to  be  the  most 
important  reasons  for  their  anomalously  high  thermal  conductivity 
enhancement.  The  new  approach  was  taken  for  the  modeling,  the 
novelty  of  which  lies  in  the  prediction  of  the  thermal  behavior  of 
water  based  as  well  as  oil  based  CNT  nanofluids,  which  are  quite 
different  from  each  other  in  thermal  characteristics.  The  results  of 
this  model  successfully  match  the  experimental  data,  varying  over 
two  orders  of  magnitude  in  their  thermal  conductivity  enhance¬ 
ment  values.  The  effective  thermal  conductivity  equation  is  shown 
as  follow: 


fceff  =  k] 


1  + 


kser  i 

/<i(l-e)rs 


(8) 


kj  is  thermal  conductivity  of  liquid  medium;  ks  thermal  conductivity 
of  particle;  e,  particle  volume  fraction;  r,  is  the  particle  radius. 

Xue  [145]  presented  a  novel  model  of  the  effective  thermal 
conductivity  for  nanofluids  by  considering  the  interface  effect 
between  the  solid  particles  and  the  base  fluid  in  nanofluids.  The 
formula  of  calculating  the  effective  thermal  conductivity  of  nano¬ 
fluids  given  based  on  Maxwell  theory  and  average  polarization. 
The  theoretical  results  revealed  that  the  effective  thermal  con¬ 
ductivity  of  nanotube/oil  nanofluid  and  Al203/water  nanofluid 
were  in  good  agreement  with  the  experimental  data. 


\  ke—km  V 

ke-kcji  4 

ke  key 

)2ke+km  X 

Ke +^2(^c,x  —  ke) 

2/Ce  +  (1  — B2  ,x)(kc,y  —  ke)_ 

0) 


where  kcx  and  kc%y  are  the  thermal  conductivity  components  of  the 
complex  elliptical  particle  along  the  x  and  y  axes,  respectively. 
v  and  u/X  are  the  volume  fractions  of  nanoparticles  and  complex 
nanoparticles,  respectively. 

Xu  and  Yu  [146]  reported  a  new  model  for  predicting  the 
thermal  conductivity  of  nanofluids  by  taking  into  account  the 
fractal  distribution  of  nanoparticle  sizes  and  heat  convection 
between  nanoparticles  and  liquids  due  to  the  Brownian  motion 


of  nanoparticles  in  fluids.  The  proposed  model  is  expressed  as  a 
function  of  the  average  size  of  nanoparticles,  concentration  of 
nanoparticles,  fractal  dimension,  temperature  and  properties  of 
fluids.  The  results  of  this  model  showed  that  the  reasonable 
dependences  of  the  thermal  conductivity  on  the  temperature  of 
nanofluids,  nanoparticle  size  and  concentration.  The  predictions  of 
this  model  were  in  good  agreement  with  the  available  experi¬ 
mental  data.  The  total  effective  thermal  conductivity  of  nanofluids 
is  written  as: 

kj>f  I +  i +  kp  +  2kf—20(kf—kp)  Nu  /  dj(2  Pf  — 1]21 

kf  ~  5  c  ~ kp+2kf+2e>(,kf-kp)+  pr  (l-D;)2  (W4nta)2“D' -1 

(10) 


Df  =  d—. - 
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Mirmasoumi  and  Behzadmehr  [147]  presented  numerical  study 
of  fully  developed  mixed  convection  of  a  nanofluid  ( water/ A1203). 
Two-phase  mixture  model  was  utilized  to  analyze  the  effects  of 
nanoparticles  mean  diameter  on  the  flow  parameters.  The  results 
showed  that  the  convection  heat  transfer  coefficient  significantly 
increases  with  decreasing  the  nanoparticles  means  diameter. 
However  it  does  not  significantly  change  the  skin  friction  coeffi¬ 
cient.  The  thermal  conductivity  correlation  that  used  in  this  study 
as  follow  [178], 

Behzadmehr  et  al.  [148]  presented  numerical  study  of  turbulent 
forced  convection  heat  transfer  in  a  circular  tube  with  a  nanofluid 
consisting  of  water  and  1  vol%  Cu.  Two  phase  mixture  model  has 
been  implemented  to  study  such  a  flow  field.  The  results  showed 
that  the  mixture  model  is  more  accurate  than  the  single  phase 
model.  However,  the  accuracy  of  mixture  model  it  seems  could 
improve  by  using  suitable  effective  physical  properties  for  nano¬ 
fluid  instead  of  volume  weighted  average  of  particle  and  fluid 
properties.  The  results  demonstrated  that,  adding  1%  nanoparticles 
(Cu)  increases  the  Nusselt  number  more  than  15%  while  it  does  not 
have  any  significant  effect  on  the  skin  friction. 

Ben-Abdullah  [149]  presented  numerical  simulation  using  the 
Landauer-Buttiker  theory  to  calculate  the  thermal  conductance 
associated  to  plasmons  modes  in  one  dimensional  array  of  nano¬ 
particles  closely  spaced  in  a  host  fluid.  The  results  showed  that  the 
contribution  of  near-field  interactions  to  the  transport  of  heat  in 
weakly  concentrated  nanofluids  when  the  average  separation  dis¬ 
tance  between  particles  is  small  in  front  of  plasmons  wavelengths. 


6.  Theoretical  studies  on  nanofluids 

Perez-Madrid  et  al.  [150]  analyzed  heat  transfer  between  two 
nanoparticles  separated  by  a  distance  lying  in  the  near-field  domain 
in  which  energy  interchange  is  due  to  the  Coulomb  interactions.  The 
thermal  conductance  was  computed  by  assuming  that  the  particles 
have  charge  distributions  characterized  by  fluctuating  multipole 
moments  in  equilibrium  with  heat  baths  at  two  different  tempera¬ 
tures.  This  quantity  follows  from  the  fluctuation-dissipation 
theorem  for  the  fluctuations  of  the  multipolar  moments.  The  results 
revealed  that  the  analysis  of  the  heat  interchanged  between  two 
NPs  separated  by  a  few  submicrons  agrees  with  the  explainations  of 
the  rapid  growth  of  the  conductance  even  when  the  NPs  are  in 
contact. 

Putnam  et  al.  [151]  described  an  optical  beam  deflection  tech¬ 
nique  for  measurements  of  the  thermal  diffusivity  of  fluid  mixtures 
and  suspensions  of  nanoparticles  with  a  precision  of  better  than  1%. 
The  proposed  approach  was  tested  by  using  the  thermal  conductiv¬ 
ity  of  ethanol-water  mixtures.  The  results  revealed  that  there  was 
no  significant  enhancement  in  the  thermal  conductivity  of  fluids 
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that  are  loaded  by  small  volume  fractions  of  nanoparticles  and  the 
largest  increase  in  thermal  conductivity  observed  is  1.3%  +  0.8%  for 
4  nm  diameter  Au  particles  suspended  in  ethanol. 

Eapen  et  al.  [152]  presented  an  investigation  of  the  conduction  in 
nanofluids.  The  classical  Maxwell  theory  has  two  limiting  bounds 
that  correspond  to  two  geometrical  nanoparticle  configurations 
in  the  colloidal  state.  For  the  first  bound,  nanoparticles  form  the 
dispersed  phase  and  the  fluid  medium  acts  as  the  continuous  phase. 
Whereas,  for  the  second  the  roles  are  reversed  with  the  nanopar¬ 
ticles  and  the  fluid  medium  constituting  the  continuous  and 
dispersed  phases,  respectively.  The  results  showed  that  almost  all 
the  reported  thermal  conductivity  data  are  enveloped  by  the  upper 
and  lower  Maxwell  bounds.  Moreover,  for  a  homogeneous  nano¬ 
fluid  system,  these  bounds  organize  the  limits  for  thermal  con¬ 
ductivity  enhancements.  The  thermal  conductivity  for  the  upper 
bound  is  maximally  biased  toward  the  thermal  conductivity  of  the 
nanoparticles  and  vice  versa  for  the  lower  bound. 

Moghaddami  et  al.  [153]  presented  analytically  the  effects  of 
adding  nanoparticles  on  the  entropy  generation  of  water-Al203 
and  ethylene  glycol  -A1203  nanofluid  flows  through  a  circular  pipe 
under  uniform  wall  heat  flux  thermal  boundary  condition  in  both 
laminar  and  turbulent  regimes.  The  results  obtained  for  water- 
A1203  nanofluid  illustrated  that  adding  1%  to  the  nanoparticle 
volume  concentration  for  laminar  flows  decreases  the  entropy 
generation  about  3.6%.  The  entropy  generation  at  the  Re= 853 
minimized  for  each  specified  nanoparticles  volume  concentration. 
The  results  showed  also,  at  the  Re=  853  and  0=5%,  maximum 
reduction  in  the  entropy  generation  occurs. 


7.  Single  phase  and  two  phase  models 

Liu  et  al.  [154]  presented  an  experimental  investigation  of 
single-phase  forced  convection  of  Al203-water  nanofluid  in  a 
circular  minichannel  with  a  1.09  mm  inner  diameter.  The  Reynolds 
number  was  varied  from  600  to  4500,  covering  the  laminar, 
transition,  and  early  fully  developed  turbulent  regions.  The  results 
revealed  that  in  the  laminar  region,  the  nanofluids  exhibit  pro¬ 
nounced  entrance  region  behaviors  possibly  due  to  the  flattening 
of  the  velocity  profile  caused  by  the  flow-induced  particle  migra¬ 
tion.  Furthermore,  the  convective  heat  transfer  of  nanofluids  is 
enhanced  in  the  laminar  flow  with  the  penalty  of  increased 
pressure  drop.  The  increasing  in  both  pressure  drop  and  heat 
transfer  were  proportional  to  the  nanoparticle  volume  concentra¬ 
tion.  The  results  revealed  that  the  nanofluids  should  be  used  in 
either  the  laminar  flow  or  the  fully  developed  turbulent  flow  with 
sufficiently  high  Reynolds  number  in  order  to  yield  enhanced  heat 
transfer  performance  for  engineering  applications. 

Kumar  et  al.  [155]  presented  numerical  simulations  for  analysis  of 
flow  and  thermal  field  in  a  thermally  driven  cavity  using  nanofluid. 
The  controlling  parameters  considered  for  this  investigation  were 
Grashof  number  (103  <  Gr<  105),  solid  volume  fraction  (0  <  <f>  <  0.2) 
and  empirical  shape  factor  (0.5  <  n  <  6).  The  results  showed  that  the 
dispersed  thermal  conductivity  is  more  in  the  region  near  the  walls 
where  the  magnitude  of  velocity  is  more  compared  to  the  center  of 
the  domain.  The  results  also  revealed  that  maximum  value  of  the 
dispersed  thermal  conductivity  increases  when  the  hydraulic  dia¬ 
meter  of  the  particles  is  increased. 


8.  Mechanism  of  heat  conduction 

Keblinski  et  al.  [156]  presented  a  critical  analysis  of  the  experi¬ 
mental  data  in  terms  of  the  potential  mechanisms.  The  considera¬ 
tions  of  this  investigation  strongly  suggest  that  the  apparent 
disagreement  between  the  experiment  and  the  effective  medium 


theory  is  simply  an  artifact  of  focusing  on  a  particular  effective 
medium  theory,  namely,  Maxwell's  theory  of  well-dispersed  parti¬ 
cles.  The  results  illustrated  that,  by  accounting  for  linear  particle 
aggregation,  the  well  established  effective  medium  theories  for 
composite  materials  are  capable  of  explaining  the  vast  majority  of 
the  reported  data  without  resorting  to  novel  mechanisms  such  as 
Brownian  motion  induced  nanoconvection,  liquid  layering  at  the 
interface,  or  near-field  radiation.  The  results  also  showed  that  the 
effective  medium  theories  are  competent  of  predicting  thermal 
conductivity  enhancements  well  beyond  the  Maxwell  prediction. 
The  results  revealed  to  that  the  significant  aggregation  is  required  to 
obtain  substantial  increases  in  thermal  transport  has  an  important 
consequence  on  the  potential  application  of  such  fluids  in  flow 
based  cooling. 

Keblinski  [157]  presented  a  study  of  potential  mechanisms  that 
were  put  forward  in  order  to  understand  nanofluid  thermal  con¬ 
ductivity  and  demonstrate  that  particle  aggregation  is  currently  the 
only  physically  reasonable  mechanism  that  can  explain  the  majority 
of  the  experimental  data.  The  results  revealed  that  significant 
aggregation  is  required  to  obtain  substantial  increases  in  thermal 
transport  has  an  important  consequence  for  the  potential  applica¬ 
tion  of  such  fluids  in  flow-based  cooling,  which  is  the  most 
important  benefit  from  the  technological  point  of  view. 

8.3.  Brownian  motion 

Brownian  motion  is  one  of  the  parameters  that  affect  the  flow 
and  heat  transfer  of  nanofluids.  It  is  basically  the  random  dynamic 
mode  of  particles  in  a  liquid  where  the  particles  collide  between 
each  other.  It  was  determined  that  the  thermal  conductivity 
increases  with  an  increase  in  temperature,  but  also  it  was  shown 
that  nanofluids  compose  of  smaller  particles  experienced  greater 
enhancement  than  with  larger  particles.  Since  temperature  repre¬ 
sents  the  overall  kinetic  energy  of  the  particles,  an  increase  in 
temperature  will  cause  increase  in  the  particles  motion.  It  is  easier 
for  smaller  particles  to  move;  therefore,  smaller  particles  will 
display  a  higher  level  of  Brownian  motion  than  larger  particles. 
Some  researcher  conducted  studies  in  this  motion,  Xuan  et  al. 
[158]  applied  the  theory  of  Brownian  motion  and  the  diffusion- 
limited  aggregation  model  to  simulate  random  motion  and  the 
aggregation  process  of  the  nanoparticles.  By  considering  physical 
properties  of  both  nanoparticles  and  the  base  liquid,  as  well  as  the 
structure  of  the  nanoparticles,  aggregates,  a  theoretical  model  has 
been  developed  to  predict  the  thermal  conductivity  of  nanofluids. 
The  results  obtained  from  this  model  were  in  satisfactory  agree¬ 
ment  with  the  experimental  results  especially  when  the  effect  of 
nanoparticle  aggregation  is  taken  into  account.  A  theoretical 
model  is  developed  to  predict  the  thermal  conductivity  of  nano¬ 
fluids.  Comparison  between  the  theoretical  and  experimental 
results  shows  the  validity  and  accuracy  of  the  theoretical  model. 
The  results  revealed  to  that,  the  thermal  conductivity  of  a  nano- 
fluid  increases  when  the  fluid  temperature  increase. 

Evans  et  al.  [159]  presented  investigation  of  using  simple 
kinetic  theory  based  analysis  of  heat  flow  in  fluid  suspensions  of 
solid  nanoparticles  (nanofluids)  to  demonstrate  that  the  hydro¬ 
dynamics  effects  associated  with  Brownian  motion  have  a  minor 
effect  on  the  thermal  conductivity  of  the  nanofluid.  The  results  of 
kinetic  theoiy  argument  and  the  results  of  molecular  dynamics 
simulations,  revealed  to  that  the  thermal  conductivity  of  a  nano- 
fluid  with  well  dispersed  nanoparticles  is  well  described  by  the 
effective  medium  theory  and  did  not  show  any  significant 
enhancements  due  to  effects  associated  with  Brownian  motion 
of  nanoparticles. 

Nie  et  al.  [160]  considered  several  mechanisms  that  might 
account  for  the  high  enhancement  of  thermal  conductivity  observed 
by  many  groups  and  presented  the  exact  result  of  enhancement  of 
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thermal  conductivity  due  to  the  Brownian  motion  and  estimated  it 
is  order  of  magnitude.  The  results  revealed  that  the  possibility  of 
convection  caused  by  Brownian  particles  was  found  to  be  unlikely. 
Analyzing  the  possibility  of  the  micro  convection  in  the  nanofluid 
concluded  that  this  phenomenon  cannot  happen  because  of  the 
quick  decay  of  the  hydrodynamic  field  around  the  Brownian 
nanoparticles.  The  phonon  mean  free  path  and  phonon  transport 
speed  in  nanofluid  due  to  the  existence  of  suspending  nanoparticles 
also  discussed.  The  authors  found  that  only  a  thin  layer  of  one 
nanometer  may  form  around  the  nanoparticle  and  most  of  the  base 
fluid  molecules  will  remain  unperturbed;  hence  there  is  no 
improvement  on  the  structure  of  the  fluid  phase.  The  exploring 
the  compressibility  of  the  nanofluid  with  single  spherical  nano¬ 
particle  or  aggregates  of  nanoparticles  suspended,  revealed  to 
negligible  changes  in  compressibility,  which  indicates  the  speed  of 
phonon  transport  found  to  be  negligible  due  to  the  existence  of 
nanoparticles  in  the  low  volume  fraction  limit. 

8.2.  Nature  of  heat  transport  in  nanoparticles 

Heat  is  transferred  via  conduction  implies  that  heat  is  transferred 
internally,  by  vibrations  of  atoms  and  molecules.  The  vibrations  of  the 
atoms  which  are  connected  together  as  give  rise  to  the  vibrations  of 
the  whole  crystal,  that  is,  to  lattice  vibrations.  The  energy  of  the  whole 
vibrating  system  is  quantized,  and  the  quantum  of  thermal  energy 
absorbed  or  emitted  by  an  atom  is  called  a  phonon.  Therefore,  a 
phonon  is  essentially  a  quantized  mode  of  vibration  occurring  in  rigid 
ciystal  lattice  and  plays  a  major  role  in  material's  thermal  conduction. 
Buongiorno  [161]  presented  development  of  two-component  four- 
equation  nonhomogeneous  equilibrium  model  for  mass,  momentum, 
and  heat  transport  in  nanofluids.  Nondimensional  analyses  of  the 
equations  proposed  that  the  energy  transfer  by  nanoparticle  disper¬ 
sion  is  negligible,  and  thus  cannot  explain  the  abnormal  heat  transfer 
coefficient  increases.  Furthermore,  a  comparison  of  the  nano¬ 
particle  and  turbulent  eddy  time  and  length  scales  clearly  indicates 
that  the  nanoparticles  move  homogeneously  with  the  fluid  in  the 
presence  of  turbulent  eddies,  so  an  effect  on  turbulence  intensity  is 
doubtful.  Thus,  there  was  alternative  explanation  proposed  for  the 
abnormal  heat  transfer  coefficient  increases.  For  a  heated  fluid,  these 
effects  can  result  in  a  significant  decrease  of  viscosity  within  the 
boundaiy  layer,  thus  leading  to  heat  transfer  enhancement.  The  results 
showed  that  the  Brownian  diffusion  and  thermophoresis  have  been 
identified  as  the  two  most  important  nanoparticle/base-fluid  slip 
mechanisms. 

Keblinski  et  al.  [162]  presented  investigation  of  development  the 
fundamental  understanding  of  heat  transport  in  solid  nanoparticles 
colloids  under  stationary  conditions  and  evaluating  the  extent  to 
which  four  specific  mechanisms  could  contribute  to  the  thermal 
conductivity.  The  four  possible  explanitions  for  this  anomalous 
increase;  Brownian  motion  of  the  particles,  molecular-level  layering 
of  the  liquid  at  the  liquid/particle  interface,  the  nature  of  heat 
transport  in  nanoparticles  and  the  effects  of  nanoparticle  clustering 
were  explored.  The  results  revealed  that  in  a  particle  with  a 
diameter  smaller  than  35  nm,  the  heat  transport  is  not  diffusive, 
but  heat  is  transported  ballistically.  Although  this  fact  prevents  the 
application  of  conventional  theories  for  the  modeling  of  thermal 
conductivity  of  nanofluids,  the  results  illustrated  that  the  ballistic 
heat  transport  still  cannot  explain  the  anomalous  thermal  conduc¬ 
tivity  enhancements,  and  that's  attributed  to  the  temperature  inside 
the  nanoparticles  is  closely  constant  and  this  fact  does  not  depend 
on  whether  heat  is  transported  by  diffusion  or  ballistically. 

8.3.  Nanolayer  thickness 

There  are  some  studies  declared  that  there  is  a  liquid  layering 
on  the  nanoparticles,  which  improves  the  enhancement  of 


nanofluid  heat  transfer  properties.  Yu  and  Choi  [163]  presented  a 
modification  of  Maxwell  model,  which  includes  a  nanolayer.  The 
modified  model  showed  the  ability  of  measuring  the  increases  of 
effective  volume  fraction  and  subsequently  the  thermal  conduc¬ 
tivity  of  nanofluids,  especially  when  particle  diameter  is  <  10  nm. 
The  results  of  this  model  revealed  to,  when  the  nanolayers  are 
accounted  for,  there  is  potential  for  up  to  an  eight-fold  increase  in 
the  thermal  conductivity  of  nanofluids  over  the  enhancement 
predicted  by  the  Maxwell  model  without  the  nanolayer. 

Xie  et  al.  [164]  proposed  a  model  including  the  consideration  of 
the  role  of  nanolayer  at  particle-fluid  interface  for  predicting  the 
effective  thermal  conductivities  of  nanoparticle-fluid  mixtures. 
The  results  revealed  that  the  nanolayer  thickness,  volume  fraction, 
nanoparticle  size  and  thermal  conductivity  ratio  of  particle  to  fluid 
have  effects  on  the  enhanced  thermal  conductivity  ratios. 

Xue  et  al.  [165]  examined  the  effect  of  nanolayer  by  using  non¬ 
equilibrium  molecular  dynamics  simulations  in  which  a  tempera¬ 
ture  gradient  is  imposed.  The  simulations  results  of  a  simple 
monoatomic  liquid  revealed  that  the  layering  of  the  liquid  atoms 
at  the  liquid-solid  interface  does  not  have  any  significant  effect  on 
thermal  transport  properties. 

Tillman  and  Hill  [166]  presented  investigation  of  determine 
the  nanolayer  thickness  and  the  thermal  conductivity  profile 
within  the  nanolayer.  The  nanolayer  of  a  nanoparticle  can  be 
also  considered  as  a  monolayer  from  the  adsorption  of  liquid 
molecules  on  the  particle  surface.  The  accuracy  and  robustness  of 
the  present  method  are  verified  by  comparing  the  thermal  con¬ 
ductance  of  the  solid-fluid  interface  calculated  from  an  inverse 
estimation  of  the  nanolayer  thickness  obtained  from  observed 
experimental  values.  The  results  of  molecular  dynamic  simulation 
showed  that  the  thermal  conduction  in  the  nanolayer  can  vary 
greatly  depending  on  the  type  of  bonding  between  the  liquid  and 
the  solid. 


8.4.  Effects  of  clustering 

Clustering  of  particles  is  one  of  the  major  problems  that  affects 
of  the  enhancement  of  thermal  conductivity  by  creating  path  of 
lower  thermal  resistance.  It  is  noticed  that  as  the  particles 
agglomerate  increases,  the  size  increases  and  becomes  denser 
where  they  clustered  down,  and  end  with  low  thermal  conductiv¬ 
ity  enhancement.  Karthikeyan  et  al.  [167]  presented  investigation 
of  the  parameters  influencing  the  thermal  conductivity  enhance¬ 
ment  in  water  and  ethylene  glycol  based  nanofluids  of  CuO 
nanoparticles  of  average  diameter  8  nm.  The  results  illustrated 
that  the  thermal  conductivity  enhancement  observed  in  water  and 
ethylene  glycol  based  nanofluids  with  1  vol%  CuO  nanoparticles 
loading  are  31.6%  and  54%,  respectively.  The  finer  particle  size  and 
monodispersity  of  nanoparticles  are  the  main  reasons  attribute  to 
the  large  enhancement  in  thermal  conductivity.  The  results 
revealed  that  the  thermal  conductivity  of  the  nanofluid  increases 
nonlinearly  with  the  volume  fraction  of  nanoparticles.  The  time- 
dependent  thermal  conductivity  in  water  based  CuO  nanofluid 
showed  that  the  thermal  conductivity  decreases  with  elapsed 
time  due  to  clustering  of  nanoparticles  with  time,  as  confirmed 
microscopically. 

Evans  et  al.  [168]  proposed  investigation  of  analyzing  the  role 
of  aggregation  and  interfacial  thermal  resistance  on  the  effective 
thermal  conductivity  of  nanofluids  and  nanocomposites.  The 
results  showed  that  the  thermal  conductivity  of  nanofluids  and 
nanocomposites  can  be  significantly  enhanced  by  the  aggregation 
of  nanoparticles  into  clusters.  Moreover,  the  results  also  revealed 
to  that  the  conductivity  enhancement  due  to  aggregation  is  also  a 
strong  function  of  the  chemical  dimension  of  the  aggregates  and 
the  radius  of  gyration  of  the  aggregates. 
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9.  Heat  transfer  enhancement 

Sommers  and  Yerkes  [169]  investigated  the  thermal-hydraulic 
performance  of  dilute  suspensions  of  10  nm  aluminum  oxide 
nanoparticles  in  propanol  (0.5  wt%,  1  wt%,  and  3  wt%).  The  changes 
in  thermal  conductivity  with  particle  concentration,  density, 
specific  heat  were  measured  and  found  to  be  linear,  whereas 
changes  in  viscosity  were  nonlinear  and  increased  sharply  with 
particle  loading.  Nanofluid  heat  transfer  performance  data  were 
generally  commensurate  with  that  measured  for  the  baseline.  The 
results  also  showed  that  for  the  1  wt%  concentration,  a  small  but 
pronounced  enhancement  in  the  heat  transfer  coefficient  was 
recorded  for  1800  <Re<  2800,  due  to  an  earlier  transition  to 
turbulent  flow.  However,  for  the  1  wt%  Al203/propanol  nanofluid 
for  ReD  <  2100,  the  pressure  drop  was  found  to  increase  from  400% 
to  600%.  In  the  case  of  high  particle  loading  (i.e.  3  wt%),  the 
thermal  performance  was  observed  to  deteriorate. 

Mansour  et  al.  [170]  analyzed  the  effect  of  utilizing  different 
correlations  for  the  determination  of  the  thermophysical 
properties  of  nanofluids  on  the  results  of  the  thermal  per¬ 
formance  analyses  of  nanofluids.  For  both  laminar  and  turbulent 
flow  cases,  fully  develop  condition  was  examined  and  nanofluids 
were  considered  as  single  phase  fluids.  Different  expression  used 
for  the  determination  of  specific  heat,  dynamic  viscosity  and 
thermal  conductivity  were  presented  and  the  discrepancy 
between  them  was  demonstrated.  Two  different  sets  of  correla¬ 
tions  utilized  for  performing  the  analyses  of  pressure  drop 
and  heat  transfer  of  a  simple  configuration  for  the  determination 
of  the  mentioned  properties.  The  results  showed  that  for  a 
different  particle  volume  fractions  two  cases  provided  different 
results.  The  results  illustrated  that  for  the  laminar  flow  condition, 
the  results  were  contradictory  not  only  quantitatively  and 
qualitatively. 

Guo  et  al.  [171]  presented  investigation  of  homogeneous  and 
stable  magnetic  nanofluids  containing  y-Fe203  nanoparticles  which 
prepared  by  using  a  two-step  method,  and  their  thermal  transport 
properties.  The  proposed  technical  route  included  preparing  stable 
nanofluids  composed  of  y-Fe203  nanoparticles  and  the  mixture  of 
deionized  water  (DW)  and  ethylene  glycol  (EG)  (DW-EG)  as  the 
base  fluid.  Sodium  oleate  was  used  as  surfactant,  and  it  was  proved 
to  be  beneficial  to  the  dispersion  of  the  nanoparticles  in  the 
nanofluids.  The  results  showed  that  the  viscosity  of  the  y-Fe203 
nanofluids  fits  Newtonian  behavior  and  strongly  depends  on  the 
temperature  and  the  volume  fraction.  However,  the  viscosity 
measurements  revealed  to  that  the  nanofluids  illustrated 
Newtonian  behavior  and  the  viscosity  of  the  nanofluids  depended 
strongly  on  the  tested  temperatures  and  the  nanoparticles  loadings. 
The  results  demonstrated  also  that  the  thermal  conductivities  of  the 
nanofluids  were  measured  to  be  higher  than  that  of  base  fluid,  and 
the  enhanced  values  increase  with  the  volume  fraction  of  the 
nanoparticles. 

Mosavian  et  al.  [172]  studied  in  this  investigation  the  heat 
transfer  enhancement  of  laminar  flow  in  circular  tube  utilizing 
different  nanofluids  including  A1203  (20  nm),  CuO  (50  nm),  and  Cu 
(25  nm)  nanoparticles  in  water.  Constant  wall  temperature  was 
adopted  as  thermal  boundary  condition.  The  results  indicated  that 
the  heat  transfer  enhanced  with  increasing  nanoparticle  concen¬ 
trations  and  Peclet  number.  The  heat  transfer  coefficient  enhance¬ 
ments  for  all  nanofluids  were  higher  than  the  prediction  of  single 
phase  heat  transfer  correlation  with  nanofluids  properties  (homo¬ 
geneous  model).  The  results  also  demonstrated  that  the  nanofluids 
that  contain  metal  nanoparticles  (Cu/water)  exhibited  more 
enhancement  compared  to  oxide  nanofluids  (Al203/water  and 
CuO/water).  Furthermore,  the  CuO/water  nanofluid  showed  less 
enhancement  in  heat  transfer  coefficient  at  the  same  nanoparticle 
concentration  and  Peclet  number.  However,  that  is  attributed  to 


the  large  particle  size  and  high  viscosity  of  CuO/water  and  high 
thermal  conductivity  of  Cu/water  nanofluid. 

Abu-Nada  [173]  presented  an  investigation  of  Heat  transfer 
enhancement  in  horizontal  annuli  using  variable  properties  of 
Al203-water  nanofluid.  Different  viscosity  and  thermal  conductivity 
models  are  used  to  evaluate  heat  transfer  enhancement  in  the 
annulus.  The  base  case  uses  in  this  investigation  the  Chon  et  al. 
expression  for  conductivity  and  the  Nguyen  et  al.  experimental  data 
for  viscosity  which  take  into  account  the  dependence  of  these 
properties  on  temperature  and  nanoparticle  volume  fraction.  The 
results  showed  that  for  Rayleigh  number  Ra  >  104,  the  average  Nusselt 
number  was  reduced  by  increasing  the  volume  fraction  of  nanopar¬ 
ticles.  Accordingly,  the  average  Nusselt  number  increased  by  increas¬ 
ing  the  volume  fraction  of  nanoparticles  for  Ra=  103.  The  results 
revealed  that  for  Ra  >  104  and  ^  >  5%,  the  average  Nusselt  number  is 
reduced  by  increasing  the  volume  fraction  of  nanoparticles. 

10.  Rheological  behavior  of  nanofluids 

Nanofluids  are  expected  to  be  used  under  flow  conditions  and 
the  flow  of  suspension  which  are  sometimes  drastically  different 
from  that  of  most  common  heat  transfer  fluids  that  have  New¬ 
tonian  characteristics,  so  it  is  essential  to  have  the  rheological 
properties  of  nanofluid  to  use  practically.  Rheological  properties 
can  provide  the  knowledge  on  the  microstructure  under  both 
static  and  dynamic  conditions,  which  are  important  to  understand 
the  mechanism  of  heat  transfer  enhancement  using  nanofluids. 
Chen  et  al.  [174]  performed  an  experimental  work  and  theoretical 
analyses  on  the  rheological  behavior  of  nanofluids  containing  rod¬ 
like  nanoparticles.  Ethylene  glycol  (EG)  and  titanate  nanotubes 
(TNT)  are  used  as  the  base  liquid  and  nanoparticles,  respectively, 
to  formulate  the  nanofluids.  The  effects  of  particle  concentration 
(0.5  wt%,  1.0  wt%,  2.0  wt%,  4.0  wt%  and  8.0  wt%)  and  temperature 
(20  °C,  30  °C,  40  C,  50  C  and  60  °C)  were  examined.  The  results 
show  that  a  very  strong  shear  thinning  behavior  of  the  TNT 
nanofluids  and  big  influences  of  particle  concentration  and  tem¬ 
perature  on  the  zero  shear  viscosity  (ZSV)  and  high  shear  viscosity 
(HSV).  Therefore,  a  higher  particle  concentration  gives  a  stronger 
shear  thinning,  a  higher  zero  shear  viscosity  (ZSV)  and  a  higher 
high  shear  viscosity  (HSV).  The  results  presented  also  that  the 
temperature  has  a  strong  effect  on  the  rheological  behaviur  of  the 
EG-TNT  nanofluids.  However  the  higher  shear  viscosity  increment 
is  independent  of  the  temperature. 

Tseng  and  Lin  [175]  investigated  the  rheological  behavior  and 
suspension  structure  of  anatase  titanium  dioxide  (Ti02)  nanopar¬ 
ticles  dispersed  in  pure  water  over  a  range  of  volumetric  solids 
concentrations  (^=0.05-0.12)  and  shear  rates  (y=101-103  s~’). 
The  suspensions  of  nanoparticle  generally  exhibited  a  pseudoplas¬ 
tic  flow  behavior,  indicating  an  existence  of  particle  aggregations 
in  the  liquid  medium.  The  suspensions  became  clearly  thixo¬ 
tropic  as  (j>  was  increased  above  0.1.  The  results  also  revealed  a 
pronounced  increase  in  the  degree  of  particle  interactions  as  the 
volumetric  concentration  increased.  Fractal  dimension  (Dy)  was 
estimated  from  the  suspension  yield-stress  and  c/>  dependence,  and 
was  determined  as  Dys; /1.46-/1.78  for  the  flocculated  nanopar¬ 
ticle  suspensions. 

Chen  et  al.  [176]  presented  experimental  work  and  theoretical 
analyses  with  aims  of  more  fundamental  understanding  of  the 
rheological  behavior  of  nanofluids  and  interpretation  of  the  dis¬ 
crepancy  in  the  reported  rheological  behavior  of  nanofluids  in  the 
recent  literature.  The  experiments  were  carried  out  on  EG-based 
Ti02  nanofluids  containing  0.5%,  1.0%,  2.0%,  4.0%  and  8.0%  by 
weight  spherical  nanoparticles  at  20-60  °C.  The  theoretical  ana¬ 
lyses  were  carried  out  on  the  high  shear  viscosity,  shear  thinning 
behavior  and  temperature  dependence.  The  results  revealed  to 
that  the  EG-based  nanofluids  showed  Newtonian  behavior  under 
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the  conditions  of  this  work  with  shear  viscosity  as  a  strong 
function  of  temperature  and  particle  concentration.  Furthermore, 
the  relative  viscosity  depends  on  particle  concentration  in  a 
nonlinear  manner  but  is  independent  of  temperature.  The 
Krieger-Dougherty  equation  can  be  use  to  predict  the  high  shear 
viscosity  if  the  solids  volume  fraction  is  replaced  by  the  volume 
fraction  of  nanoparticle  aggregates.  The  results  showed  that  for 
the  spherical  nanoparticles,  an  aggregate  size  of  approximately 
3  times  the  primary  nanoparticle  size  gives  the  best  prediction  of 
the  experimental  data. 

Chen  et  al.  [177]  presented  experimental  work  on  the  rheological 
behavior  of  ethylene  glycol  based  nanofluids  containing  titanate 
nanotubes  over  20-60  °C  and  a  particle  mass  concentration  of  0-8%. 
The  results  revealed  to  that  EG-TNT  nanofluids  are  non-Newtonian 
exhibiting  shear-thinning  behavior  over  20-60  °C  and  a  particle 
mass  concentration  range  of  0-8%,  in  contrast  to  the  Newtonian 
behavior  for  EG-Ti02  nanofluids  containing  spherical  particles. 
Furthermore,  the  non-Newtonian  shear-thinning  behavior  becomes 
stronger  at  higher  temperatures  or  higher  concentrations.  The 
results  showed  for  a  given  particle  concentration,  there  exists  a 
certain  shear  rate  below  which  the  viscosity  increases  with  increas¬ 
ing  temperature,  whereas  the  opposite  occurs  above  such  a  shear 
rate.  The  results  also  presented  that  the  normalized  high-shear 
viscosity  is  independent  of  temperature. 

11.  Challenges  of  using  nanofluids  to  enhance  the  heat 
transfer  for  MT  and  MC 

The  use  of  nanofluids  seems  attractive  in  a  broad  range  of 
applications  as  reported  in  the  previous  sections.  But  the  devel¬ 
opment  in  the  area  of  nanofluid  application  is  hindered  by 
different  reasons  like  the  lack  of  agreement  of  results  obtained 
by  different  researchers;  the  poor  characterization  of  suspensions, 
and  the  lack  of  theoretical  understanding  of  the  mechanisms 
responsible  for  changes  in  properties.  This  review  concludes 
several  important  issues  that  should  receive  greater  attention  in 
the  near  future  such  as  the  thermal  conductivity,  Brownian  motion 
of  particles,  particle  migration,  and  thermophysical  property 
change  with  temperature  [190],  However,  recent  studies  indicate 
particle  aggregation  and  deposition  in  MT  and  MC.  Further  study  is 
required  in  these  areas  to  identify  the  reasons  for  and  the  effects  of 
particle  deposition.  Finally,  there  appears  to  be  hardly  any  research 
in  the  use  of  nanofluids  as  working  fluids.  Nanoparticle  disper¬ 
sions  in  two-phase  heat  transfer  applications  can  be  studied  to 
explore  the  possibility  of  improving  the  heat  transfer  character¬ 
istics  of  nanofluids  at  different  application.  Applied  research  in 
nanofluids  which  will  define  their  future  in  the  field  of  heat 
transfer  is  expected  to  grow  faster  in  the  near  future  [191,192], 
The  challenges  of  using  the  nanofluids  are  presented  as  follow: 

11.3.  long  term  stability  of  nanoparticles  dispersion 

Long  term  physical  and  chemical  stability  of  nanofluids  is  an 
important  practical  issue  because  the  nanoparticles  always  form 
aggregates  due  to  very  strong  van  der  Waals  interactions,  so  the 
preparation  of  homogeneous  suspension  remains  a  technical 
challenge.  To  get  stable  nanofluids,  physical  or  chemical  treatment 
have  been  conducted  such  as  (i)  an  addition  of  surfactant; 
(ii)  surface  modification  of  the  suspended  particles;  (iii)  applying 
strong  force  on  the  clusters  of  the  suspended  particles.  Generally, 
long  term  stability  of  nanoparticles  dispersion  is  one  of  the  basic 
requirements  of  nanofluids  applications.  Stability  of  nanofluids 
have  good  corresponding  relationship  with  the  enhancement  of 
thermal  conductivity  where  the  better  dispersion  behavior,  the 
higher  thermal  conductivity  of  nanofluids  [194].  However  the 


dispersion  behavior  of  the  nanoparticles  could  be  influenced  by 
period  of  time.  As  a  result,  thermal  conductivity  of  nanofluids  is 
eventually  affected.  Eastman  et  al.  [195]  reported  that,  thermal 
conductivity  of  ethylene  glycol  based  nanofluids  containing  0.3% 
copper  nanoparticles  is  decreased  with  time.  In  their  study,  the 
thermal  conductivity  of  nanofluids  was  measured  twice;  first  was 
within  2  days  and  second  was  two  months  after  the  preparation.  It 
was  found  that  fresh  nanofluids  exhibited  slightly  higher  thermal 
conductivities  than  nanofluids  that  were  stored  up  to  two  months. 
This  might  due  to  reduced  dispersion  stability  of  nanoparticles 
with  respect  to  time.  Nanoparticles  may  tend  to  agglomerate  when 
kept  for  long  period  of  time.  Lee  and  Mudawar  [196]  compared  the 
A1203  nanofluids  stability  visually  over  time  span.  The  results 
revealed  that  nanofluids  kept  for  30  days  exhibit  some  settlement 
and  concentration  gradient  compared  to  fresh  nanofluids.  It 
indicated  long  term  degradation  in  thermal  performance  of  nano¬ 
fluids  could  be  happened.  Particles  settling  must  be  examined 
carefully  since  it  may  lead  to  clogging  of  coolant  passages.  Choi 
et  al.  [197]  reported  that  the  excess  quantity  of  surfactant  has  a 
harmful  effect  on  viscosity,  thermal  property,  chemical  stability, 
and  thus  it  is  strongly  recommended  to  control  the  addition  of  the 
surfactant  with  great  care.  However,  the  addition  of  surfactant 
would  make  the  particle  surface  coated,  thereby  resulting  in  the 
screening  effect  on  the  heat  transfer  performance  of  nanoparticles. 

33.2.  Increased  pressure  drop  and  pumping  power 

Pressure  drop  development  and  required  pumping  power  during 
the  flow  of  coolant  determines  the  efficiency  of  nanofluids  applica¬ 
tion.  It  is  expected  that  coolants  with  higher  density  and  viscosity 
leads  to  higher  pressure  drop  and  pumping  power.  However,  several 
literatures  showing  significant  increase  of  nanofluids  pressure  drop 
compared  to  base  fluid,  Lee  et  al.  [129]  investigated  viscosity  of 
Al203-water  nanofluids  and  ZnO-ethylene  glycol  nanofluids. 
Results  clearly  show,  viscosity  of  nanofluids  is  higher  than  basefluid. 
Vasuet  al.  [198]  studied  the  thermal  design  of  compact  heat 
exchanger  using  nanofluids.  In  this  study,  it  is  found  that  pressure 
drop  of  4%  A1203-H20  nanofluids  is  almost  double  of  the  basefluid. 
Pantzali  et  al.  [193]  reported  that  there  was  substantial  increase  of 
nanofluids  pressure  drop  and  pumping  power  in  plate  heat 
exchanger.  About  40%  increase  of  pumping  power  was  observed 
for  nanofluids  compared  with  water.  Peng  et  al.  [199]  reported  that 
the  frictional  pressure  drop  of  refrigerant-based  nanofluids  flow 
boiling  inside  the  horizontal  smooth  tube  is  larger  than  that  of  fluid, 
and  increases  with  the  increase  of  the  mass  fraction  of  nanoparti¬ 
cles.  The  maximum  enhancement  of  frictional  pressure  drop  can 
reach  20.8%  under  the  experimental  conditions. 

33.3.  Thermal  performance  of  nanofluids  in  turbulent  flow  and  fully 
developed  region 

The  weakness  of  thermal  performance  of  nanofluids  in  turbulent 
flow  and  fully  developed  region  is  an  important  issue  that  must  be 
addressed  carefully.  Recently,  there  was  inconsistence  of  results 
reported  by  the  researchers.  Kim  et  al.  [200]  revealed  that  no 
convective  heat  transfer  improvement  was  noticed  for  amorphous 
carbonic  nanofluids  in  turbulent  flow  despite  8%  improvement  in 
laminar  flow.  However,  Duangthongsuk  and  Wongwises  [201] 
reported  that  heat  transfer  coefficient  of  Ti02-water  nanofluids  is 
higher  than  basefluid.  This  property  increases  with  the  increase  of 
Reynold  numbers  and  particle  concentrations  ranging  from  0.2%  to 
2%.  Although  the  study  revealed  that  26%  enhancement  can  be 
observed  for  nanofluids  with  1%  of  Ti02  nanoparticles,  it  showed 
opposing  results  at  2.0%  volume  fraction.  The  study  indicated  that 
heat  transfer  coefficient  of  nanofluids  at  this  condition  was  14%  lower 
than  basefluid.  Pantzali  et  al.  [193]  added  substitution  of 
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conventional  coolants  by  nanofluids  seemed  beneficial  for  laminar 
flow  compared  to  turbulent  flow.  Moreover,  Ding  et  al.  [202] 
reported  the  weakness  of  thermal  performance  for  nanofluids  at 
fully  developed  region,  they  found  that  convective  heat  transfer 
coefficient  of  nanofluids  at  low  Reynold  number  has  the  highest 
value  at  the  entrance  length  of  the  tube,  starts  decreasing  with  axial 
distances  and  eventually  reaches  constant  value  in  fully  developed 
region. 


31.4.  Higher  viscosity  and  lower  specific  heat 

The  viscosity  of  nanoparticle-base  fluid  suspensions  increases  in 
accordance  with  increasing  particle  concentration  in  the  suspension. 
So,  the  particle  mass  fraction  cannot  be  increased  unlimitedly 
[129-134].  In  addition,  the  ideal  heat  coolant  fluid  should  possess 
higher  value  of  specific  heat  so  the  fluid  can  exchange  more  heat. 
Previous  studies  revealed  that  the  specific  heat  of  nanofluids  is  lower 
than  basefluid.  Namburu  et  al.  [203]  reported  that  CuO-ethylene 
glycol  nanofluids,  Si02-ethylene  glycol  nanofluids  and  Al203-ethy- 
lene  glycol  nanofluids  exhibit  lower  specific  heat  compared  to 
basefluids. 


31.5.  High  cost  of  nanofluids 

Nanofluids  can  be  produced  by  either  one  step  or  two  steps 
methods.  Both  methods  require  advanced  and  sophisticated  equip¬ 
ments.  This  leads  to  higher  production  cost  of  nanofluids  which  is 
one  main  reason  that  may  hinder  the  application  of  using  nano¬ 
fluids  in  industry. 


11.6.  Difficulties  in  production  process 

The  previous  efforts  to  manufacture  nanofluids  have  often 
employed  either  a  single  step  that  simultaneously  makes  and 
disperses  the  nanoparticles  into  base  fluids,  or  a  two-step 
approach  that  involves  generating  nanoparticles  and  subsequently 
dispersing  them  into  a  base  fluid.  Using  either  of  these  two 
approaches,  nanoparticles  are  inherently  produced  from  processes 
that  involve  reduction  reactions  or  ion  exchange  [190].  Further¬ 
more,  the  base  fluids  contain  other  ions  and  reaction  products  that 
are  difficult  or  impossible  to  separate  from  the  fluids.  Another 
difficulty  encountered  in  nanofluid  manufacture  is  nanoparticles' 
tendency  to  agglomerate  into  larger  particles,  which  limits  the 
benefits  of  the  high  surface  area  nanoparticles.  To  counter  this 
tendency,  particle  dispersion  additives  are  often  added  to  the  base 
fluid  with  the  nanoparticles.  Unfortunately,  this  practice  can 
change  the  surface  properties  of  the  particles,  and  nanofluids 
prepared  in  this  way  may  contain  unacceptable  levels  of  impu¬ 
rities.  Most  studies  to  date  have  been  limited  to  sample  sizes  less 
than  a  few  hundred  milliliters  of  nanofluids.  This  is  problematic 
since  larger  samples  are  needed  to  test  many  properties  of 
nanofluids  and,  in  particular,  to  assess  their  potential  for  use  in 
new  applications  [204].  As  an  undisputed  view,  the  nanofluids 
have  more  advantageous  than  disadvantegous  for  usage  as  cooling 
fluid.  This  calls  for  a  more  intensified  effort  in  the  research  on 
nanofluids.  Research  needs  to  examine  closely  a  variety  of  issues, 
such  as  synthesis,  characterization,  thermophysical  properties, 
heat  and  mass  transport,  modeling.  Hence,  a  multidisciplinary 
approach  comprising  researchers  such  as  thermal  engineers, 
chemical  technologists,  material  scientists,  chemists,  and  physi¬ 
cists  needs  to  be  undertaken.  Only  such  an  approach  can  ensure 
the  future  of  nanofluids  [205], 


12.  Conclusions 

In  this  article,  a  comprehensive  review  of  previous  efforts  is 
presented  for  different  convective  flow  regimes  and  heat  transfer 
through  microtube  (MT)  and  microchannel  (MC).  The  effects  of 
several  parameters  in  geometry,  boundary  conditions,  and  types  of 
fluids  were  extensively  introduced  and  investigated.  In  addition,  an 
extensive  review  for  preparation,  parameters,  mechanisms,  charac¬ 
teristics,  convective  heat  transfer  enhancement,  applications,  safety, 
and  environmental  impacts  for  nanofluids  and  challanges  was  also 
reported.  It  is  found  that  nanofluids  will  solve  many  technological 
problems  and  space  limitation,  and  will  be  used  commercially  such 
as  other  ordinary  coolants  in  the  future.  It  was  observed  from  the 
literatures  that  many  works  have  been  done  using  conventional 
fluids  to  study  the  heat  transfer  and  fluid  flow  characteristics  in  MT 
and  MC.  Thus,  more  work  using  nanofluids  is  needed  to  investigate 
the  heat  transfer  enhancement  in  both  MT  and  MC.  The  outcomes 
from  this  literature  are  summarized  as  follow: 

•  It  has  been  found  that  Si02-EG  nanofluid  has  the  highest  Nusselt 
number,  followed  by  ZnO-EG,  CuO-EG,  A1203-EG,  and  lastly 
pure  EG.  The  Nusselt  number  for  all  cases  increases  with  the 
volume  fraction  but  it  decreases  with  the  rise  in  the  diameter  of 
nanoparticles. 

•  It  was  also  found  that  triangular  MC  has  the  highest  thermal 
efficiency,  the  trapezoidal  MC  was  inferior  to  triangular  but 
superior  to  rectangular. 

•  Greater  enhancement  was  observed  when  nanofluids  were 
used  in  MT  and  in  MC,  as  well  as  other  applications. 

•  Theoretical  and  experimental  research  investigations  are  needed 
to  comprehensively  understand  the  heat  transfer  mechanism  in 
nanofluids. 

•  The  improved  thermal,  rheological,  and  heat  transfer  perfor¬ 
mances  for  nanofluids  has  been  reported  is  some  articles  without 
correlating  these  performances  with  the  specific  applications. 

•  The  production  cost  and  nanofluids  stability  are  the  major  factors 
that  hinder  the  commercialization  of  nanofluids.  Solving  these 
challenges,  it  is  expected  that  nanofluids  can  make  substantial 
impact  as  coolant  in  different  applications 
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